Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


Étude des problèmes aux limites de l'équation de Laplace en configuration cartésiennes trois 
dimensions. 


Typiquement on doit résoudre un problème suivant : 
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La séparation des variables dans les trois coordonnées cartésiennes introduit deux variables de 


séparation a, et où : 


Ӧ?Т (х, y, z) А Ӧ?Т (х, у, z) E OT(xy,z). 
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ах dy dz 


Par principe de superposition ce problème peut se décomposer en 6 sous-problèmes chacun où 
dans une des dimensions х,у ou z, les conditions aux limites sont inhomogènes sur une seule des 
six faces. Commençons donc par le simplifier sur le probléme inhomogène dans la dimension x : 
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Traitons les conditions aux limites dans les dimensions homogènes, à partir d'une solution 
générale développée en série, on remarque que la variable de séparation sur l'équation 
différentielle ordinaire en x est la somme des carrés des variables de séparation sur les 
coordonnées y et z. D'autre part il faut également tenir compte de la présence de la solution de 
valeur propre nulle. Elle apparait lorsque les conditions aux limites homogènes peuvent impliquer 
une indépendance simultanée de la solution en y et z, dans ce cas la solution linéaire A+Bx doit 


également étre retenue, cela équivaut aux valeurs 
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Par contre si l'une ou l'autre des valeurs propres est nulle, alors il s'agit simplement d'une 
réduction du probléme de 3 dimensions vers 2 dimensions, soit déjà envisagés dans le 
développement en série, et ne changeant nullement la forme hyperbolique de la fonction en x. 
Comme on peut le voir en injectant une solution de valeur propre nulle dans les conditions aux 
limites mixtes homogènes, c'est seulement dans le cas où les conditions sont strictement de 
Neumann de part et d'autre que l'on peut trouver une solution non nulle et qui s'avére également 
étre trivialement une constante : 
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Solution A, = B, =0 sauf si =0= Bila; + B71.)+ai B? = 0 B; = B; =0 








0 
Вг аг Bi 
Вг = В; =0> aj B. =a В, = 0 B, = 0 7(2) = A, 
Cela traduit bien le fait que dans се cas la solution пе dépend plus de la variable dans Іа direction 
de laquelle les conditions aux limites sont homogènes, lorsque les conditions aux limites 


inhomogènes sont constantes dans cette direction. Conséquemment le problème est à chaque 
réduit d'une dimension. Pour des conditions aux limites inhomogénes quelconque, le système 


orthogonal de fonctions propres complet devient alors : 
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La composante de valeurs propres simultanées nulles `” n'admet comme solution 


qu'une fonction dépendant de x, de solution triviale T(x)2Aoo*Boo x, et alors le développement en 
série admet la forme : 


пу=+00 nz= *» Cos(A;, ,J)Cos(A, AE 








(AX, + В ох)+ тоу! x( P. „Сов (ил, „+ В, „5їпй( и}, ol 


mets тезе (AY Cos (A2, y) + B}, Sin(Æ;, y) N AZ Cos(4;,z) + BZ Sin; z))x 
F 2 x (A3, „Сози, „х) + Bï Sinh(u; „х)) 
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La condition aux limites dans la direction y donne : 
T(x, У, 2) Е X ous YY OZ. (z) 


К j OY, (y) 3 
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La condition aux limites dans la direction z donne immédiatement par transposition des notations 
deyenz: 
Z (z) = ас A7, Cos(A 2) + Bs Sin( X; 2) 
2 
Equation transcendante А, tq (5) оса; — p, Ds sinc; L)2A, (о; D + aj B; \Cos(% l.) 
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Les conditions inhomogénes en x donnent : 
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D'autre part les normes se calculent ainsi : 
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Il nous reste à calculer la composante de valeur propre nulle dans le cas Во = Вг = By =Й =0 : 
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x-0 
x x x OX (x) x x 
=> Eau (х) = Coo (o; + f х) = = m Сю 
LL 
ау | dz fò (y,z 
x OX, oo (x) x x De x x REY x x | y| fi o 
a, << h X, ou (x) - f 0.2) > Chola; Bi + Bras + 868, 1)- 11 
x-l, Yy z 
1, L 
[dy] 770,2) 
>C 0 o et X, 097 Chola + Bax) 


M ДЕДЙ а ВВ) 


Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


Tout bien réfléchi et avec cet excès de notation qui nous caractérise, on parvient à la solution : 
Y, (y) = a$ A; Cos(A;, y) * By Sin(Z,,y) 
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Prenons maintenant le probléme complémentaire sur la face opposée : 
Ӧ?Т (х, y, z) E Ӧ?Т (х, y, z) А Ӧ?Т (х, у, z) 
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La solution se construit де la même façon, mais еп développant avec les fonctions en x suivantes : 
X „o Œ) = Ау „Сози „(1, — х))+ Bs, ти», A, — х) 
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La composante de valeur propre nulle se calcule ainsi : 
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Ce qui donne la solution à la condition homogène sur la face opposée ajoutée à la première 
solution : 


(lei A, Cos) Віту) 2, tq (а, ага? – prp Jin) = 2, (a? Bt +a? В} Cos) 
Z (z) = аќ А.Соѕ(А.2)+ ВіЅіт(А с) AL tq (у aja! — B; Bi bai L)- AL (n B; + ai B; Cos.) 


Hia = үй, AS > Ron) аби", Cosh(us „ (1, — х))+ BrSinh(u;, „ (L, —x)) 


bs 
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+ ai B 
ou 
yo =Á Is Cos(4),1,) - A, Sin(X! IU „e| + (2, E 2 1 
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ои | 
zo = AL Cos(Æ.L)- А„«{5їп(А* 1.) JE =| Ze | 

Сы + Вх) Саг + 870, х) Z e — o E M ees (^ i »e 

ny=0 nz= y 
T'(x,yz)=isi Bi-B;-Bi-Bi-0 avec А, = T E. ds = 2 
3 
V Cs E E AE OZ 
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La solution compléte du probléme général présenté en tout début de chapitre, est une 
superposition des solutions pour chacun des axes. Et chaque solution sur l'axe se déduit de la 
solution calculée auparavant, par simple permutation des coordonnées, comme ceci : 


T(x,y,z) T" (x,y,z) - T" (x,y,z) - T' (x, у, 2) 

T?” (x,y,z) 2 T'(x y, y х,2) 

T'(x,y,z) 2T' (xz, y,z > х) 

Pour illustrer cette construction et c'est le seul exemple que l'on donnera, prenons un probléme de 


Dirichlet avec 6 faces du cube à des températures différentes. Soit le probléme suivant : 
AT(x,y,z) 0 


aj =a} =a, =a =a} =aj =0 
Po = By = Bo = Вг =Й; = Bí =1 
Т(х,у,2) a = O,2D=T ТО, у,2) ADT 
а= JU. eT 






































Tanda KEDR Tara A CAN 
Les conditions aux limites homogènes impliques les équations transcendantales suivantes : 
2 l 2 | 
0) = Sinay) Z.G)-SmQ2 [ON =. O => 
Equation transcendante &, tq Sin(4,l,)=0— А, ; AA 
у 





nz z nz z 


Equation transcendante À. tq Sin(A71) 202 А = < "1, Cos(Z; 1.) = (-1)* 


2 


Хоа) = Sinh(u3 Mall A Q9 = 5їп(иў, x) 





jo Sin(4;, y) =l, ro ја Sin(A; z) =|, аш ы 
N.T 


лл 2 


2п +1 | + a 
2n 1 

E 0 DH D DH X y 2 

O,ny,nz l,ny,nz 51 п, pair ou п, pair Hny+1,2nz+1 7 | | | 


X — 
Ly nz л 











=> C$ 
L, L 
d 167; e 167" 
Cô, 2ny+1,2n2+1 = ` = С onis F = | x 
л (n, È DOn, + DSinh( илл, m л (2п, + Dn, + DSinh(u,,,. 211.) 
8 ov (r; Sinh( 13... 2nz+1 (I. —х))+ ТЁ$їп(ИЎ ll Sin (2n, + рл у Su (n, + Dr J 
ny-0 nz-0 Qn, + DOn, + DSinh(u,, zech, l 


2 
Аи passage on remarque qu'il ne peut y avoir de solution de valeur propres nulles simultanément 
dans les deux dimensions y,z. 











Ts (x, ўа) 
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Avec la substitution suivante, il vient : 
T'(x,y,z) 2 T' (x y, y > x,z) 


Е mal (ai 
Hä) äi “Л E + ] 


x 2 


TO WC (T? Sinh(u?,..s s, E y)^ T? Sinh{u nv) Sin Qn, + 0л x Sin (2n, T 0л z 
л° nx-0 nz-0 (2n, + DOn, ЗЕ DSinh(u. 3 zuel.) 5 l 


T (x,y,z) 2T' (x z, y,z > x) 


f Zu AU (2n, +1 
Hie äu = T E 5а + ] 


Т? (х, у, z) = а UN (TESI H už о (I. = 2))+ ЛЕДЕН) Su (2n, + Dr < Jo (2n, T Эл у 
Л ^ w-0 ny=0 Qn, + (On, + DSinh(u, 3 zait, І 


J 
D'où la solution totale du problème : 
T(x,y,z) - TG, y,z)+ T” (x, y, z) + T* (x,y,z) 
E nz=+œ (ТИЗ о (7. - x)}+ TSinh(u all Sin (2n, ris j Su (2n, +1)л J " 
ny-0 nz-0 (2л, =t D(2n. + DSinh(u,,. ul.) l l 


› 
16| "EZ (Т Sihh зс (1, us y))+ T? Sinh All Su (2n,+D7 Jd Qn, +1)л J F 
l 





=> T?” (x, y, z) = 

















x 





2 


T(x,y,z) = | + 
ее Z 2 Qn, * DOn, + DSinh(us,, als) 


à N "V (TE Sinh( tees any (1. = :))+ T7 Sinh(t cl si (2n, +1)л Jed Qn, * Dr ) 
l l 
2 


пх=0 пу=0 (2n, + DOn, + DSinh( Hinnale) x у 





х 2 
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Étude des problèmes aux limites de l'équation de Laplace en configuration cylindrique polaire et 
axiale trois dimensions, cas du cylindre complet plein ou creux de hauteur z, cas d'un secteur de 
cylindre plein ou creux de hauteur z 


Typiquement on doit résoudre un probléme suivant pour un cylindre creux : 
O'T(r,O,z) lOT(r,OO, 1 OT(r,O,z) O'T(r,O,z) 
+ + d = 
or? r or r? дө? oz? 
GL. 
T(r,z) fini 
aiT,(r,0,z2)- BIT(r,0,7) _ =-FH(0,23) — ŒT,(r,0,2)+ ETEB a. = f! (0,z) 





0 


a$T,(,0,2)- BiT(r,0,2)) _ =-fi (0) aiT,(r.0,z)+ B;T(r,0,2)) _ = DÉI 
Et pour un cylindre plein : 
O°T(r,0,2) 10T(r,O,z) 1 © T(r,0,2) O°T(r,0,z) 
2 + +> 2 + 2 = 
Ór r Ór r 00 Oz 





0 


C.L. 
T(r,O,z) fini 
a,T.(r,0,2)+ B,T(r.O.z) | = f, (0.2) 


a;T,(r.0,2)- BiT (0,2), — 7:050) — aiT. (0,2) В:Т(,Ө,гу _ = fi (0) 


La séparation des variables dans les trois coordonnées cylindriques introduit également deux 
variables de séparation a» et оз : 


ô*T(r,0,z) , 1 ôT(r,0,z) , 1 &*T(r,0,z) O°T(r,0,2) _ 














0 
ôr? r дг r^ = de? дг? 
=> T(r,0,z) = R(r)O(0)Z(2) 
d^R(r) ldR(r) (a d'or d?Z(2) 
jj) + Per „2 +a, |R(r) 20 SCH + a,0(0) = 0 S +a,Z(z)=0 


Là encore on décompose le problème en plusieurs sous-problèmes par le principe de 
superposition en exhibant sur chacun des problèmes une dimension où les conditions aux limites 
sont homogènes. 


Le Jacobien en coordonnées cylindriques et la métrique sont : 


Ox X 0 
26 0 Cos(Q) -rSin(0) 0 adi аар 
J=det Z L Die del Sin(0) rCos(9) 0|-2 2 5 © -= (Cos? 9. Sin^9)- ғ 
ôr 00 À ^ || ôr ô@ 000 
Qr. 0 1 


ds? = dr? +r? d? + dz? 
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Singularité du système de coordonnées cylindriques 3D et régularité du gradient 


Le Jacobien s'annule en r=0, ce qui implique que /e système de coordonnées cylindrique comporte 
une singularité. Cela se traduit par l'annulation du Jacobien en r=0 et la singularité "apparente" du 
Gradient de toute fonction solution de l'équation de Laplace. La solution doit donc neutraliser la 
singularité apparente du gradient en assurant une valeur bornée à ce dernier soit : 
Grad(T(r,0,z)) = GP S082) I, + оо) I, + ото.) І. (r,0,z)eQ/r=0 
Or r 00 Oz < 
eren m een) у SC 
Or r? 00 Oz 
Singularité du Jacobien des =r =0 condition Grad(T(r,0,z)) reste bornée Vr,0,z 


1 zZ 2 aa 

r? 00 Ce 00 

Cette contrainte est de facto remplie par la contrainte de finitude de la solution: Exemple : les 

solutions admissibles par séparation des variables pour un cylindre complet (voir plus loin) sont du 

type: 

T(r,0,2)= У УІ, Ar) A,Cos(n,0)+ B!,Sin(n,O)\ A; Cos(4;,z) + Bi Sin(A; 2)) > R(r)=1, (Ar) 
ng-0 nz=0 

Le respect de la finitude de la solution pour les valeurs r=0 suffit pour assurer la condition de 

finitude du gradient soit : 


? (r An) zo Y += 2 „ү 
r 2 DIE Am, AU, 2 


r k=0 


2ng 2k 
At +оо 1 Д2 
= pot) ^nz uz > Cste si n, 21 
" 2 £xT(k n, 41K 2 











Grad(T(r,0,z)) = | 


oc r? Cste 








> Cste =| 





г=0 














1 (0Т7(,0, 2) 
? 00 


E 
Dans le cas d'un secteur cylindrique 








Ha > „9 =-7- Өл 2А, 2021; > Cste si, 21 
Д r 00 
Lorsque l'angle d'ouverture du secteur cylindrique est plat (a=n), les valeurs propres sont entières. 
Il reste le cas où l'angle d'ouverture du disque est plus grand que т. Dans ce cas, on fait encore 
appel au raisonnement impliquant la construction de la solution sur le probléme complémentaire 
du secteur au cylindre entier. Appelons problème n°1, le probléme sur le secteur d'angle obtus. Il y 
a existence et unicité de la solution du probléme de Dirichlet intérieur construite par les fonctions 
de l'angle obtus. Cette solution s'étend mathématiquement à tout le disque, donc sur les frontiéres 
du secteur d'angle aigu. Cela induit donc un probléme aux limites de Dirichlet n°2 sur le secteur 
d'angle aigu. Ce dernier se construit également à l'aide des fonctions propres qui posséde la 
régularité du gradient requise sur le système de fonctions et valeurs propres 


2k 
Ka 
24,6 + CX E 
= prof À. 2 
r=0 


пл 1 (2022) 
2  T(k+n,+1)k! 





B=2n-0, А, = I m лт,8 <n = А„> п >1. Et par existence et unicité de la solution induise 


le respect de la régularité du gradient de la solution dans son domaine complémentaire qui doit 
être identique à la solution du problème n°1. Ce type de raisonnement s'étend à d'autres types de 
configuration avec des domaines complémentaires. 
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Orthogonalité bidimensionnelle des fonctions propres sur un problème homogène cylindrique 
dans les deux dimensions radiale r et angulaire 0 


On se propose de montrer l'orthogonalité des fonctions propres : 
T(r,0,z2) -U(r,O)Z(z) U(r,0)=R,,, (r)O,, (Ө) 












































2 2 
Avant la séparation des variables r,0 — 080 + н А nee) = (4, FU(r,0) 2 20) = (4, P ZŒ) 
r Or Ог r 00 Oz R 
Deux solutions 
1 д 2,0, п (7,0) 1 OU, п, (7,0) 
1 — + 1—0 +(4, FU ,0)=0 U ,0 
0 5 A - z =>. S an O x U, uL (0,0) 
1 Ô dr, т (7,0) 1 О, т (7,0) 
2 = + D +(4 FU .0)=0 x U ,0 
o - é( < = = E EU snb) sa (+8) 
1 a{0rU, „ (r,0) 1 OU, „ (r,0) 
= c + ct +, JU ,0) |U ,0)- 
t é( Or r? дө? ( т; n, ng (r ) m, mg (r ) 
1 8 ОО m (r,0) 1 eU, m (0) 
= + 1770 +414, FU ,0) U, „ (7,0) 0 
É al - Ае EE 
Or ‚| Or ,0 
кок, nn (50) U. të 1 ô ( OrU „©“ IA, U, um 
r Or Or Sch r òr Ог ® 
1 ÔU, n (0), 1 ӨШ „ОЪ 25 
+ р? 20? DU. ‚тө (r, 0) r? 20! U ‚пө (r ,0)= (a - ES cae (SEU m, m (r,0) 
1 д п, ‚п, (r, 0) OU, JH (r,0) 
=> (3 2 r, 0) — —— = ol 
(3) EE zs U m. m, r0) 3; Un n (Fs ) 
1 д n. yt (r, Ө OU n т, (r,0) 
+ 2 ,9--_*”2"” "y 9 - V OU n „(т.б 
r? SE дө ан (r ) 30 p. ‚лө (r, ) = ((à if ( 22 uM (r ) Е SE ) 
% n 9, oU . (r,0 ôU. (r,0 S 
Intégrons  [dOfdrr(3)= | do| r Un n (ON y „ (r,0)- 206) р „ (1.0) || + 
9 n 9 дг КЫ дг Ut n 
ro 0 QU ‚Ө 9% 9, Dy 
jar P E sa 00) dua Go ne J U, (r, d cl, P -( PI d6 dr rU, „ (r.0)U„ „,(r.0) 
9 9 n 
U, п (r, 0) OU, п (r,0) é 
C.L.— OU OO у U m. m r0) - —— оди) =0 
И ôr эле òr S , 
Toi (7,0) OU, „ (r,0) " 
C.L. | _ U n m (5,0) - —————U, ,(r,0)| 20 
| 00 dii 00 „се z 





(4, Y -(a„ f j aof drrU, „ (OU, . (r,9) =0 


9 n 


en (rU. m (r,0)=0 Vm,,n,/m,#n, 


5 [dol art 
9 


ñ 
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Le cas particulier d'une valeur propre nulle, permet de démontrer l'orthogonalité bidimensionnelle 
de la solution aux valeurs propres nulles du problème aux limites avec les fonctions propres de 
valeurs propres non nulles. Et lorsque les deux valeurs propres sont nulles alors le problème aux 
limites se ramène à une probléme réduit d'une dimension, le cas radial (rz). L'orthogonalité des 
solutions est bien établie également dans ce cas (soit par la théorie de Sturm-Liouville, soit en 
calcul directe). 


L'orthogonalité bidimensionnelle permet de construire facilement les solutions aux limites par 
double intégration sur les deux variables, avec un poids radiale (r) bien établis, sous une forme du 
type : 


T(rüm-» Y A, al (r»O)Z„ (2) 


n, ng 





9, 2 
[do] drr rte, OÙ, „,(r.0) 
cp A, = 





n, ng 


[ao drr(U, ,,(r,0)) 
8 


п 


Urn Ç> O)" = [49] arr, 09) 
9 n 
[aol arr ro ou. (r,0) 
>A = 


n, ng 





Un, C0) 
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Le cylindre plein complet : 





Compte tenu de la géométrie du solide, il ne peut y avoir de problèmes aux limites inhomogènes 
angulaires car les limites du domaine ne comporte pas d'iso-surface à angle constant. Les seules 
iso-surfaces sont à rayon constant et à hauteur constante. Par conséquent, les seules problémes 
aux limites entrent dans les deux catégories suivantes : 


- le probléme homogène axiale (en z) , il correspond donc à un probléme inhomogène radial 
- le probléme homogène radiale (en r) , il correspond donc à un probléme inhomogène axial 


Problème homogène dans la dimension axiale z : inhomogéne radial r 


Commengons par le probléme suivant sur le cylindre plein et complet oü les conditions sont 
homogénes sur l'axe z : 
AT(r,0,z)=0 T(r,0,z) fini 


CL. a,T(r,0,z)+ B,T(r,8,2). = f.(8,z) 
aiT,(r,0,z)- PiT(r,0,2)) _ —0  aiT,(r,9,z)+ B;T(r,0,z)) _._ =0 


Les équations séparées de base sont alors : 


d'R(r) | 1 dR() Jr re | ко)=0 LEO A rem 20 Tan 








dr? r dr г? ES 


Selon les principes de la séparation des variables, la solution se développe en partie en série, sous 
la forme : 


"et nes (49, Cos(n,0)+ B? Sin(n,0)J A; Cos(4;,z) + BZ Sin(Æ, 2))x 
T(r,9, z) = Kä к, Gr) 


no=0 п2=0 ^ (4 „1, (Qr) + В, 
Condition de finitude de T(r,0,z) => B, „=0 


ng ,nz 
ng nz 


ng =+00 nz=+0 


T(r.0,2)= У SL (A2,r)(A2,Cos(n,09)+ B? „Sin(n,0) AZ Cos(A?.z)+ B; Sin(A;,z)) 

ng-0 nz=0 
à laquelle il faut ajouter le développement lié à la valeur propre nulle sur la composante axiale 7. 
selon un raisonnement analogue à celui réalisé dans le cas cartésien, le terme axial de valeur 
propre nulle admet une solution non nulle uniquement lorsque les conditions aux limites 
homogènes sont de Neumann de part et d'autre des sections du tube. Et de plus la composante 
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axiale admet uniquement la solution constante. 
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1⁄2. =0 A e í А , : A 
Comme "e , la partie radiale n'est plus une fonction de Bessel uniformément nulle mais la 


solution du problème de Laplace bidimensionnelle polaire, soit en ajoutant à la série précédente 
cela donne : 


T(r,0,2)= У У, Qr X5 Cos(n0) + B^, Sin(n,8) V A; Cos (47.2) + В Sin(Æ.z))+ T,(r,0) 
ng-0 nz=0 
T,(r.0) = Y (Ayr? + Bir"? \A2,Cos(n,0) + B°,Sin(n,0)) 
по=0 
Condition de finitude de T, (r,0) => B,, = 0 


Dg =+9%0 nz=+0 


T(r,0,2)= У Y, (2.04, Cos(n,0)- В°,5іп(п,0) (A,Cos 45, 2) + B} Sin(Æ.2))+ T,(r,0) 


ng-0 nz=0 


Б: eg (49,Cos(n,0) + B°,Sin(n,0)) 


ng-0 


Les valeurs propres et fonctions propres dans la dimension axiale (z) sont les suivantes : 
Z, (z) = оф А. Cos (45,2) + B; Sins, 2) 


nz 


Equation transcendante А. 14 (a. Jojo; - B; B; sina.) = 2: (a; B; + a B; Cos(A5..) 
Valeur propre nulle Z, (z) «1 
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Pour ce qui est de la dimension angulaire, les valeurs propres sont toutes entiéres et se 
développent suivant une décomposition par symétrie de la condition aux limites sur la surface du 
cylindre, le cas a déjà été étudié à deux dimensions. Voici le Tableau synoptique des fonctions 
propres suivant la symétrie des conditions aux limites. Cette décomposition par symétrie polaire 
s'applique aux deux développement : 

















Propriétés de symétrie Symbole de|Fonctions propres utilisées 
symétrie dans le cas d'un disque 

complet en angle creux ou 
plein 

f.(0,2)- f(x -0,z) et f.(0,z)= f,(-0,z) |Y+,X+ Cos(2n0) 

f.(8,z)= f(x -0,z) et f.(0,z)-—f,(-0,z) |Ү+,Х- Sin((2n +1)Ө) 

f.(0,2)=-f.(m-0,z) et f.(0,z)= f,(—0,z) |Y-,X+ Cos((2n +1)0) 

f,((8,z2)-—f,(m-0,z) et f.(0,z)=-f.(—-0,z|Y.,X- Sin(2n0) 

















f, (8,2) = fiy x (6,2) fiy x (0,2) fy x:(0,2)+ foy- x (0,2), sachant que f, (8,2) = f. (2x +0,z) 
Ўл 6,2) = ARE UR A E) d 
J.(9,2)+ f. (t -0,z)- f. (70.2) - f, (9 — 7,2) 
4 
J,(0,z)— f, (m 0,2) + f. (70.2) - f, (9 — 7,2) 
4 


Ууу 0,0) L9 f. (x —0,z) Al 0,z)+ f.(0-x,z) 


A partir de cette décomposition, on calcule les coefficients du développement en série ainsi : 
Symétrie Ү+,Х+= C.L. f y, x, (9,2) 


Z (2) = ас A;,Cos(A, 2) + B; Sin(A7 2) 


nz 











Jy: (0, z) а 





Jr ze (9, z) = 

















zo| = Ir 2Соз(А1.) + B; Sin(Æ.L)) [ ыш Jr p dE 
a;(x.) о; 
ои 


z-o As ? Соѕ(2:1,)- А,а5т(:1.)) |1 H Zen? KE UE 


Ng =+% nz=+œ0 


(+д„) T(r.0,2)= У Y. Am nln, Sir) CosQn,0)Z, (z) 


ng-0 nz-0 








@„, (9) =Cos(2n,0) | 





Í do | dz f, y. x (0,2)C08(2n,0)Z,„ (z) 
Di + дьо |; 


2л L л L 
De plus |40 | dz f. y, x. (6, z)Cos(2n,8)Z,, (2) = 2| a0 | dz f. y, x.(0,2)Cos(2n,0)Z,. (z) 
0 0 0 0 





2ng nz 


Z CN GAL ' ALI) + BL, AL) 
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Ainsi que les termes issues de la valeur propre nulle axiale : 
Symétrie Ү+, Х+= СІ. ў. у. x. (9,2) 


ng =+0 п2=+00 


(+д„) T(r.09,23)= У Y A, о" Cos(2n,0) 


ng-0 nz=0 


o. (9) =Cos(2n,0) | 





2л L 


л L 
1, [a0 dz f, у. x, (0,2)Cos(2n,0) j^ [a0 dz f... (0,2) Cos2n,0) 





EP = 


0 0 = 0 0 

a 1,1" (1+8, Ka, 2n,* BA) т 1/"(+0,.).(a,2r,+B,1.) 
2л 1, л LU 

De plus |40 | dz f, y, .(0,2)Cos(2n,0)= 2 | d0 | dz f. y, ,.(0,2)Cos(2n,6) 
0 0 0 0 


A partir de cette expression on peut étendre à tous les autres cas de symétrie et donner la solution 
du probléme en conditions aux limites axiales homogénes comme ceci : 


J,(0,2) = feas x (0,2) * fuy. x -(0,2)+ y. x. (9.2)+ foy- x (0.z), sachant que J, (09, z) m f. 27 +0,z) 
/.(0,2)+ f. (x -0,z) * f,(-0,z)+ f. (0—z,z) 




















frx (0,2)= S 
о 400 (0-06) LCS) 002) 
faa у с _—_ 
ае) L 0,2)+ f,(0=7,z) 


Z,.(2) = а À, Cos(A5,z) + E Sinti) 
х. q (каа om sina) = az (ar Bi + a B cosi) 
L 
б = | dO Cos(2n,0) | dz Zo (Ð) frx (02) cli, = Í аө Sin((2n, +1)0)[ dz Z„ (2)... 0.2) 
0 0 


2n,+1,nz 


T 


л 1, 
ie = | d0 Cos((2n, +1) j а Z„ (D) fr x (0,2) су. = | dO Sin(2n,0)] dz Z,.(2)f.y x. (0,2) 
0 0 


Cong +1,nz 2ng ,nz 
0 
L 
g H е | do Cos dd кед ду ups j аө Sin((2n, +1)0)[ az f. 5.» (0,2) 
0 0 


I. л L 
= j аө Cos(Qn, +1)0)[ dz fs x,(0,2) сш = |40 Sin(2n,0)[ dz f. y x,(0,2) 
0 0 0 0 


Cas pj = В, =0 


Z,,(2) = Cos(#,2) = Е d 


2 





nz T 


X. dq Sin(æÆL)=0>%, = | Z (3) == 





et 








nz z 
2 
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Ce qui donne les développements en série suivant : 
Si B, f, 2-0 





wer d'Zeit Үгү" 
ng=0 ( T до Xe, 2n, + BI ) l, 
ES ME 
T (r 0 z) = 21, ng-0 (a, (2n, +1)+ B,1,) 1, 
zhN 9V > T ү ng =+0 с Cos((2n, 4: 1) y 2ng4l 
£X («„(2л„+1)+ 8,1.) (1, 
esc EE | r | 


| ng=0 (a, 2n + B,1,) L, 








~x 
N 

















Y+,X+ 
Cn, nz 1 


ng =+% nz=+0 , 2ng 


EEN É E | 


2 


2, H ( 2 | OZ AX nz T 
1 T до A „Anlon 1 l, + Blan, 1 l, 


2 


бым e Ju, + Okt 











z 








+Y У z + 
Lu a Al gb E l, + BL +1 с 1, 
0 l, 0 l, 
"ei 
Tt, 
ng =+% nz=+œ ER жын | - а oan + | T я | 
Z 2 + 











+ 
2 Z 2 | NZ T Nz T 
QA, l l, + BL, +1 l l, 


XI (5E нал 7 


ng =+œ0 nz=+œ ` Zo MZ ` 2Ng 
2 2 


+ 
2, Z A Nz m Nz T 
|е, | l 1 + P.L. | 1 l, J 


Sinon — 

















Peg 


2ng,nz * 2ng (4,.r)Cos(2n,0)Z,. (2) 
= e (1+0, JO GAL OL BE, AL) 
Cere «Шы CAE 
ng-0 nz=0 Znz (z) | la ALL, +1 ' (LA) T Ds. +1 (Al, )) 
+ "R Cu ER ES +1 (А r)Cos((2n, + 18)Z,. (z) 
ng=0 п2=0 2,3 (z) f ERNE | (LA) ad DI GA )) 
Hg too nz o Cia ns ER (A, r)Sin(2n,0)Z, (z) 


+ 2 
ng=0 п2=0 Z. (z) (0,1. | (GA) E Blan, (42.1, )) 


Пу=+%® nz=+00 








аа 
Л 





Ү-,Х- [ 
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Problème homogène dans la dimension radiale r : inhomogéne axial z 


Passons maintenant au problème homogène dans la dimension radiale : 
AT(r,0,z)=0 T(r,0,z) fini 


CL. a,T.(r,0,z)+ B,T(r,0,z) _ =0 
aiT,(r,0,z)- BéT(r,0,2) =-fé(r,0)  afTi(r,0,2)+ B;T(r,0,2)| = f;(r,0) 
Il se décompose en deux sous-problémes, le premier sous-probléme présente les conditions aux 


limites : 
a,T,(r,0,2)+ B,T(r,O,z| =0 


aiT,(r,0,z)- BéT(r,0,7) | —0  afT.(r,0,z)+ B;T(r,0,z)| = f(r.0) 


Et le second sous-problème: 
a, T:(r,0,z)+ B,T(r,0,z) _ =0 


aiT,(r,0,z)- BéT(r,0,2) =-fé(r.0)  aiT,(r,0,z)+ ВгТ(",Ө,2) =0 


Commençons раг le premier sous-probléme : 
a,T,(r,0,2)+ B,T(r,O,z) =0 


%Т,(т,Ө,®)— BiT (r,0,z) =0  a;T,(r,0,z)+ BFT(r,0,2)| _ = f(r,0) 


Les équations séparées de base sont alors : 


d^R(r) 1dR(r) |@„ („ү _. POO) 3 a. d’Z(2) („ү _ 
л 1 6.) оо E + (п, Ÿ @(8) =0 -(x f zo 


La solution se développe sous la forme d'une première série à laquelle il faudra ensuite rajouter les 
termes de valeur propre nulle radiale : 

` ei n+ (48, Cos(n,0) + B°,Sin(n0)Ÿ A: Cosh(X,,., 
T(r,0,z)= | 
ng-0 nr=0 x (4; Jg J n, (4; r) + B, Y, (Ал, r)) 


nr ,ng n, Ng ^ ng 


Condition de finitude de T(r,0,z) = В, , =0 











„2)+ В; Sinh(A,., z))x 


pr 


ng =+0 nr=+00 


T(rG,z)- У VJ, Qr Mo Cos(n) + BlSin(n,0) (A; Cosh( 4 


r nr ,n 


.z) + ВЕ Sinh(X,, „ z)) 


А 
ng-0 nr=0 
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Termes de valeurs propres nulle 
Les termes de valeur propre nulle appréhendent les solutions possibles du problème aux limites 


lorsque les valeurs propres définies plus haut sont nulles, soit quand la solution ne dépend plus 
d'une ои de plusieurs des dimensions aux conditions aux limites homogènes. On peut dire que le 
problème dégénère en problème de dimension inférieur. Ces termes représentent également la 
part de la solution de dépendance dimensionnelle réduite. Lorsque c'est uniquement la valeur 
propre radiale qui est nulle, Abas = 0' la solution radiale est de la forme, 





d'R() | 1 dR() _%) 
dr? r dr 
Condition a,T.(r,0,z)* B,T(r, 0,2) SCH 


R(r)=0 R(r9A,r"-B,r'  finitude— B, =0 


> А, o ad" + p," J- 0 saufsi п, =0е B, -0— R(r)- A, 


ce qui confirme qu'il doit y avoir simultanéité de valeur propre nulle radiale et angulaire. 


: . &.=0 et n,=0 
Les termes de valeur propre nulle radiale et angulaire ^” 9 sont de la forme 


R(r)=A+B Log(r). Compte tenu de la condition de finitude c'est une constante, qui n'est pas nulle 
uniquement dans le cas oü les conditions aux limites homogènes radiales sont de Neumann : 
a =1 et f,20— T.(r.0,z) = 0. Cela conduit évidemment à des termes de série indépendant 


de r et comme dans la séparation de variable cela conduit aux trois équations différentielles 


ordinaires suivantes : ZRC) , 1 dR(r) _ 0:8) аӨ(Ө) | 0 dZ) 
Уз 


= 0 H 
dr? r dr аө? 


еі 





On en conclut que le terme global de valeur propre nulle est linéaire enz: 
T,(r,0,z)= Aj + Bj z 


ajT,(r,0,z)- BiT(r,0,z)| _ =0=> ai Bj = B; 4; 
T,(r,0,z) = Aila; + Biz) 


Dans ces conditions la forme du développement est bien la suivante dans le cas particulier des 
conditions homogène de Neumann : 


=+00 nr =+œ0 


T(r,0,2)= А (ai + Biz) у V I, Cur Ad, Cos(n0) + Вё„5їп(л„ө)\ Aï Совй(3, „ z) + Bi 8їпһ (А! 2)) 


ng-0 nr=0 


A, tq Ju (2, nr ails )= 0 


nr ng 
Dans le cas général de conditions homogènes mixtes, l'indice angulaire est entier de part la 
condition de périodicité 2л imposée à la solution. Et l'indice radiale est déterminé раг une 
équation transcendantale obtenu en appliquant la condition homogène sur la surface du cylindre, 
comme ceci: 


(r) =J„, (4; 0 


Ran wr) a,T,(r,0,2)+ B,T(r,0,2) _ =] 50Sa Annn ‚(бл L) + BI n, Au 1.) = 


As tq a, À, nr,ng J n, De ng 1.) + B,J по (4, p ) FE 0 
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Pour ce qui est du développement angulaire, on se conforme à la décomposition de la fonction 
limite suivant les 4 types de symétrie angulaire établies : 





Propriétés de symétrie Symbole de|Fonctions propres utilisées 
symétrie dans le cas d'un disque 
complet en angle creux ou 








plein 
(7,0) = Р (,л-0) et (7,0) = fi(r,-O) |Ү+,Х+ Cos(2n0) 
f (r,0)=fF(r,n-0) et f'(r,0)=-f;(r,-0) |Y+,X- Sin((2n +1)8) 





f(,0)=-ff(r,x-0) et (0,0) = 30,0) |Ү-,Х+ Cos((2n +1)0) 





f(r,0)=-f;(r,x-0) et } (т,0)=-—/ Di |Ү-,Х- Sin(2n0) 














Prenons comme exemple la composante de symétrie Ү+,Х+, il vient : 
Rs. (r) = Jy, (Aron, r) e, (0) = Cos(2n,0) > Ryo, '(r) E oos: 2ng (EM ,2ng r) | 


Del 





0%7,(7,0,2)- BéT(r,0,2) fe Z,,, (X) = a, 


nr,2ng 


Cosh(A,, 2n, 2) + Bo Sinh(A,, an z) 


lnr 


TOODAB TEED = сог, | dO fj, х,т,Ө)Сов(2п„бӨ) = 2 | dO 52, x, (r,0)Cos(2n;0) 


1, 
2 | dr A ant) | dO fi. x. (r.9)Cos(2n,0) 


Y+X+ _ 
=> C rano = 





B ta; (| le. ө) (4 i nr an COSh(A n. an, L Ye a; B; += Bo )+ Sinh (А. ,, | Ж a V (2, nr m + Bo By ) 


Р ' 2 
М | MM ш | + Boso (L, y pe Qn, ) 
A 2n; 1, sl » ,2ng ) 











2 
Ron, Gi = 








2 
2 
М BEE an, l. Wuel ©) 





2 
Roos Gi = 


Pour le problème symétrique : 


2л л 
соТ, (r,0,2)- BiT(r.,0,2)) = — fas (70) | dO fiy. х. (r.0)Cos(2n,0) = 2 | dO fiy, x.(r,0)Cos(2n,0) 








2 


aiT.(r,0,2)+ B;T(r,0,z)) | =0=> Z, „(х) =j Ar, „Соз, au, (l, — z) + Br Sinh au, (l, - z)) 


L,nr 1 ^"nr,2ng 


2 | des Jour) | dO fs, x+ (r,0)Cos(2n,0) 
Y+,X+ 
=» O,nr,2n0 ` 





К, (Ө, glo. Соза, з. сг таги) Sinha, „ао (Ana) +8487) 
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Pour ce qui est de la valeur propre nulle, liée uniquement à une condition homogéne de Neumann, 
rappelons une propriété des fonctions propres radiales de valeur propres non nulles lorsque le 


probléme ne posséde pas de dépendance angulaire : 
2 





1, 
2n 20 RO) Л) @ JU GLD) T 0,170. [RO =f drr =” 
0 
А, "oi (A. ol, ) D r 
jo r J (Auro) = "GF =0 > orthogonalité de = R, (r) Jj (Ar) avec R, (r) =1 
Ао 
L, 2л 
Z(z)- А (ač + Biz) [arr [a0 f;(r,0) 





ai Z'(z)+ BrZG), , = 7709.0) > Aui i Pied + E=" 


127 
farr fao j; (r,0) 
(Déeg + Bros + BEBE) 


De plus {40 fë (rO) = [de (rs. tege, (r.0))= 2| dO (75.0.0947, x, (1,0) #0 
0 0 


0 





Aq; = 


2л л 2л л 

[40 070) =J dO fié. x. (r.0)+ [dO угу, ,.0,0)=2[40 угу, x, (r,9) # O 
0 0 л 0 

2л л 2л л 

[40 fé. 0,0) = |40 E x. 0,0) | dO fj, x,(r,0)=2 | dO fj, ,,(r,0) # O 
0 0 


" fane rer 8)- fé (r.-0)- 7209,0 -л)) 
En effet [dB fis, s (r,0)=" 





4 


" 1 J Rd) = fao Ji (r.0) + fao PA COR fao veel =0 


2л 


fao (f (r,0)- f (ra -0)- f#(r.=0)+ fr,- eil 


2л 
De même [40 fix 0,0) = 5 4 5e 
0 





L л 
EEN 
Aj = — E 
l nla B; + Bros + B; Bil.) 


Z(2)= Aj a; BH. ll far r fao fi(r,0) 
De méme aj Z'(z)- Bize ; =- fo (0,0); > -Aio (ai B; +B; aj + Sei => : 
ai Z'(2) + ize), = 








l'n 


r 


2f drr jet foy xs Fs 0) foy- x. s 6) 


L, рут + Pí [774 + B; B; Ï L.) 





> Áo = — 
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Il vient par superposition des 4 composantes de symétrie, ainsi que des 4 autres composantes pour 


la solution du deuxiéme probléme, la solution compléte du probléme : 
AT(r,0,z)=0 T(r,0,z) fini 


СІ. a,T,(r,0,z)+ B,T(r,0,z)., =0 
aiT,(r,0,z)- BiT(r.0, 2) __, = — fo 078) «їТ,(т,Ө,®) + BT(r,0,2) _, = } (7,0) 


Les premiers coefficients, fonctions propres et valeurs propres: 
Ans (r) = J an, (4, r) À, tq A Aar 2n Jy, Ces 1.) + B.J an, (A, on, l.) "S 0 


nr,2ng nr,2ng т 
Кон) = Jana (А, za fl ЖЕ tq A Amr ons Sa et Le. ans be) + BJ a onu) = 0 
@,*"*(0)=Cos(2n,0) 07^" (9) = Sin((2n, 10) Ө; (0) = Cos((2n, -10) 07. (0) = Sin(2n,0) 
Zar ng (z) = А, Cosh(X, 2) + Po Sin h( Mn 21, z) 


nr,2ng nr,2ng 
Zi od (z) = «А. on, 4 COSh(A,, 59, 12) + Po Sinh(A,,,, 12) 
Zo 2n À) = AFA „Созу, 2n, (1, =z2))+ Bi Sinh, (L. — z)) 
Zr ang (2) = ОТА an, СОЗА, a, +1 (I. m г))+ BESI HA, à (I. E z)) 


Jo (r,0) = foy x. O) + Ja, e (050) * foy x, (550) foy y (70), sachant que Jo (0,0) = Jo (r,27 +0) 
fo (0,0) fy (rx -0)+ fy (r,-0) f, (r,0— n) 














foy. x, (0,0) = | 

у: eo- 0 т -0)- fi (n-0)- fé 0,0 ==) 
0oYs,x- M» | 

Jo (r,0)= fo(r,0)- fi (r.a —0) fy (r,-0)- f$ (r,O =w) 
0,Y-,X+ (> 7 

Л, (r,0)= Join, тее) fg 0-0) 9 fç (rO m) 
0o Y-,x- M» S 


s (r,0) — ОВУ Зе) 0S8) fiy-x-(r.0), sachant que Jí (7,0) = Jí (r,27 +0) 
J (r.0)+ fé'(r.,a -0)+ f;(r,-0)+ f/(r,O—n) 











Дух» (r,0) = 1 

f: (бу fj (10,0) (0,0-0) – f;(r,-0)- f;(r,0-7) 
LY+,X-\ ? 4 

f: бру fí (,89)- fr (r.a -0)+ f;(r,-0)- f;(r,0-7x) 
LY-,X-N ? 4 


Л (1.0) - fr (rz -0)- fí (1-0) + fí (1,0 =w) 
4 





Te? (r,0) = 


Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


Avec les coefficients suivants du développement en série : 


l; 
[drr J ang Uy an, d 40 fj, x. (r>9)Cos(2n,0) 
0 


























Y+,X+ _ 
L,nr,2ng 
бна] ЖОО к ï=. re) 
к SE Eus) 
1, л 
[arr Dont e 40 firsa- (r,0)Sin((2n, +F 10) 
Ce wart ° { 
2 
о U жЕ DD т гё | 
l, 4 ^ m. 
1, т 
[arr ГӨ, PER à 40 Jy (r, Ө)Соз((2п, $ 10) 
Cin eum = : : 
2 
E Ud | - Sch, 
l '( nr ssl 
[aer J ang (4, 21, nf dO fj, x:(7,0)Sin(2n,0) 
Cis e = ` i 
2 
«ыз es ы] - + 
а Ly Da ie? | dO fiy. x. (r,0)CosQn,0) 
Cint =f 
ISCH IL ang aD Peas CE A5) 
i; т 
[drr AU, dO für, x-(r.0)Sin((2n, +1)0) 
CES E= ` ` 
2 
(J Jua Arr, any ihr 1)} + PN ЖО Е Sch, 
Оаа 
1, т 
[arr Jo CE || 40 Jorasa (r, Ө)Соз((2п, * 10) 
С issu = i ° 
2 
aUas al] ааа) ты ж 
l, ab nr East 
[drr J ang (А, 2л, nf dO foy:,x:(,0)Sin(2n,0) 
Cosa = : 


Le nr,2ng4l 


o UT eos) 2.) 
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La solution développée en série pour la partie radiale homogène donne, en notant : 


T(r,,z)-» T.(r,0,z). En commençant par le cas de la condition homogène radiale de Neumann : 
Si p,-0 


1 
d по E 
" Jett, LX B ai Bs) SinhA Lear (Ао, + ВВ) 
: 1 


Aini = 
(s, at, Lo BE + aF Ps) SinhA au ala Nasa; osa) + Bs ) 








1, Е 1, л 
2| arr [ав (5. x 0.0) 9 f. 0.0) 2| drr [ав (fé, x. 07.0) fos x, (7,0) 
Z- BEN 0 0 Zi ы 0 0 
Aur L alai Bi + Bias + BB) Lo шер; + Brat + AB) 
Ros, (r) = J an, (Aon, r) d con tq Jan (As 2,1.) -0 
Rosa) = Jua Sen uf Ап tq J angar оп) =0 
— 
T (r,0,z) = А (оф + Biz) AZ (aœ; + Br, -z))+ 


Ng =+0 nr=+œ0 


У z ЖЯ Р Es J an, (o r)Cos(2n$0)Z, „ол, (z) + 


ng-0 nr=0 








ng-—too nr=+0 


Y+,X- e 
+ Kä 2. SEN GË 2ng4l (GO i )Sin((2n, + LEE (2) F 


ng-0 nr=0 


ng =+% nr=+0 


Y-,X 
+ 25 > An и ол ud (Ar, an, ar) Cos((2n, + 18)Z, ana (2) + 


ng-0 nr=0 


ng =+00 nr=+00 


+ > 25 A5 Co Jan, ss r)Sin(2n,0)Z, m. an, (z) + 


4 ng=0 nr=0 





2 ng =+œ 
Бый nr=+0 
zl е 


r r Ү+,Х r 
T > >> 2 [a MUR J an, (os, r)CosQn;0)Z, mang (z) + 


ng-0 nr=0 


ng =+% nr=+0 


+ 25 У ЕКИ Сом ае (s, an) Sin((2n, + D8)Z mangi Z) + 


ng-0 nr=0 


ng =+% nr=+0 


Y-,X 
+ Y Kä асо (2,2, ar) Cos((2n, + PA (z)* 


ng-0 nr=0 





nr,2ng "Dm, Ange" 2ng Ne KT EI CHEN (2) 


+ >. 2,4, CEA J 
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Et le cas des autres types de conditions homogènes radiales : 





Sinon 
r 1 
nr,2ng 
À. 2ng Cosh(A,,. 2m! LY aj Bi + a, í Po )+ Sinh( A., 2ng L.) e oai À (2 nr m. T Bs Bí 
1 


Fr — 
nr 2ng+l — | Е 


A, an à COSh (An al. Xi B? + a7 Bi )+ то, deier (А, а) + Bébé) 


Y+,X 
У A. IN А 2ng (A. 2ng r)Cos(2nj0)Z, „зл, (z) * 


ng-0 nr=0 





по =+0 nr=+0 


> А, ,2ng as = Jua (Aar ax an) Sin((2n, + 1017 +1 (2) E 


ng-0 nr=0 


по =+0 nr=+0 


Y-,X 
Kä А, ‚2лә+1 TM ad ona (as, ur) Cos((2n, + DONZ был (2) + 


ng-0 пг=0 
"n nr--oo 
Y-,x- e 
х 4, ,2ng C, nr,2ng Ju (A. 2n ,DSinQnj8)z,, 2n, (z) + 


ng-0 nr=0 


T, (r,0,z)— 





ng =+ nr=+œ0 


4 j * 2 D А, ,2ng Lo Jy, (A. 2n ,')Cos(2n$0)Z,, os, (z) * 


ng-0 nr=0 


по =+0 nr=+0 


+ У У А, ,2ng TO Jua (А„. 2ng „)5їп((2л, ate 6) on (z) + 


ng-0 nr=0 


ng =+9%0 nr=+0 


+ Kä Kä Any E ee 2лу+1 (4, 2ng „r)Cos((2n, + 1)0)Zo mangi Z) + 


ng-0 nr=0 


=+00 nr=+œ0 


У; 2b. ,2ng (UM CE Jy, (А, 2п ,rYSin(21n,0)Z m. 2n, (z) 


ng-0 nr=0 
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Résumons donc les deux résultats, pour la partie radiale homogène ci-avant, et la partie axiale 
homogène en notant, on a T(r,0,z)-> Ta(r,9,2) : 
E cama CosQm8) _ ( r Y” " 
ng-0 a, Ya, 2n, t BL) 1, 
F пу=+о GA o Sin((2n, +1)0) r Si А 
; RZ 2 21 ng=0 (а, (2л +1)+ B,1,) r 
5 Bi - fi == de Х+ | 2n9+1 
“ К S с! ,Cos((2n, +1)0)( r ig 
ng-0 (a Dn, +1)+ £4 er E 
„у Eno OMEGA), б 0 o Sin(2n,0)( r 
пе=0 (a,2n,+B,1,) 2n, + 8,1) lL, 





= 
ч 











nz T 





Тү 


ng =+% nz=+0 2ng,nz ` 2ng 


Fo a, À, I dn l L| В, Lon, =|, 


E (9 in, +1)9) e=) 


2 


+ 
L QAI "n" l L, |+ Ва l =. Ü 


(A (eegend "E 














Hosch nz=+œ0 



































4 
ДЕ 
T 
ig —90 nz=+00 iem Se "со, + Dl са 
+) У S = A 
ng-0  nz-l L a I 2 E UEY J| f, L, : пал, 
ө E " L 
: 2 >" 27 ПІЛ | 
m L a An 1; 2n, L, 01, L, 
D (e 
Sinon — 


учы ү 


элу nz Tang (Æ.r)Cos(2n,0)Z,.(2) 

2: 25 A : 
r E= 0 dE) Ierger (PY) 
GE eom lama Ps Sinn, +10), (2) 

Z LOG Aus D + Ban AD) 
ei e Gaeren + 10)Z,. (z) 

ng=0  nz-0 Z (z) (0, О) 
Ng 90 пг=+о p sa Tar, (A; r)Sin(2n,0)Z, (z) 


+ 
F wm ZG lo ALES CD B.E GSM 


r^ "nz Lon, 


ng =+0 nz=+œ0 








2 
T,,(r,Q,2) = — 
T 
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| EE | 
La solution totale s'écrivant : (+ 0, z) = T, (r,0,z) - T, (r,0,z)- 


Pour fixer les idées donnons un exemple trés simple avec trois températures distinctes sur chaque 
surface du cylindre, pour un probléme de Dirichlet suivant: 

AT(r,0,z)=0 T(r,0,z)fini СІ. а =a;=a; =0 В. = Bj = fj; =1 

T(r,0, i = f,(0,z)=T, T(r,0,z), = fo (r>0) =T, T(r,0, = f (r,0) = 


La partie homogène radiale donne le résultat suivant, sachant que les deux conditions aux limites 
sur la section haute et basse du cylindre sont constantes donc de symétrie Y+,X+ : 


(EN l x) + BJ (4, d ) = 0 > Jon Li, ,Qug r L)- 0 











Détermination des valeurs propres о, X, 


nr,2ng Jy, nr,2ng"r 




















js r Jang ds gd dO fiy. x. (r,0)Cos(2n,0) ja r Jan, Aran] dO Cos(2n,0) 
Qu. _ 0 re 
ER (1 T Oo (CS ' Aw angl )) | ( + до Tos Qs, anb. 1) 
Ch la, T: i Se kel, ә рет 
ОТТЕ COR rade eco 
jo Yo ens d dO fj, x. (r,O)Cos(2n,0) jo r Jan, ras, nj d0 Cos(2n,0) 
Y+,X+ z ^ 
Mi ( + GT Jo J ang (E dis, À + FN (4, РТА y) ^R + Ont Jo J an, (A, anl )) T (5, 0 nr ans y] 
hi 





=ô T = д Т === 
O,nr,2n9 210,0: 0 2(J, (2 n LYS dE 2А. oJ (2, ol, ) 


J an, (Ar Aner 1.) = 0 
Autrement dit les seuls coefficients non nuls sont ceux de la valeur propre n4=0, et l'on retrouve la 
méme solution lorsque le probléme ne dépend pas de variable angulaire (dimension rz) : 


Ad tq J o^, "9 ) =з 0 
d E JA ar) Sinh(A, о) aps A or) Sin hlar, (1. — 3 


1)202 J, (À: 


nr,0 








T, (r,0,z) = — T, 0 
L, À. a gh (Ar. ol, )Sinh(A,,. ol ) nr=0 ^ar, oJ (A, ad, )5їпА( А. ol ) 
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Passons maintenant à la partie homogène radiale, dont également on ne retient que la symétrie 
Y+,X+ : 
_ MZ л 


SDT, Z„(D=SinX.2) X. tq т(,) 05А, = | 


nz z 








Ï 
2.0) = = 


2 


Deer Y+,X+ =) 
nz 


2ng,nz nz'z ng 0 2ng,2nz `` 


л 1. 1, 
ce =T, f dO Cos(2n,0)[ dz Z,. (z) = x Т,д„ „| dz 5їп(А 1) =T,ô 
0 0 0 


nz 2 





LS IUS, SIL) "P (2nz+1)x m. 
2 


л 2% Des 41), 05, uL) ECC "a 


2 


— T, (r,,z) =T, 


mr Or 








Ee A) : 
l, n„=0 Ао sol, )Sinh(Z,, ol.) nr=0 Ass oi (Un. l, ) SinR(A,, 01.) 
E. 4 ne I, (45, P) Sin AL) 





"e E Qnz* UI OS, ul) 


Ce résultat ne présente pas explicitement de dépendance angulaire du fait des conditions aux 
limites de symétrie Ү+,Х+. Dans le petit exemple suivant, on va tricher, et supposer de manière pas 
forcément réaliste que les fonctions limites du probléme de Dirichlet présentent des profils 
sinusoidales angulaires de symétries Y-,X-, comme ceci : 

AT(r,0,z)=0 T(r,0,z) fini СІ. а =а= о? =0 В. = В = В; =1 

T(r.0,z)| ,. = T Sin(2n; 0) T(r,0,z), , = Ty Sin(2n50) T(r.0,z) = T; Sin(2n; 0) 


tq JA, 4,)=0 
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Ces conditions limites injectées dans les formules établies plus hauts réduisent automatiquement 
la double sommation à une simple sommation pour ne retenir que le terme harmonique 
sinusoïdale de même fréquence, à savoir : 

Fr ' £ r 
Se tq J an, (4, ns lr )- 0 J, '(x) = СЕ SCH (x)— J vil (x) = J an, (4, no Lj- E J ang (А, 2,2.) 


Zi, nr,2ng (2) = Sinh(A,, an, 2) Zo, Ze (z) = sii À 2ng (L E z) 






































AL tq SinAL)=0S 4, =Ê et Z„(z)=Sin(A.z) |2,02) = = 
v == === јел. jeu 
X Sinh(A, nr ,2ng L) s 
Í der Jy, (A, sr) | dz Sin(2n;0)Sin(2n,0) | deed 227) 
Сол, = = T. ` r x T > On : r 
(J Jua (А, оп, LY) Г an D} 
1, 
| С А | dz fE(0)Sin(2n,0) Í drr J, (A, т) 
Y-,X- z T 0 iba: 
Сол nr,2ng m = T; 2 д лө r } 
(J Js „(4 anb. D Us E) 
1, L 
En posant Hj. 2nj - [ar 5 Jus (Ан г) NW 2nj - jar 5 J, „(А 2nj r) 
0 0 
Dran Dran 
Сы, = SCH p 20, : nr,2n; E — T; T d . nr,2ng 
2 FE (4,2 ei 2 ГЕ Ub an; L.) 
z^ 
Сэ, A z Lo. A dz Z, „(®) > p 2а = =0 с us Я 
ө» nj, lg 2 0> 0> Ла nj, o 
r r r 2 aNg T 
А anf tq Ju (A, KI Bye 0 ^, 2nj tq Ju (А, 2nj L)- 0 A. 
m- Sinh( A . .z) Tj an; (4, r) 
ТЕ "їл Ei птг,2п{ пг,2п{ Sin(2n;0)+ 
ims BILA 2nj L) = к nr,2nj l.) 
100) т Si ax re An 1 (А ) 
ё nr=+ Sin [j NA 
+ TZ nr,2nj nr,2n 2nj nr,2nj Sin(2n;0) 


пг=0 ЛАА, о L) le +1 (4, 


nr,2ng Y 
AT. | == L y Uo, ar)Sin(2n;0)Sin(45,, 17) 

p к nj. 
l, nz=0 Аа, ‚ (22,1, ) 
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Prenons le cas où toutes les harmoniques sinusoïdales sont présentent et où les conditions aux 
limites inhomogènes présentent une symétrie Y-,X- , par exemple, comme suit : 
AT(r,0,z)=0 T(r,0,z) fini CL. a, =а=оа =0 В, = В; = Ві =1 


Т, oeo Z] Т, ERI 
2 2 

T. Oe E a T. 0 E 2л! 
2 2 d 





7(7.0,2),, = 70) = 








Ту LES rs ZERI T7 eclaf =i 
T(r,0,z) ә = fo (0) = T(r,0,z)., = f;(0)= 
` 3л Е Зл z ; 3л А 
Т, er. = Г, ERI T, 0 c| r, Z - T, 


Prenons alors les termes des séries uniquement selon la symétrie Y-,X-, il vient : 
r Le V r r 
À tq J an, (A, on, 1,) = 0 J, (x) Е x J, (x) A Jua (x) = Jy, Ars 2n L) = -J an, +1 Anny L) 


nr,2ng 
r,2ng z) Zo 2n; (z) = Sinh(4^, a, (1. z z)) 


_ MZ T 


A. tq Sink) =0=>K, =Á et ZG) Sina) | 


Z, 2s, (2) = Sinh (A 


Ï 
z-o = S 








nz nz z 


2 








1 ` 1—(—1)* | ` 
EM Se dz Z =  Posons lj, —-l|drrJ, (А „г 
nr,2ng Sinh(A, an, L) | 2 nz (z) AS Ü nr,2ng | 2ng ( nr,2ng ) 

л л/2 Т? 
(dz fr (6)Sinn,0) = 27; | dz Sin(2n,0)= — (- cn") 
0 0 Ng 

л e | л ; | 2p: 
= [dz f (0)Sin(4n,9)=0 ег [а fr (6)sin((An, +2)0)= => 

0 0 ng + 


L, л л 
[arr A, Xuan) dz SEOS) D. | а уу (O)Sin(2nç0) 
22020, 0 0 

















o я 2 = ; 2 
PNE (OU Kë uon 
[^ л л 
[arr A, Xuan) dz EOSIN) D. | а f; (O)Sin(2nç0) 
C 2 2 7 = | ; 
ЖЕ ЖЕ И) 
л 1 л 
A 2 
ср = | dO f, (6)8inQn,8)| dz Z„(2) = es =0 с = r [40 j,(0)Sin(2n,0) 
0 0 2п2+1 0 
Ci =C mam Gru Si» 2p 
Y- X- E 2I Ts dnd Y _ 2T% hs 
Lar,4ng+2 `` 0,nr,4ng+2 `` 
2n, +1 ER CAE А) 2n, +1 Ge TENE 
еї A" 4Т; 
4ng9+2,2nz+1 P (2n, + 1) 
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Ce qui donne la solution : 


Р 2п2+1)л У À 
а= Crete et A, an, +2 tq Jan +2 (А, ал21.) = 0 


2 


1, 
r Li r 
TI] an,+2 n [ar r JA (A, ano) 
0 





"rere Sinh(X, p 
T. ш SE "o J ang (дп „r )Sin((4n, + 2)0)+ 
т 0 ng-0 nr-0 Sinh(A,, an, ol.) (2n, t Y Чуу аз (4, sual 1)) ° g 
Z í 2 р rl? METE Sinh(A^,. 4ng*2 (I. Ge z) Tin ,Ang*2 





+Ту 

Z nr=0 Sinh(A;, ke L) (2n, + (у JA (2. Ana tir 1.) 
16T. OAM 4по+2 (45117) 51п((4п, ч 2)0 )Sin (a2) 

T° ng-0  nz-0 (2n, + UE + 1) 4ng42 (Аза, ) 


Pour calculer les intégrales des fonctions de Bessel оп utilise des formules données en fin de section 
que l'on reporte ici brièvement. Il vient : 


(2) Í dz 2] (2) =[zJ„„(2)]+v Í dz J (2) > Í dz zJ,(Az)- = Í dé CIC) avec Q-AÀz 


1)j 5 J anys2 (i, an or) Sin((4n, + 2)0) 








(4 јело 23 nm алаа ano Data 
CCS T sU on nr жый, 


r = e 
Des Rad Í dr r Jas, +2 (А, ansia) 
0 


1 r r 
= [E Jo, 4по+2 L, Joss (А, „4ng +2 1.) + 4( (2n, + 1) D J ang+43+2k (Z u ,„4по+2 л) 


г 
= Уап, +2 ( 
mr, Ange +2 


Ре тёте 


1 , 
Ij. nf -far 4 J, „(4 2ng r) © à em LJ aU, nr,2ng 1.) ^ 4n, 30 +1+2k (X. ,2n d ) 


nr,2ng 


I r 
NW 2a -fa d J, „СА. 2п; г) = r [as 1.7, 2 „ CA nr,2nj 1.) j^ 4ni Se +1+2k (X. Ani d ) 


nr,2nj 
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Le cylindre creux complet : 





Problème homogène axial sur le cylindre complet creux 


Commençons par le problème sur le cylindre creux et complet où les conditions sont homogènes 
sur l'axe 7: 
AT(r,0,z)=0 T(r,0,z) fini 


CL. a„T;(r.0,2)-B„T(r.0,23), , =-f1(@,2) Œœ„T.(r.0,2)+B„T(r.0,23)),_, = 7,0,2) 
osT,(,0,2)- ByTQr,0,z2] ,=0  afT,(r,0,z)+ B;T(r.0,z)|_, =0 


Ce problème va également se décomposer en deux sous-problèmes, 
AT(r,0,z)=0 T(r,0,z) fini 


СІ. a„T,(r,0,z2)-B„T(r,0,3)) , =0 a„T,(r,0,z)+ B„T(r.0,2)| _ = f.,(0,2) 
aiT,(r,0,z)- PiT(r,0,2)) =0  ajT,(r,0,z)+ B;T(r,0,z)| _._=0 


et 
AT(r,0,z)=0 T(r,0,z) fini 


СІ. a„T,(r.0,z2)-B„T(r,0,3)| , =-f (0,2) — oT 0,2) B„T(r,0,z)) , =0 


aT. (r,0, z)- BsT (r,0,z). .. -0 aiT.(r,0,2)+ B;T(r.0,23)) _, e 
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Le systéme de fonctions propres axiales et angulaires sont identiques au cas du cylindre plein, la 
différence essentielle provient de la partie radiale, qui maintenant va comporter les fonctions de 


Bessel de premiére et deuxiéme espéce. 


Voici le calcul de la partie radiale pour les deux sous-problémes : 








Rua) m s, (Aa) + Bj „К, Aie) 

a„T,(r,0,z)-B„T(r,0,2)) _ =0= К? „ (r) 

2,4, ГАНА (Qd ) + Buy inz Kn, "Ah ))- Bn Li „L, (42.11) + B, nz K n (Anel m ) = 0 

> А Gal, (Ala) Е BAL, Al ))+ Бн Dok, (Qa) ES В.К, G-A, ) = 0 

I. (А „г K, (Ar 

See 2 ' al = | 2 ) | 2 ' = = 2 | 
Aril n (4,11) Ke Pal. (4,11) Ang A „К, (Aila) E В.К, (Al) 

«„Т, (r,0,z)+ B,,T(r,0,z) PRE 0> Re (r) 

I, Or) K, (Amr) 

> Ranz (r) 9 2 ' 2 , Z m 2 ' Z ч 2 

Иа (4,1,5) + BJ. GA ( nz X „o K n (4.1.2) + BK, (Ara ) 


Ces deux fonctions radiales introduisent un terme de division du fait de l'application successive des 


conditions aux limites inhomogénes еп г: 
AT(r,0,z)=0 T(r,0,z) fini 


CL. aT, (r.8,z) - B,T Qr O,z) 4 =-f„(0,2) 

















giis I„ Ar) K, Ar) 
r) = 
LÉ (0,1, (Qa) * Barn, (Al) Lok, (Us) * DK. (Al) 
Lo KA 
P. n (r) = À T 20 + =”) z 
А (Aou, (Als) Brat n, PELA) (ak, (o L4) + BK, E) 
rl 5 A Anz L m Al) n BAL, (a) „К, (Ala) = В.К, Al) 
у (0,1, (Ql) + В,1,, (Al) ak, (Qa) * В,К,, (Al) 
C.L. oT, (r,0,z)- В„Т(т,Ө,т)| = f:(6,2) 
R”? ( ) Lb (Ат) К (шыт) 
ng,nz r)= 2 ' de z 2 ' 2 2 
( nz X als, (Aila) BAL, (Al) Lok, (L4) - BAK,, (Al) 
І, К "(4 
ка = Aur) 
° Loi, (Qa) = BAL, (A, )) (а.к, (La) В.К, (A7, )) 


r2 _ a, An 1 (EWEN F Bal n, (A) a, Anz K n (EIER zs PK, (Al) 


nz^ пө 


"* Brad, (Eda) Bralna) ULeuE, (А 1„)— BK, (А 1„)) 
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Il ne faut pas oublier les termes de valeur propre nulle, dans le cas du cylindre creux. On rappelle 
qu'il faut ajouter le développement lié à la valeur propre nulle sur la composante axiale z. selon un 
raisonnement analogue à celui réalisé dans le cas cartésien, le terme axial de valeur propre nulle 
admet une solution non nulle uniquement lorsque les conditions aux limites homogénes sont de 
Neumann de part et d'autre des sections du tube. Et de plus la composante axiale admet 


2 


: А А =0 { A ; 
uniquement la solution constante. Comme ^” , la partie radiale n'est plus une fonction de 
Bessel uniformément nulle mais la solution du probléme de Laplace bidimensionnelle polaire, soit 
en ajoutant à la série précédente cela donne : 


T,(r.0) = S (Ayr + Bir"? Kat, Cos(n09)  B^,sin(n;0)) 


ng-0 


I n'y a plus de condition de finitude, les fonctions radiales et coefficients associés s'établissent par 
le respect des conditions aux limites : 


Cas (1) К, (r) = Ar Br? a aR, r) BAR: (r) 


n 





г=1,1 = 0 


=> nO QA (arro = Bb Á LiB, [p t sss 0= Aal" (n8 Q4 — Li, ) = Bua. (n8 0,1 + LiB, ) 


nO -nO 
R... (r) = Logo. +!1„ B а +(пӨа,, в.) 


Саѕ (2) А, (к) = A rr? + Bir”? => @„„КЁ„|'(г) + B. R. (r) 


rl 


= 0 


r=l,2 





=> n0 0.5 FR E BT )+ LB, ы + в") 0 — A; ol," (n0 a, T Lila) = Big (лө 0.5 EX LB, 


nO —n0 
r F 
Rar) = (n0 a, DESI *(n8a,, =L, | 


r2 1, 


Оп doit même envisager le cas de la valeur propre angulaire nulle, qui conduit à une solution 
radiale en A+B Log(r) : 


r 


Cosi R... (r) = А + ziel. 





| a, Rss (7) V В.к, (r) et SH 


rl 


Sas = D.A z 0 => LAB. М ES „В, 


r 
Ro (r) - Qa + Bod E] 
rl 
r is r ' 
Cas (2) Rao) = A, + B, roef +) > d Ras (r) t Ko GI. =0 


r2 


= QUE T B240 2x 0 = L2B,245 =a -a, B, 


l 
Rao (r)2a,, eL sto 3 


r 





Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


Par exemple les coefficients pour un terme Ү+,Х+ du développement : 
Symétrie Y+,X+=> Ө, (0)- Соѕ(2п,0) |ө, (6| =л(1+д„) 
Cas (1) СІ. а,,Т,(7,0,2)+ BT (6,2) _ = f:(0,2) 


2ng -2ng 
Ï I 
Ra (L5) v (2n, di + zs Е | + (2n, Qu TA us 
rl 











l, 1 


2ng -2ng 
2n l l 
EMU E (o. Quat zx 3 E (2n, ал dE) | 
12 rl rl 
p De 
a, 2n,| (2n, €a +L Ba j - (2n, a, -L, B, "m 4 
1 rl Fi 


r2 1 j l 2ng 1 -2ng 
+ Bola C d Se Lib: [2 + (2n, di ш La B. | S | | 
rl rl 


л 1, 
f do dz fazy, x+ (0, z)Cos(2n,0) 


r2 0 0 


nef xl (1-0 


2 














ng DH r2 


Cas (2) a, T,(r,0,2)- B,TQ,8,z) _, =-f„(0,z) 


2ng -2ng 
l l 
К, (Li) = (2n, Q + Lal n | + (2n, 9,2 FÀ MU | 








1, r2 


2ng -2ng 
2n l l 
Ra) = A C 4.5 E 1,В,2 | i- | Е (2л, а, — ызы] | 
r2 r2 72 


2ng -2ng 
a, 2n, (2л, d „2 +l. B, 2 — (2n, d „2 - LiB, 2 + 
1 L Ds 


2ng -2ng 
rl l 1 
+ puan d 2 + NS | + (2n, d „2 == Lp 2) | 
r2 r2 


л L 
„ 140 | dz fay, x. (0,2)Cos(21,0) 
rl 0 0 


An „= С^ e anh, 
RA, T L 1 + до rl 
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Les coefficients du développement pour le terme А+В Log(r) sont : 
Cas (1) C.L. «Т, (r,0,z) + B,:T(r,0,z) 4 =f„(0,2) 


LiB 


I 
R,20 (1) = a, + леа) Ro (L5) x ` 
r2 


rl 


2л T 
9 | [аө [az f,,(0,7) 
D,50 = 05, LE + pe, + „бло z J ? Ay. = — 


r2 








rl 


Cas (2) œ„T,(r.0,2)= B„T(r.0,2)) _.`==>f„(0,2) 


l 
Raio (La) = 9,5 + Län 3 


rl 





l, 
| Rao 04)7- L 


rl 


Dig = 





р Кы feo 
|-<- 


Lbr 
1 л L D, 


l 
+ а. t ab aLog D 


rl rl 


Ce qui permet d'élaborer la solution du problème homogène axial pour tout autre condition aux 
limites homogènes que Neumann de part et d'autres : 


f4(0.z) т ЛДлу+х+\0,2) + Jiya, x (0, z) + Jun, (0,2) t Jux LÉI, z), sachant que 72100,2) = J (Ол +0,2) 
J (0,2)+ f, (œ -0,2)+ f, (-0,2)+ fa (0-7,2) 




















Sy. (0,2) = : 
fures. ш с>” (972 
far. (6,2)= F(0,z)— f(x -0,z) Ju? z)- f (0-7,2) 
Fri o dae qai 6) Lal 0,2)+ f (0-7,2) 


fa (0,2) = Jans (0,2) Tyre (0,2) + fox (0,2) + fizy- x- (0,2), sachant que fa (0,2) = fa (27 +0,2) 
/(Ө0,2)+ f, (m -0,2)+ (ОЕ (0-7,2) 




















/лзу+х+\Ө,2) = ; 
dade EU cone ACHSE 
Far x.(@,z)= f(0,z)— f. (7-0,2) Ha (—0,z)- f,(0-7x,2) 
О ор ылы a 0,2)+ f (0-7,2) 


Z (z) = ad А Cos(A,,z) + Bē Sin(%.z) 


0 ^*nz 


х. tq (оза - Bi Bi sima. = lapi +a; B? Cos) 











2 
п. 


рр (оссо) + Bi тал) t +2 lier À | |+: 
2 ol » аг 


Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


x [ + IN y 3 t ao Bi d 














ou MA (Bi Соз.) 0 





We I„ (Ar) K, On 
Ru, ý (Ea Ал „1 n, (Anz L5) + or. ‚(А )) ak, (EN L5) + B.K, пө Fam) 
RANA Ln Arr) K, (Ar) 





Wad ULL Bab 1.) (оК, OLLO- BK, (А 1„)) 
ar c Au GO LOB ALL). 5 OLEO BAKE, AL) 
"n Gran LL) Bal, C) an, KR 0L) 
BE ULL BAL EE 
n Ta at Al BA, HAN (оК, TO В.К, GE 








lz lz 
= | dO Cos(2n,0) | d ZO Gb = | аө Sin((2n, +1)0)[ dz Z„ (2) fy (0,2) 
0 


T 


T 1, 
nan = | dO Cos((2n, +1)9 j ez aheno d = | 40 Sin(2n0)] dz Z, (2) fuy (0,2) 
0 


0 






































слее Z fao coseno] dz Z Хе) 0:2) ie = fao Sin(( (Ол, +1)9 je 2 (2) far e (0.2) 
0 0 
gu FA Cos((2n, +1)0 j dz Z„ (z) all Gun = fao Sin(2n,0) | dz Z. (2) fiy x.(0,27) 
^ 0 0 
ree ez Cie 0) 
no-0  nz-0 (1 + д, „М, 2ng nz 
| ee Sin((2n, -1)8)Z, (2) R5, aue " 
acht nz=0 Z, (z yr DU nz 
+ b» Tex Cm x Cos((2», + 1)9)Z„, (DR a: nz 
= dm ZN Das 
кзы; е E Sin(2n,0)Z,. (2) Ge 
"AA Kat. 
T,(r,0,z) = — 


Л | по=+0 nz=+œ 248 Se ^*Cos(2n,0)Z,. DR, „(т r) 


ng-0 п2=0 (1 + д JS. D nz 
ng =+% nz=+00 с” FA Sin((2n, + 18)Z,. (2) 85, a, nz 
+ 






































М 2лу+1,лг 
пе=0 nz=0 zo DE a nz 

ES E Солу) (RS au 
mm ZE) D as 

NI пг=+® С. 12  Sin(2n,0)Z„ CRE Se 
пе=0 nz=0 zi GT ID nz 
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Auquel il faut rajouter les termes de valeur propre nulle, dans le cas des conditions aux limites 


—n0 
EAM 
1, | 


homogènes de Neumann de part et d'autre. 


nO 
R. (r,n0) ES (n6 o. tL. D L (n0a,, ES 


r2 


nO —n0 
r 
R, (r, n9) = (n8 4.1 $i LiB [z + (n8 4.1 ч LiB [t 
rl rl 
M e m d 
„п (п, а, +13, 1. -(n, а, lb, T. + 
r2 r2 
D, (n0) => | ў 
In La i La И 
+ Bol, (n; а, * Lip, i + (л, 0, — LP) T 
r2 r2 
5 * 1,5 d 
Q jns (л, а, LiB, P -(n, a; - LB) [о + 
1 rl rl 
D, (n0) = 1. j ER 
ra 1,5 i 1,5 | 
+ Bl; (n, d +1168, d * (n, Q - lf, Ca 
rl rl 
л L T L 
ANIS = [dof dz far, x, (0,2)Cos(2n,0) А10 = [ав[а E (0,2) Cos(2n,0) 
0 0 0 0 
л 1, л L 
Аю = [a0|dz fay, x-(0,2)Sin((2n,+1)0) Aus = f do| dz fr. x (0,2)Sin((2n, +1)0) 
0 0 0 0 
л L л L 
Anas = [d0| dz fay- x. (6,2) Cos(Qn, +1)0) As = [a0|dz f.a y. x,(0,2)Cos((2n, +1)0) 
0 0 0 0 
л 1; л 1, 
Ang" = fdo] dz fay- x (0,2)8n(2n,0) А4521 = f dO | de fy- x (0,2)Sin(2n,0) 
0 0 0 0 


l r 
Rao (r) = 9,5 + ab, Log 2) Rao (r) — di + Ld Z) 


rl 


Dj, = 91 йа а, + ЖС) Ро 70, =”. S sl, + Lb tod L, | 
1 


l, rl r2 rl 


an = аеја f.(6,z) А? - [ae] a ICH 
0 0 0 0 
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КЕЕ Е nz T 2 L 
Si B; = В, =0= Z. (z) = co[ = Iz. (=) = 5 
T. ТО MERC. = > o. +L,B, Lo Ө ПЕС |, +1 B.Logl = | |+ 
zh = Zä 2л LD r2 r2F r2 g r 2л LD, rl rl rl g L, 
n= ДЛ Cos(2n,0 naro ДХ Sin((2n, +1 
Y 2ng,0 ( 0 Is (r,2n,) E: Y 2ng41,0 ( 0 y) Rm (ғ,2п, + 1) 3 
45 (+д„„)р„0п,) Zi Dn) 
ту=+ю ДА+ Cos((2n, +1 noto А2 Sin 2n,0 
2по+1,0 (( 0 y) R. (r.2n, +1)+ >. 2ng 41,0 ( 0 le r2) 
Ex n9 -0 D, Qn, +1) ng -0 D, (2n,) 
T l. | чә Aja o” Cos(2n,0 "C Aus Sin((2n, +1 
ES ` 2ng,0 ( 0 Is (r,2n,) Y 2n,+1,0 (( 0 y) В, (r,2n, +1)+ 
Z «0, „)Р„(2п,) = D ,(2n, +1) 
nct» ДУУ Cos((2n, +1 чеке ДУ Sin(2n,0 
2ng41,0 (( Ө y) R, (r.2n, X dr Y 2ng 41,0 ( 9 ) R (r2n) 
SC D,,(2n, +1) FE D, (2ng) 
eege Е ` Je. (r) 
< + 
ng-0 nz=l (i6, , um 
E rl,Y+,X- с: rl 
ME ER Con +1 nz Sin((2n, +1)9)Z„ (2) Ron, de 
+ > 2. D” 
ng-0  nz-l 2ng-l,nz 
ng =+% nz —oo cn Cos((2n, + бсо = а je. +1,п2 
+ < + 
2 п2=1 De 
с BEEN = 2 Je. (r) 
+ < + 
2. nz-l D d 
" 4 
TL TT Cos2n)Cos ka 
+ < + 
ng-0 п2=1 (1 сЕ Ô о Jp. 
ics Ren, Je Jc 
+ < + 
2, nz-l D pe 
Tig =+% nz=+0 ae Elte, gel jrs 
T Y r2 : T 
ng-0  nz-l 2n,+1,nz 
| «ez Sn соң er " 
Hg =+0 nz —--oo 0» 1, 0» 
+ 2: D? 


ng-0  nz-l 





2ng nz 
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Problème homogène radial sur le cylindre complet creux 


Passons maintenant au problème homogène dans la dimension radiale : 
AT(r,0,z)=0 T(r,0,z) fini 


CL. a„T,(r,0,z)-B„T(r.0,3)) _ =0  a,T(r,0,2)+8,,T(r,0,7)., 
aiT,(r,0,z)- BiT(r.0,3)| = озб) aiT;(r,0,z)+ B;T(r,0,z)| _ =f (r.0) 


Cette fois-ci, ce sont les fonctions propres radiales qui vont être substantiellement modifiées. 
L'indice angulaire est entier de part la condition de périodicité 2n imposée à la solution. Et l'indice 
radial est déterminé par une équation transcendantale obtenu en appliquant la condition 
homogène sur la surface du cylindre, comme ceci: 


oT, (r,0,2)- B„T(r.0,2)| _ = 
no Anr) Y, (Ar) 

„АЈ. ng (Ad ad Bad, (Ar [у a, M, Ï, ng '(2,, LJ: Bar, (4 nr ml) 
mo. + B,,T(r,0,2) E 
— a, D ns. sng“ пө "A, ‚лө 15) + Bd, (А, ‚лө La) a, Aarm Y, 519 ng "A, Jg Lj) + DH, (4. Jg 12) 

a, PI ng“ ng "A, ‚лө E E B.J. (2, nr mln) a, Vor n Ï, Mio" ng "A ata) Bi BY, BEN nr ng 11) 
a tq a, 2r ns. ng“ ng ' (A. ng r 12) + Bd, (А, ‚лө la) _ a, in io" ng (4. К r Lj) + BY, GR e L3) 
"s a, А ‚лө Ј (E ‚лө у Bad, (А, Jg 11) a, Va wl io" ng (E e ln) - Ba Ya, А лө 11) 


REH A Datt Ү, (А) 
‚п E = 
im. Qa Ao ns. Jug" ng (EN ng b Bua n, (А, ‚лө La) KEN nr,ng^ ng Y. EI ‚лө L)- BY, А ‚лө 11) 





=> К, (V r) = 














Termes de valeurs propres nulle 


Comme nous le répétons, les termes de valeur propre nulle appréhende les solutions possibles du 
problème aux limites lorsque les valeurs propres définies plus haut sont nulles, soit quand la 
solution ne dépend plus d'une ou de plusieurs des dimensions aux conditions homogènes. On peut 
dire que le problème dégénère en problème de dimension inférieur. Ces termes représentent la 
part de la solution de dépendance dimensionnelle réduite. Lorsque c'est uniquement la valeur 


r 


А . Au „ÍO =0 : A 
propre radiale qui est nulle, " a , la solution radiale est de la forme, 
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жа (n) R(r)=0 R(r)- A, r" +В r^ 
r? 


dr? r dr 



































У er rl rl 
Condition = 
r2 


angl, (4, p = В D js Ba (4, Lr + B, Ж )= 0 

a, ny, (A, E e В, b n B. (4, E Б В, n )= 0 

A, en" Lh = Babi" )-B, „la Baha Sba Lu 

A, (a, Jub LU" )-B, о> spaa” )= 0 

(„лы um Bu = o, nl La Wi Bala _ 

a, n l, "Ba" fols аот ' 

= adr I '- Bala” в," – E на") (o, mut a + Bala ” Lk angl i THE 
= а, В," Lo De Baur ала, т "hug ; Noc Чы" = 


no-l] nọ- —ng— пө ng ng-l -no ng 
+а,10:,5п А 12 ` Fd „aB. nl rl 15 + Bad, tgl. "1° +818, ‘ls 


по 1 ng-l пө пө 
2| , ai) EM l l 
= QA „UX 579 | la | 2 | -Lla |= + 21 j] |+ 
rlr2""0 | rl É | r2 H | sell m | H | 
ta, nil, DIE + а + Bond: hi + Le 
L, La 15 L 


Sin, оа (5 (I ]] аео, =0 ou B,,=0 






































rz 


=> 


А = 











la) a 
Une solution radiale non triviale n'est donc possible que lorsque simultanément il y a valeur propre 
nulle radiale et angulaire. Les termes de valeur propre nulle radiale et angulaire 
A, 79 et n,=0 


nr,ng 


. La solution radiale est de la forme R(r)=A+B Log(r), injectée dans les 
conditions aux limites, la solution radiale est constante, et les calculs simples conduisent à la 
nécessité de conditions homogénes radiales de Neumann de part et d'autre : 


d'St), 1 ак) 





=0 R(r)-A,- B,Log(r) 





dr? r dr 
Deux cas В. =0 ou B,,=0 
a  R'(r) к=! 
(D cas В =0 > À =0 


a, R'(r) + BRO), 





>a, =0>38,=0 et В„4,=0= B, =0= R(r) = A, 


rl 



































(2) cas B,,=0— e rt 


r2 











dë e 205820 er PASS VSA, 
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Cela conduit donc à des termes de série indépendant de r. On en conclut que le terme globale de 
valeur propre nulle est linéaire: 
T,(r,0,z)= Aj + Bj z 


aiT,(r,0,z)- BT (8,2), _ =0 = ai Bi = B; 4 
T,(r,0,z) = Ajo + Biz) 
Cela donne deux termes de valeur propres nulle : 
Si Ba = В, E 0 
Jarr fao rre ког 


2 stoe + Brai + BERE) “9 121, stoe DECH 
Terme i - + ps z)- Ai (а; + Pí (I. - z) 








2 
А; = 


Revenons aux termes де valeurs propres non nulles. Les normes se calculent ainsi : 























d Ju eas І Li (cial) 
аА, оТ, Lë Sala) Bad ang Aar Sala) Ar oaa Car zua) BA (оь) 
SG CN СА 
аА, an d ang Amr angl) Ва, (ьа) a 2ng Lang (А2) Baba (пл) 
te Jul Sea) | La Air Sala) 
QA, ољ Јн (4. 2л Г) Ba ança A, nen) la) vu ЕЗИ (4, Anen) lg) = Baba A nr 2ng la) 
AN == Jina Arnab) Tal Sanel 
QA, ui aat (4. Zuel ln)- Badan (A nr 2ng H 11) Ea Au Tang (4, unt 1)- Bab, aU; rng 11) 
Riran (a) DE Gi Jan, ay an li) Yan, s osi) 
d. л scho UR EE 
Ry, 052) 2р" (L.)= 1, 22) Yon, (пі) 
m TONERS OX PERPE Ee 
Roma (а) pt ens Jua (һј) Ya Wa zwa) 
r 


r nr,2ng4l r r r r 
Аъ CAT „219 Jua (4. mala) Badan i (A4 nr,2ng4l Li QA, ,2ng +l Dad (A, aula) - Bas avos ala) 


К, ,2ng4l iU n) = DE (1 )= Jua (4, ,2n9+1 15) Tun (4. ,2n9+1 la) 
r nr,2ng4l V r2 r ei r 
A Set QA, ,2ng Ja "A, mal) Bad angi Ar, mala) QA, ,2ng4l Ya (4. ,2n9+1 LI: BE (A, ,2ng4l li) 





2 


2 
LA o of Ra р | oS af (a a 1 lr e 
r2 | пг e l, ( = - 











LR = ; 
2 2 
le, ES 1 Daf Ile, E lr Dal, 
| of ES 1, ck nr nw l, A nr ins 
nr Aner) I 
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Voici les fonctions axiales, ainsi que la décomposition symétrique des fonctions limites : 
Zar ng (z) = Ag an, Cosh(A, 2n z) + BoSinh(Z,,.3,,2) 

тарау (2) = AG Apr зп, a COSA, 21, az)* Po Sinh(A,, зл, 12) 

Zo,nr, 2n (z) m агн Cosh(4^, a, (I. Б г))+ Bi Sinh, au, (1. e z) 

Zona (2) = AF An ar COS all, = 2))+ B Sinh, —2)) 

Jo (r,0)= Йод. (r,0)+ Ja, (r,0)+ É xs (r,0)+ Jazz (7,0), sachantque Je (7,0) = fy (r,27 +0) 


fo (7,0) fy т,л-0@)+ fy (r,-0) * fj (r,O =w) 














fissa (r8 = - 
и Jo (0) + fo (r. — 0) =k (r,.-0)- fó 00 -л) 
fr Je (0) - fo (r.m — 0) RE Jo (r.9 -7) 
fé, 00) = Jo (7,9) - (57-0) n (r,-0)* fy (r,0 7) 


Л (r,0) = firex (r,0) + Ire (r,0) + Руха (7,0) + fix (r,0), sachant que Л (r,0) B Л (7,27 +0) 
fi (0,0) - f (r,na -0)+ f;(r,-0)+ f;(r,0-7) 











fires (r,0) = 4 

en. frëen бз) (n Ө) E Dë të 
LY+,X-\ ? à 

ЖИ ут AL) NC A) A à 
LY-,X-N ? ? 


Jí (1.0) fí (r.a -0)= fí (r. -0) + fí 00-7) 
4 





fir- (r, 0) = 
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Les divers coefficients des séries sont : 


r 


nr,2ng = À » z Dz z oz А r 2..2 r 2 z AZ 
Gg Cosh(A,, з, L Yo; В, > а; Вг )+ Sinh(A,, 2n, l, Jaia; И ) + Bs В, ) 


1 
e". ng p 
"TT Jr COS am al. Xe Br +a BE )+ Sinha, Nasa (Xom) + Bébé) 


Ainsi que les termes intégraux liés aux conditions aux limites : 
























































Lo л 1.5 л 
[drr Ruan, (D dz fiy. x. (r.9)Cos(2n,0) [drr Ruan, (OÍ dz 7, s. (1.0) Sin(2n,0) 
CY+X+ — ln 0 Cr-Xx- da 0 
lnr,2ng 2 L,nr ,2ng 2 
(1 + a] [m (r) | К,ол, o| 
Lo л 
[arr Ry, a0 [ dz für. x-(r.0)Sin((2n, +1) 
(THX= = 11 0 
Lnr,2ngl 2 
Ка) 
DR А 
[arr К,а) [d fir- x. (r.0)Cos((2n, +10) 
CX 2h 0 
L,nr,2ng+1 2 
[ЖАК ()| 
12 л li л 
[drr Ra, Of dz fé,,x,0,0)Cos(2n,0) [drr Ryan, O| dz fiy, x. (7, 0)SinQn,0) 
Се. Е la 0 Y-X- _ la 0 
O,nr,2n9 2 0,nr,2ng 2 
Is, ui [Ruan CY] 
i л 
[arr Rascal dz fy, x, 0)Sin((2n, + 16) 
Crex- z la 0 
O,nr,2n9+1 2 
Raana O| 
ПУ Е 
[drr Ri, aO [dz fes x. (r.9)Cos((2n, +1)0) 
crx [ 0 





O,nr,2ng4l — | 


2 
Ron, +1 GH 
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Ce qui donne la solution du problème homogène radial : 


r 





wn, ЇЧ 
a, А p (E ‚2ле Lj)* Bd (4, EI la) _ a, А 2m Ton An 219 L5) + 8,7, 2n (4,,. ,2ng 12) 
Œ aAA 2n 2n, "A, EI la) — Brit an, (А, POA la) 9414, sal (EM 2ng L1) - Bar 2ng (Q5, ST l4) 
Se tq 


' L 
a, 27 "nr mama Ja (А2, „2ng alia) + Bd (Аһ опу) _ 


A 
а, nr,2n9 +1 Yn (EN nr,2ng4l L4) B, Put Ao 2n, ll 














Q AI 2n9 +1 Jangar (Аһ, sb = Baad ang (Aron ala) A AA, Е Р '( ae ра CRE 
Ron (r) = z х Jam (Au EN 2. = - | Yang (4, EI r) => 
14, EI J an, (4, EI? ks Burt an, (Air, 2ng Lj) a, 2ng Yan, (4, EI? Li)- Bar 2ng (4 nr 2mobrt) 
Rona) = S Jangi ass al) Ү, aas, al) - 
Q Any ән, Ju (A; nr опна) BJ a0 nr Zen abr n) 914, m Yu (A, 2091111) — Bab Li Sen il) 


T (0552 
T 





ng =+% nr =+œ0 


Y+,X 
>. У А» 2n C, ros, R nr,2ng (r)Cos(2n,0)Z, nr,2ng (z) + 
ng-0 nr=0 
ng —t9 nr=+00 6? 
ay ë 
+ 2, > Ay ,2ng4l 7 Lnr,2ng Jod (r)Sin((2n, Ed NNO NZ ina (z) + 
ọ=0 nr=0 
ng =+9%0 nr=+œ0 PE 
-,X+ 
+ 2 > Ay ,2n9+1 Lar. On A ewe +1 (r)Cos((2n, > DEZ à (2) + 
ọ=0 nr=0 
ng =+% nr=+0 ci- x- 
+ Z У А, ,2ng C; r,2ng E ns (r)Sin(2n,8)Z, „an, (z) + 
ө=0 nr=0 
ng =+9%0 nr=+0 Ys 
TAX 
T 2 > Am ,2ng Сол nr,2ng T (r)Cos(2n,0)Z 2, (2) + 
ө=0 nr=0 
Jig =+% nr=+0 C! X 
ENE 2 
+ 2, Y А, 2ng4l C, ,nr ,2Ng К,ол (r)Sin((2n, + DA end (z) + 
ọ=0 nr=0 
пе=+%® пг=+00 
C ,X+ 
+ SS Kä A;, ,2ng9+1`0 Conr e ot o dins (r)Cos((2n, + Zu (2) La 
ng-0 nr=0 
gs nr =+0 
CT ,X- H 
У E 279 C, ,nr ,2ng Ron, (r)Sin(2n,0)Z,,,, an, (z) 


ng-0 nr=0 
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Lorsque les conditions homogénes radiales sont de Neumann de part et d'autre , alors il faut 
ajouter à cette solution les termes de valeur propre nulle : 


Si Ba = В, =0 








12 2л 12 2л 
[arr [a6 j;(r,0) [arr [a0 f;(r.0) 
AZ = Li 0 APL Li 0 
0, =, `Z. 2 z Z A 2. 0,0 Z Z. 2 32. Z 2 Ж 
DIE Nai B; + Bras + Bi BIL.) LJ -1 Jai Bi + Bras + B; BFL,) 
r Jy, (EM ,2ng la) Yan, (E ,2ng 12) r Jua Gr, 2n9+1 15) _ Longs Gr, 2ng+1 15) 








tq 
nr,2ng q ' ' nr,2ng4l ' ' 
Jy, (A, ,2ng m Yan, (4, ,2ng La) Jua (А, ‚2пе+1 11) Yan, +1 (A, amail) 


J an, (Aar 2ng r) Yang (А 2n" r) R у= J an, a эл1) Yn (4, anal) 
nr,2ng4l m 
Jy, E ‚2пө P 4 (EW ,2ng е SE? Ја Arr, 2n9+1 la) Te 2ng4l La) 
D HN (Lj) = D, E (1,5) = D, n (1) = D. NM (L5) = 0 








Rp ax C ) = 








2 
lo 1, (R Ж ОО, D =La L (R ax Gf | 1- 7 
uf: E e 
2 
PI MER) е SM Жу (N _ 
r2 | nr,2ng4l Sc nr,2ng4l 





2 
Re +1 (GH z 


T4 (r,9, z) = A; (a; + Bj z )- Ai lo + B; (L. -z)« 
>. NA, 2ng o doc К, ,2ñ9 (r)Cos(2n,0)Z, эл, (z) + 
ng-0 nr=0 





ng =+%0 nr=+œ0 


+ x > Ay P GM Rp „2ng +1 (r)Sin((2n, 18)Z, oa (z) * 


ng-0 nr=0 


ng =+0 nr=+0 


+ Z 3 A nr stati d nr,2ng4l (r)Cos((2n, +1)9)7 un (z)+ 


-0 nr-0 


по =+9%0 nr=+0 


T Z Ул, nr,2ng Су, Rw ,2ng (r)Sin(2n,0)Z, „. an, (z) + 


=0 nr-0 


+ 
à |t 


ng =+œ0 nr=+0 


T 2 » A nr,2ng Ci R,, ,2ng (r)Cos(2nç,0)Z Zu (2) + 


-0 nr-0 


ng =+0 nr=+0 


+ > >. M CR B nr,2ng4l (r)Sin((2n, i DEZ. ind (z) 25 


ng-0 nr=0 


пе=+%® nr—too 


RE Z 3 A; nr TUNE ORAE: nr,2ng4l (r)Cos((2n, + П рай (2) T 


-0 nr-0 


=+00 nr=+00 


+ 2 > A nr,2ng Gun R. ,2ng (r)Sin(2150)Z, „m. an, (2) 


-0 nr-0 





La solution du probléme est toujours donnée comme superposition du probléme homogène radial 


et du probléme homogène axial : T(r,0,z)=T„(r,0,z)+T„(r,9,z) 
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Secteur cylindrique plein 








La section est définie par le domaine produit cartésien des trois coordonnées du systéme 
cylindrique oü typiquement on doit résoudre le probléme suivant : 
ОТ (ғ,0, 2) 10T(r,0,z) 1 OT(r,0,z) OT(r,0,z) 
2 + + 2 2 * 2 
Or r Or r 00 Oz 
Q=fr/0<r<l}x{0/0<0<0,<27}x{z/0<z<1} 


CL. a,T.(r,0,2)+ BT(r,0,2) = 7,0,2) 


=0 T(r,0,z) fini 





asT,(r,0,z)- ByTQ8.z), =- Cz) — a T,(r,0,2)+ BT(r,0,7),, = / (uz) 


aiT,(r,0,z)- BiT(r.0,23)) _,=-fF(r.9)  aÿT.(r,0,2)+ BiTQr 0. z) | = 370.0) 


Probléme qui se développe selon la méthode classique de superposition en 5 sous-problémes. Le 
seul avantage de la recherche des solutions c'est que l'on a plus à se préoccuper de la symétrie 
angulaire des fonctions limites, car la partie angulaire de la solution n'a plus la contrainte de 
périodicité 2л. On va séparer l'étude selon l'axe d'inhomogénéité des conditions aux limites : axial 
z, radial z, angulaire Ө. Nous noterons les solutions des trois sous-problémes comme ceci : 

T (r,O,z) inhomogene radial 

Tj, (r,.0,z) inhomogéne angulaire 

Т (7,0,2) inhomogene axial 

T(r,0,z)  T,4(r,0,z2) - T3, (7,0, 2) - T., (7,0,2) 
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Solution du problème inhomogène axial pour le secteur de cylindre plein: 


Le probléme inhomogène axial est le suivant : 
={r/0<r<1,}x{0/0<0<0,<2r}x{z/0<z<1} 


CL. a,T,(r,0,2)+B,T(r,0,2)| = 
cT, (r,O,z) — féT(r,O,.z), -0 aÿT.(r,0,z)+ BIT(r,0,z),., = 
aiTi(r.0,2)- BéT(r,6,2)  =-fé(r.0)  aT.(r.0,2)+ BÉT(r,6,2) = fé (r,0) 


Le ES EEN en série est de la forme : 
=+0 nr=+0 


T(r,0,z) = 2 Y (45, Cos(4! 0) + B?,Sin(72 


no 


6) (A; Сон, „ z) + Bi Sinh(X, D Va On 


=0 nr-0 


Les fonctions propres radiales se présentent sous la forme suivante, pour assurer la première 
condition aux limites radiale homogène de Robin, comme suit : 


Rar ‚пө (r) = Ju (2, пө r) 


o RO) ERO 0 ALS м aL de ut) B. Kg 


2 
Цина E] 











Run, CJ = ; 
A 
De plus J „ (s, „1)= das urge Uem e aor) 
ө К 
C H = Le L, y р, d " Jg L, } -2 i Jo Eë 1, у, -1 Ge. l, ) 
= Is, of = —— 





2 
Les fonctions propres angulaires sont les suivantes : 


®,,(0)= a228 Cos(2? 0)+ вет? ө) 


пө ng 


A. valeurs propres solution de Sin, ES 0 IER An, J ala? - pe pe )= A, Cos(4*. 0, ler as + Blat) 


kum Cos(22.0,)+ B Sin(22,6,)f + (929 sin 0,)- в°Соз(2° ef | 
Je, c» = 


2 exe od 4* Cosa 0,)+ вод ө, {оо Sinf 6, )- Bt Cos(42 0, )+ a Be 
Les fonctions CH deviennent : 
a;T. (r,0,z)— B;T(r,0,z) [ups 0 


=0 


m" | | Z nrn (Z) = XG Anr n COR, 
aiT (r,0,z)+ B; T(r,0,z) Ks = f, e 


nr ,ng nr ng 


z) + Po Sinh(A;, „ z) 


aT. (r,0,z) — BT (r,0,z) Sch => f, (r,0) 
aiT.(r,0,2)+ pjT(r,O,z) | «0 


z-l, 


=> Zon n (2) = AFR, „Соз, „(L =z))+ Br Sinha, (L, — z)) 
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Et les termes intégraux du développement en série prennent la forme : 


r 


1 
B eei Xo; а?) тњо, „l, Yei ai (X... Y + 6; B; )) 





j drr Kb, (r) Í dO f;(r,8)O (Ө) Í drr R„ (r) Í dO f; (r,0)0, (0) 


Warm — 2 2 2 2 
Is, C 6,4 (0) |... | 8,4 (0) 
Terme A, (Coun ом, 2) + Cora Zinn Ras, (P)O „9 (0) 


La solution de valeur propre nulle angulaire est constante et coïncide avec des conditions de 
Neumann homogènes de part et d'autre des parois latérales du secteur cylindrique : 


X... =0 > = 2 | | | 
nr,ng By = В . Les fonctions radiales sont les fonctions de Bessel d'ordre nul, qui 
sont déjà intégrées dans le développement en série de la solution (cas ne=0). Il reste à envisager la 


; 33 А„=0 et A., =0 et s e os ce 
simultanéité des valeurs propres nulles : ` ”? ш B,=B = В, . Soit le cas 


de conditions homogènes radiales et angulaires de Neumann de part et d'autres. C'est un résultat 
similaire au cas cartésien puisque sont exigées des conditions de Neumann dans les dimensions 
homogènes du probléme. Il vient alors une solution purement linéaire en z, indépendante de z et 
0, soit : 


Cas: Jed eB =0 








L,nr „ng 














i ER 

aT. (.0,2)- B;T(r.0,2) _ =—/у (1.0) 2| der |40 ra TE 
| G > À = — Z(z)= А2 (а? + В2(1 2 
aiT.(r,0,2)+ В'Т,Ө,2) =0 US Ol lar Bi + Bros + ВВА.) NO 
0, 

ajT,(r,0,z)- BiT(r,0,2) _ = | afars I fr (r,0) PS 

Pu 7 z > А = 2(2) = Aj la; + Biz) 
aiT:(r,0,2)+ BiT(r,6,2) = 720,0) Q4 Tai pi + Brai + BB.) 
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La solution du probléme inhomogéne axial est donc, hors du cas de Neumann bilatéral homogène : 
Va,,as,a?,B,,B$,B" incluant lecas В = ff =0 et œ =a; =1 
sauf DER Sel et a =af =a? =1 

dE r P Fr ' P 
Rd (r) =з Ju Eë r) los Iq UA o Ju Ж L )+ B.J 9 


2, 


А ГИГА „(е А) iz (4) 
Dé nr,ng r A, V Water К 


[Ran = 


DNA 








2 


À 
e Si s L, y b. б. l, ) B LA. 





2 
ou Is. (r)] = 





6,,(8) = a2 Cos(4* 0)+ P Sin(a2 Ө) 
A, valeurs propres solution de Sin(4? 6, DÉI одо? - pz: |+ a. Cos(4? 6, (Bla? + Blat) 


|'® (аел Соз(4,,0,)+ BE Sin(22.0,)f + (6029 Sin(22 Ө„)— Bà Cos(^? ө, y) 

„ге = 

18,50) =; Fa => oU, Cos(a2.0,)+ BE Sin(22 0, Næh, „тіле 0,)- Bt Cos(^? 0, ))+ o B? 
Ж (2) = AFA, oeh, „ (L, — :))+ Br Sinha, „ (1, -z)) 


Z rin, (2) = AG Army Cosh(4,,.,,2) + Po Sinh(A,, „2) 





1 
{ш Совт, „L Nas Br + a Bs) Sinh, „1, ро; (A, Y + В; В) 








jo ER S Zr f; (r.0)6,,(0) far Жын Zr f (r.6)0,,(0) 
Co, =2 2 2 = Сор = : 2 2 
Is, C) (Ө, (Ө) Is, ПӘ, (Ө) 


T(r.0,2)= У УА, (Conn Zomn E+ С (EUR, (reg 
ng-0 nr=0 
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La solution du probléme inhomogéne axial est, dans le cas de Neumann bilatéral homogène : 


Si Rebel uu. a, =aj =aj =1 
| | D 
2 r ni `, 
L, м, (A, n, L) i p Ee у 


2 





Кл (r) = Jo Cus r) Aa tq J (Ал, l.) = 0 Roaring 2 = 


























6,,(6)- Cos(2 0) et A tq Sin(22.0,)=0= 45 = 2а ег |Ө„(Ө) = gd 
0 
Zorn E) = AF, n, COSA, „ (1, —z)) B; Sinha, (L, - z)) 
Zi nr (Z) = A0 Anr n COSA, „2) + Bo Sinh(A,, , z) 
үз 1 
"> {ж Got? „L. as BE + œ? BE )+ тт, Loses Qs...) + BE B: )) 
T drr Je (A, оја Е J d drr J, (Am) f dz f: | “ш J 
——. 0 $ 0 0, С = 0 d 0 Ө, 
0,лг,пу 2 Lnr ng 2 
eub (2, 1) TN eub (27 3i T 
Or Na, Lara: St - | ofr V go Vnrnotr PU. y 
Ï. E 1, % 
2[агт [a0 550,0) 2[агт [a0 f;(r,0) 


2 


А = 0 0 
"^ OU" (a; B; Bras + BE BL.) 


2 


Ao, 








e. 0 0 
0,1" (o; B; + Bro; + B; Bil) 


T4 (r.0,z) = Aj (œ? Br (1, -z))+ (а + Béz)+ 


F Kë У. Anin (Conn, Zi (z) + Cu Zn EVA (Anr n осо is ] 


0 





ng-0 nr=l 
Solution du problème inhomogéne radial pour le secteur de cylindre plein: 


Le probléme inhomogéne radial est le suivant : 
CE a,T,(r,0,z)+ B,T(r,0,z) _ =Jf.(0,z) 


ait, (r,0,z)-BéT(r,0,z), —0  aÿT;(r,0,z)+ ВРТ, 0,2), =0 
%Т.(т,Ө,®)— B;T(r,0,2) =0 а; Т,(,0,2)+ B;T(r,0,z). =0 
Les fonctions propres axiales sont les suivantes : | 
Z,.(2)= à, Cos(4,z) + Po Sin(A;.z) 
A. valeurs propres solution de Sin(A I Xa. Ï aza: - B: Daf A; Cos; LYX(B;or; + Ba) 


n. nz z 


Lenz Cos.) Вто.) + (as зто) В; Cos.) 


1 
215: => (ai ti. Соз 1.) + BE Sin(A2.1.) Yo A; Sin; L) - BECos(A2.1.))+ a; BE 


nz z nz z nz z 


ZG = 





nz 
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Les fonctions propres angulaires sont les suivantes : 
Ө (Ө) = a? Ai Cos(? Ө)+ P Sin(a22 ө) 


ng 


Ав valeurs propres solution de Sin(2 Ө, (D J ata? - pe pe )= A Cos(4^ 6, YB?a? + Blat) 


ng 


MC (sx Cos(4? 6, )-- BE Sin(22 of + (a22? тле 6, )- B{Cos(22 ө, y) 
P0) "H, а costo, psu 0, ett sun o) gelt ape 
Le développement en série de la solution prend la forme générale suivante : 
Т(7,0,2) = Ka "S (4^, Cost, 0) + B^ Sin(4* 0) А: Cos(A;,z) + B; Sin(A; z)JI, (Xir) 

ng-0 nz=0 Ө 
Dans laquelle la partie radiales introduit un terme de division du fait de l'application de la 
condition aux limites inhomogènes en r: 


CL. a,T,(r,0,2)+ B,T(r,O,z) | = f.(0,2) 
Dy, ‚п2 = a, À, 1 ô (ENN E BI, (4,.1,) 


Ce qui donne Іа solution du probléme inhomogène radial: 
Z (2) = a, A; Cos(A;,. 2) + Bs Sin(A; 2) 


À. valeurs propres solution de 5їп(А 1 [E Jaja; - pip; )= A Cos( A; 1 X í Ug + B;a; ) 


nz nz z nz nz z 


VL (ал: Cat?) Bi Sins Ln Sim Best 
1 


2.0 => 
S 2 = lai. Cos(x 1.) + B; Sin(A,L) асл Sin) - BiCos(Æ.L)}+ o B; 





nz z nz z nz z 


6,,(0) =a? a2,Cos(22,0)+ B? sin(4*,0) 

M, valeurs propres solution de Sin(4*,6, Yos, fat a? - Bi B? P 4^, Cos(47,0, (Bla? + Bi at) 
‚16 ((a222,Cos(a2,0,)+ Be Sin(22,0,)f +(a2a2,Sin(22,0,)- в°соз(л,9,)) 

[ө®„@ 75. = (a2a2,Cos(a2,0,)+ 8 Sin(22,0, Na! 2°, Sin(a2,0,)- Bt Cos(4*,0, lu et Be 


n 


ө, 1, 
ст = | 4Ө®„„(Ө)| dz Z„ (z) f.(0,z) 
0 0 





bats улы M STEEN 
qoo H Babel ' Lë rr 3,1„ Al.) les (0) IZ, Œ) 








I| ne faut pas oublier les termes de valeur propre nulle, liés soit à la valeur propre nulle sur la 
composante axiale z, soit à celle sur la composante angulaire. Comme on développe un 


2 


raisonnement assez similaire au cas cartésien, une valeur propre axiale nulle An = D n'est possible 
que lorsque les conditions homogènes axiales sont de Neumann de part et d'autre et cela implique 
nécessairement une fonction propre axiale constante, soit 1. Dans ce cas les fonctions radiales ne 
sont plus valides pour représenter les solutions du problème aux limites, car soit elles s'annulent 
soit elles présentent une singularité : 

(Amr) =0 Ит. _ K (А.т) =œ 


Lim Ar ol Aj. A20 
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Cela coincident avec le fait que l'on a en quelque sorte réduit le probléme de la dimension axiale, 
avec un développement de la solution avec les mémes fonctions propres angulaires de la forme : 


T(r.0)= S (ио Br" )Ә, 0)- Y` 4/9 "0,5(0) 
ng-0 ng-0 


Les fonctions radiales et coefficients associés s'établissent par le respect des conditions aux 
limites : p (r) 2 r^ On doit méme envisager le cas de la valeur propre angulaire nulle, qui 
Fr 


conduit à une solution radiale en A+B Log(r), soit avec la condition de finitude, une constante 


Rao (r)=1 А; 


. Par exemple les coefficients du développement lorsque ^: 7 0 deviennent : 


O L 
fao | dz f,(0,2)6,,(0) 
CL. a,T(r,0,7)+B.T(r,0,7)_ = f.(0,2)=> A; q = 0 — 
- L |O, O) s BLAST 
Z o H 
Et les coefficients du développement pour le terme A+B Log(r), soit lorsque A, =0 et A i 
sont : 





2 








[aof dz f.(0,2) 


0 
9, L B, 


On rappelle également que le dernier cas A, =0 et Ang = 0 implique des conditions aux limites 
homogéne de Neumann de part et d'autre dans les deux dimension axiale et angulaire. On peut 
maintenant écrire la solution compléte dans le cas oü les valeurs propres nulles doivent étre 
envisagées, 





CL. a,T(r,0,2)+ B,T(r,0,z) _‚ =f.0,2)> 4 = 


Si Po =P; =0= Z, (z) = Соѕ(А,.2) À, = 





_ HZ T D oF = L 





@,9(9) ai el le BSin(2,6) 29, tq Sin(12,0, Jt ata? - pep? )=22,Cos(22,0, (Bat + Ba?) 
a, (la 222 ,Cos(22,0,)+ B? Sin(22,0,)f +(a2 4°, Sin(22,0,)- B!Cos(22,0, y) 
2 +7 (ar Cos(20 0,)+ B? Sin(2°,0, ot 29, sin(42,0,)- B? Cos(42,0, ))+ at p? 


1 
Je, y => 


@ E d 46 Ө (0) | dz Cos(A,z)f (0,2) А o= | dO j dz f.(0,2)0,,(0) 








ng—too À m © (0) 
T... (r,0, z) = ] | "0 o, а, ву Je, c» 
HN SR DER T 
L И ут Casel a, (А.Р) Ө „(Ө)Соз(А„2) 


Ach 0 0,0,1, UL) + BIS Q1) |0,,(6) 
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Et lorsque toutes les conditions homogénes sont de Neumann, de part et d'autre, la solution 
s'écrit : 


























Si EE Z.-CmULO AE |Z. (o)f =" 
@„9(9) =Cos(a2,0) A2, = [© (Ө) = 
"E -| d0 Cos(2°,0 j dz Cos(À, z) f. (0,z) 
L L 
X s | dO | dz f,(0,2)Cos(4,0) 45, = | d6 dz f.(0.2) 
0 0 0 0 
1 Au, X A, o 
T, (r.6,2) = Cos(4! + 
' 9L B, TA, LE т | "е r c 
4 ds C, z^ 48 (Ar) 
+ => Cos(42,0 JCos(A Ты? 
9, Je 2, Соз )Соз(А„г)— аа OAN pela z5) 


Solution du problème inhomogène angulaire pour le secteur de cylindre plein: 


Le problème inhomogène angulaire est le suivant : 
C.L. oft zt B,T(r,0, z) _ = 
aT. (r,0,z)— BIT(r.0,3)|,., --fé(rz)  afT.(r,0,z)+ BIT(r.0,3)|,., = f (r, z) 
aiT,(r,0,z)- BiT(r,0,2)) | =0  aiT,(r,0,z)+ B;T(r,0,z)| _ = 


La séparation des variables dans les trois coordonnées cylindriques donnent trois équations 

différentielles ordinaires : 

@T(r,0,2), 1 OT(r.0,z) | 1 &^T(r,0, z) | eT(r,0,z) 
or? r Or tp 00? дг? 


Dao tdo E far 90) (и Jomo 17€. ( Tac 
r 2 





-02 T(r,0, z) = R(r)G(0)Z(z) 











dr? r dr ? 
Les solutions radiales sont alors dans la catégorie des fonctions modifiées de Bessel d'ordre 
purement imaginaire, et dont la partie imaginaire de l'ordre est non entière : 


2 
Posons u= À et у= A. >r? - de m + (v? - ur? Ro) = 0 





Solution du type R(r)= A, J (i uu r)+B,Y, (iu r)- Ad; (u ce RK, Lu r) 


r IM Tv 


Fonctions de Bessel modifiées d'ordre imaginaire 


Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


Hélas sur une géométrie de cylindre plein, ces fonctions ne peuvent constituer une base de 
fonctions propres. Ce n'est pas tant en raison de la finitude de la solution, car ces deux fonctions 
restent bornées lorsque l'on s'approche de zéro (voire les propriétés des fonctions de Bessel de 
premiére et deuxiéme espéce et d'argument purement imaginaire, développées plus loin dans le 
document), mais surtout par leur caractére rapidement oscillatoire à l'approche de la singularité 
г=0, on obtient par exemple pour l et K : 

Développement autour de zéro de I, (u r) 





2k 2k 
: e (F) ue Ur 
z T |œ 1 z u r IV +оо 2 ivLog 5 ivLog(r) +оо 2 
L(z)- 1. meg ЗШЕ. 0. I N g MR un EN 
(=) (2) Sr) =н] | | Zen R < Triv De 
Développement autour de zéro de K, (u r) 


K,(z)= SS D r(k 4 + DA! (z) ` Е > г(к+ : + DA! où 
= Ки r)- s) > Tit < + DA! | eil (8) » r (k + > +1}! | T 




















2k 2k 
Qang Ur) dos Ur) 
—ivLog| S | _. + 2 ivLog| Z | ; = 2 
K. = 2 ivLog(r) _ 2 ivLog (r) 
= K, (ur) 2Sin(izv) < * атату ги Ё < ZT(k riv De 
ur 2k ur 2k 
ZA K, (u r) = | іл EET < (5) ES ele) puit) < (5) 
2Sinh(xv ) £xT(k iv +1)! £xXT(k— iv +1)! 


On constate une rapide oscillation qui méme si elle reste bornée rend impropre ces fonctions pour 
caractériser tout forme de solution comprenant la singularité г=0, a fortiori de fonctions propres 
radiales. I| est à noter que le caractère oscillatoire est identique pour les fonctions Bessel de 
premiére espéce et d'argument purement imaginaire. 


Pour le probléme général inhomogéne en 6, « c'est raté » ! sur un domaine bornée. Par contre sur 
des domaines non bornés radialement sur la base de cette construction par séparation des 
variables, on peut utiliser une transformation intégrale dítes de Kontorovich-Lebedev pour 
résoudre, par exemple, des problémes de Dirichlet de Neumann ou de Robin aux conditions 
inhomogénes sur des secteurs d'angles. Nous aborderons ce point dans un chapitre spécial dédié à 
ces transformations. 


En revanche pour un secteur d'angle cylindrique creux, soit entre deux rayons intérieur et 
extérieurs, l'utilisation des fonctions de Bessel modifiées d'ordre imaginaires 1 et K comme 
fonctions propres est possible car le systéme de Sturm-Liouville présente un caractére régulier et 
aboutit à une série de valeur propre discrète. Il en est de méme pour un secteur d'angle creux 
illimité au seul rayon intérieur. Dans ce cas, la fonction K à ordre imaginaire (dites de Mac-Donald) 
est la seule utilisable car finie à grande distance. Elle présente en outre des oscillations selon le 
graphe reporté selon l'ordre, signe numérique de la présence d'un spectre discret de valeur propre, 
qui confirme sa possible utilisation dans une série de Fourier généralisée. Nous aborderons un 
exemple de ce type dans la suite de l'ouvrage. 
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Ce qui ne veut pas dire que l'on trouve pas par séparation des variables des solutions particuliéres 
à des problémes particuliers inhomogénes en 0 plus simples. 


Illustrons cette recherche par une tactique d'évitement, qui permet de ramener un problème 
inhomogéne angulaire en la résolution d'une suite de probléme inhomogéne dans les deux autres 
dimensions. Cela permet également d'illustrer des solutions concrétes aux deux autres catégories 
d'inhomogénéité axiale et radiale. 


Prenons le probléme de Dirichlet suivant sur le secteur de cylindre plein : 
C.L. T(r,0,z). = Т(т,Ө,®)| o = TOS 2) =0 
T(r,0,z),, , = T(r,0,z),, ,, = Т, 
Il est facile de montrer дие la solution du probléme se décompose en deux sous problèmes chacun 
inhomogène axial et radial, comme suit : 
T(r,0,z2) 2 T, -T(r,0,z)-T,(r,0,z) 
T,(r,O,z) probléme inhomogène radial 
C.L. T,(r,O,z), , = Т, 
TF 02) Bé = T.(r,0,2)},., =T;(r,9,z) 





= T;(r,0, z) 





z-l, E 0 


z=0 

T,(r,O,z) probléme inhomogène axial 

CL. LEE = T,(r,0, z)|, , = Т,(т,Ө, z), =0 
Tir, 2). = T,(r,0,2)., = Т, 


Prenons maintenant la solution du problème inhomogène radial, comme énoncé précédemment : 


n, T 


Z. (z) = Sali zk si 


n, T 


= A tq Sin. L)=08 А = | n, =1,2,.... 


2 








2 


По 


6,,(0) = Sin(22,0)= SEA Xo tq Sin )=0S a n =1,2,.. 
0 


0 


à 

















l 2 0 

Z,@ =< => 
| nz (2) 2 {© (Ө) 2 

, 0 sin,oun, pair 

0 2 Т. 
с T, | dO 59029,9) dz ask ——(- C)" 1—(—1)” )= 4T,1,0, : 

"m as. 2 sinon 
0 0 (2n, +1)2n, +1) 7 
(2n, -1)z 
ng =+% n, =+00 Le. Ï r (2 1) (2 + 1) 
T (nb A 2 Y Е г Sin VIE |с Se * 77 g 
z^ Z £4 (n, n, +1) E +Ü) w; | i o 
Анор" 


РВЕ a ваш 
Ө, 
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Développons également la solution du probléme inhomogéne axial : 
6,,(8)- Sina? Ө) 4,14 Sinf.) =0 < 20, = s n, =12,... | 


0 


Ө 
0,009) => 





REI БЕ Жу 


l jp. es L, | i. -1 E L, | 


2 r 
к, „(9 = > ; 











г) А l 


S ` Sali at. 


Zom „y (2) = Sinha, „ (L. —)) 21, (2) = Sinh Qu 


nr ,ng 


KEE Wi drr J „ЛД, r) 




















= __21, zs Е 
Conr. = Суза, "d (ЖЕ "i an — C uda; =0 
fe 
А " _ 87, Jar r Jo. и) : 8T, ege | 
is men (2n, m 1) Ln DE DIR J (2n, i 1) Ln Us on al ) 
l, 
avec Jan = ү rJe Hour) 
T,(r,0,2)= > SAC Zn D CAC) OO 
ng=0 nr=1 
>T,(r,0,z)= 
(2n, +1) x 
ES 81, ng =+0 nr=+00 T cds aJ (a. 2ng D 0, Sinh, " E И ГА a (L _ z) 
12л ng-0 gc (2n, = a Á A pU y Sinh. ne ) 


Ce qui donne la solution suivante pour le problème de départ inhomogène angulaire : 
Gaz „ | 


nr,2n9+1 
0 


T(r,0,z) =T, — 


0 
Аы = 


r == = c2 r 
tq Te. GG. Ak 0 n, = 127243 T od E far r Je. (Aan, ar) 
0 








ў [= r] 
FEA 167, чу Ww 1 Ze L Sin (2n, + l) z z Sin (2n, SS Us 0 
TY É < (2n,+1)(2n, +1) , (6742, ) 

Ang: 1 2 


2 


. ( (Qn, +1) x 
BT, "enge )sin мыш Sinh(A,, an, 2) + Sinh A, a, (L. —2) 


L^z ng-0 nr-l (2n, + W ii (sciL ji Sinh(4;,. 2ng +1! 2 


On voit tout de suite que même dans les cas simples comme celui-ci, la solution est relativement 
complexe à exprimer. 


y 
Low ada E ТА 
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Secteur cylindrique creux 





La section est définie par le domaine produit cartésien d'intervalles des trois coordonnées du 
système cylindrique où typiquement on doit résoudre le problème suivant : 


©T(r,0,z) lOT(r,OO,) 1 G'T(r,jO,) OT(r,O,z) 
4 4 + 
or? r or r? 00° oz? 
O-ír/L,xrst,)x(09/0x0x0, «2z]x(íz/0xzxL) 


CL. aT, (r,8,z) - B,T(r.O,z) 4 77 fa (0.2) 
oft (7,0,2) - BaT (r.O,z)|, = — fo (r,2) 


aT. (r,0,z) - BiT(r.0.z) _ ——f,(r,0) 





=0 T(r,0,z) fini 


a T.(r.0,2)+ В„Т(т,Ө,2)] ‚ = 7200,2) 
a1 T, (r,0,z)+ BrT(r,0.z), , = fi (rz) 


aiT,(r,0,z)+ B;T(r,0,2)| = f7 (r,0) 

Probléme qui se développe selon la méthode classique de superposition en 6 sous--problémes. Le 
seul avantage de la recherche des solutions c'est que l'on a plus à se préoccuper de la symétrie 
angulaire des fonctions limites, car la partie angulaire de la solution n'a plus la contrainte de 


périodicité 2л. On va séparer l'étude selon l'axe d'inhomogénéité des conditions aux limites : axial 


z, radial z, angulaire Ө. Nous noterons les solutions des trois sous-problémes comme ceci : 
T(r,O,z) inhomogene radial 
Toin (r,0,z) 


inhomogène angulaire 
Т (7,0,2) inhomogène axial 


T(r,0,z) - T,,(r,0,2) - T4, (r,0,Z)+T „(r,0,z) 
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Solution du problème inhomogène axial pour le secteur de cylindre creux: 


Le probléme inhomogène axial est le suivant : 

X T(r.,0,z) | L'OT(r,0,2) | 1 O'T(r.O,z) , OT(r,0,z) 
or? TS Ór Up 00? oz? 

Q=fr/l rz1,]x(9/0x0x0, «2z]xlz/0xzxLj 


CL. aT, (r.8,z) - BT Qr 6,z) _ =0 aT, (r,0,z)+ B,:T(r,0,2) , 





=0 T(r,0,z) fini 


aT; (r,0,2)- BÎT(r,0,2), =0 ai T;(r.0,2)+ Be T(r.O.z), - 
a;T (r,0,z) - BéT(r,0,2) _, =-/f,(r,0)  ajT,(r,O,z)* BrT(r,0,z) _, = (7,0) 


Le développement en série est de Іа forme : 

ZER SU - (A5, Cos(22,0)+ B° Sin(2 0) JA; Cosh(A; 
T(r,0,z) = 
ng=0  nr-0 x (4 ng ,nz Ju, (Anang P) + Bu, VE nz ng (A, ‚лө r)) 

Les fonctions propres radiales se présentent sous la forme suivante, pour assurer la première 
condition aux limites radiale homogène de Robin, comme suit : 


Ror (r) = s nz Jo Ar „+ SC? nz А (4, nz ,ng r) 


a Кл, '(r) = BR, "ы Е, 


г) + В; Sinh(X, „ z) )x 


nr ,ng 


в (б?) Fa (Аш?) 

AA A en Bia ЖИИ AN RETZ NC 

йы, OR Boom il. 

mM rs i oua) бы» (а) As Yu Diesel? Ваў Dale) 
ü aAa Vol) m ВаЈә (һа) GA aal Ba (а) 


Et l'expression des normes des fonctions propres radiales (l'ordre des fonctions de Bessel est non 
entier)devient : 





= Rum f r)= 




















R (1 ) di (А, nr ,ng la) Yo sul „ng rl |) 
Ge a À nr Jg di Li, Jg la) Bad (Жл, 11) Q AA, Í, Jg a a. "E BY, (1 Jg rl |) 
R (1 ys Jy s, Jg 15) e (^, дајә) 
"ne Q4 nr А, аә 7a Jg L)- Bad y Q5, Jg la) Es ле d 7a Jg ln) - BI, Ace Jg la) 
R. Tora, Ка? '( ai = = рв (1 ) _ Jy Mo Jg 11) - Y, Qi, mln) 
- RDUM ад nr ,ng Jo 73 Jg 11)- Bad j ‚б ле li) аА, л, Ї, ле Da a ле L)- BY, ‚б! ,ng rl |) 
R. Rar А22) '( UR R (1 ) = PAC Jg la) o. Am ле 1) 
Aarm SE аА nr ,ng d a Jg La) ВЈ ‚(б min) E 76 ле X, 7А d) BY (Ж ,ng rl |) 
А.) A 
L| (p^ nr ng (LJ «(n nr ng (LP l= @ | Jp? nr ng (Lf + (r nr ,ng (LP 1- KI 
2 l, LG La (v, z) 
Rs D = 
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Les fonctions propres angulaires sont les suivantes : 
Ө (Ө) = a? 2? Cos(? Ө)+ P Sin(a2 ө) 


ng 


As. valeurs propres solution de Sin, DN 0 Jon) As, J ala? - pi p? E À Cosa? 2 9 ler ai + Blat) 


пө 


, Ө, Jun costa ө „)+ 82 Sal? 2,0 y (ota Sin(a° An, 9, 8,)- В, ^ Cos(4 Ө | | 
Io, => L (a84 Cos(4? 6, )« pe sin(z*, o, ai 3 glat ө,)- В?Соз(д° ө, baier 


пө mI 
пө 


Les fonctions axiales deviennent : 
aT. (r,0,z)— BiT(r,0,2) -0 


=0 


TEN ë " Zi, nr ng (z )= AA Jg Cosh(4;,. ng z) + Bi Sinh(A,, , Z ng z) 
aT, (r,0,z) * Dj T(r,0,z) Ka = f, me 


OT. (r,0,z)- BT(r,0,z) =- fq (r,0) 
| : DE S Zama, (2) = Au Col, (I =2))+ В.З", „ (L. — z)) 


L "nr,ng 





a;T, (r,0,z) B/T(r,0,z) = 0 


=, 
Et les termes intégraux du développement en série prennent la forme : 


r 


Аа = 
(ar, Совт. A ESCHER ln Gre 





|же Rua, (P) Í d0 fé(r,0)6,, (0) [ass Ru (0) dO fF (r.0)8,9(0) 


0,nr ng | 








Lnr „ng | 


2 2 2 2 
Run, (|Ә, (Ө) К, „(| 10,90) 
Terme А), ‚лө з d; (z) + Ci arn Zi nrn (GAS, . Ono (Ө) 


La solution de valeur propre nulle angulaire est constante et coïncide avec des conditions de 
Neumann homogènes de part et d'autre des parois latérales du secteur cylindrique : 


й„„=0 > ?=0 i | | | 
He Bi = В . Les fonctions radiales sont les fonctions de Bessel d'ordre nul, qui 
sont déjà intégrées dans le développement en série de la solution (cas ne=0). Il reste à envisager la 


simultanéité des valeurs propres nulles : A "CH eb AL, =0 et В, = В = Bo =Й =0 . Soit 
le cas de conditions homogènes radiales et angulaires de Neumann de part et d'autres. C'est un 
résultat similaire au cas cartésien puisque sont exigées des conditions de Neumann dans les 
dimensions homogènes du probléme. Il vient alors une solution purement linéaire en z , soit : 


Cas Pa В, Bi В’ 0 

















aiT.(r,0,2)- Tt = 70.0) f id | „у 
й > Ао = 2 
OU Ne E D Var pë + p 1 ai + B; BL.) 





aiT.(r,0,2)+ Tél = 
Z(z) = Aj „la + B; (L. -:) 


0, 


oT (,0,2)- BsT(. 0,2). =0 2 | ОСО 


2 
0,1 


, > 4; = 
aiT.(r,0,2)+ B;T(r,0, 2) _ = f; (r,0) 0, St "Le + Bras + Bo Bil L.) 


Z(z) = A; (ai + Bj z ) 
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La solution du probléme inhomogéne axial est donc, hors du cas de Neumann bilatéral homogène : 
0 0 0 0 H 0 0 0 0 
X 450,250, а EE Boos В, incluant lecas y mp0 et ор=о =l 


P o 8 o 
sauf f4,-B,-p,-p,-0 et о =a,, = a; -a, =1 























R ys Jo (Au, r) Е Y, aco r) 
T a, Ann I En 2 (E Jg La) T Bady (A En La) а TTE d лө А "A, лө La) js BY, o Qs, лө La) 
А а, M „ng Jo (EM e 1.) t Bia (4, лө La) a, D ND d En А (EM Jg 1,5) + B; Y, o Qo, Jg 1.) 
Ж tq PY r ш r (зг r 
a „Аз ‚лө ы (Ar лө La) = Bad (A, 4 лө л) а ТТ En D (4 ‚лө 11) y BY Und ag rl |) 
(1 ) PAC лө La) Y'A лө La) 
P SE ЖИР ‚пө Ае a En 11) DE Bad (Ar Jg ln) Q AA, „Ў, лө po al En ln) gr BY, , Las лө a) 
2, Us ‚лө L5) Y, (А, Jg 1,5) 
D 19 ( l.) = 





A AA, wl лө А, al ‚лө 11) a BJ (A, s лө a) а, Va n, T; En А, a io Li) 7 BY, (2, ‚лө La) 


ln. 04: осе "ln. 0.) + EE x. 





ms dis: - 
6,,(8) =al% Cos(4* Ө)+ Bt sin(4? Ө) 


пө 








2 


A^ valeurs propres solution de Sin(2?. ө Di Fala! – B? p? È 29 Cos(4* 0, (Blat + Blat) 
Ө daks Соз(4,,0,)+ BE Sin(22.0,)F + (222 sin(4* 0,)- Bt Cos(42,6, y) 

ds I (a222 Cos(4^ ө,)+ B Sin(22.0, Nata, gä a, LgrColät 0,))+ eg 
ды O) = aj pelt. (Le rimi, A 

Z wn, (Z) = 904, n COSR(A,, n, 2) + Bo Sinh An n, 2) 





r 


"SUM x sa, COSS, nl 1.)(ai 8; ta; p; )+ Sinh(A;,, L, ea (a; De: + B Br ) 


% 


| drr Ru „(| de f$ (,0)0,0(6) [TP 








0,nr ,ng 





L,nr „ng 








Ru, (Ol le 





nO (ө)? 











Р, „Ө, „Ө 


T(r, 0, z) = 3 У EN ‚лө (С C, ,nr „ng ZA ens (z) + Саад, Zin Ron (r)O,, (Ө) 


ng-0 nr=0 
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La solution du probléme inhomogéne axial est, dans le cas de Neumann bilatéral homogène : 














| 8 P 9 8 
Si B,-B,-p,-p,-0 et a,-2a,,-a, =a; zl 














= (CE ng r) Y, (4 a лө r) 2 ‚лө lya) 6 ‚лө 12) 

R. n (r) = r рү À, пө tq "TZ Ve at 
CS 2, (і ле SUM A Dar по rl 1) MC ng La) o ` (An пе 11) 
2, ur ng 11) Ү (4, ле 11) (А ле L3) Ү (4, ng 1,5) 

Ro (la) = Шр ' r К. ( 122 Ve at "Zr 
MCA jg Sg Y А, (A, Jg ly ) MG o SN Y Dé (4, ‚пе ly ) 


D „ng (Li) = 0 Ds Jg (L5) = 0 


A CA) 
pete c x rl (к nr t "n » Ay 
l.a EG ) ЖОЕ, F ‚лө 


R Ge 
wa, O| > 














et |Ө,„00)] - 





Ө,(0) = Cos(420) et A tq sin(40,)-0— 2° = ; 
0 


Zo, nr ,ng (2) = a, А ‚лө Cosh(4;,.,. U пө (1, - 2 г))+ В; Ѕіп(д: nr ,ng (L = z) 
Z wa (2) = ASA, „ COSh(A,. „ 2) + By Sinh(A,. „n, Z) 


nr ,ng nr ,ng 
r 


A 77 = | А 
( / Cosh(,, vs Yo B; +a By )+ Sinh(A,, onis ia: (ai nr m, + B; Br ) 


nr ,ng 





4 | drr R„„ (r) | dz Ду (т, nod пт " +) 


Со, = 
,nr ,ng 2° P À, 
OI Kean viir 2a) 


1.5 Do 
4|drrR r)|dz f;(r,0 RS 
| s 0| fÀ E | 














Larang ` 
B ER ER 
2 _ na ane _ ii; 
0, Lo De c p ES y La (Run, Ok 1 js ES E у 
а fao reo АС 
40 = ў 


DIS Ss ger. + Bras + Bo Bil 1.) ^u = a? = Tape + Brai BERI .) 


T4 (6,2) = Aj ler Br (1, — z))+ Aj (е + Biz) 


Р ` Y rom Ж DN (2) FC LL авы OR. (Ca al 


ng-0 пг=1 
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Solution du problème inhomogène radial pour le secteur de cylindre creux: 


Le probléme inhomogène radial est le suivant : 
O'T(r,O,z) lOT(r,O,z) 1 O'T(r,O,z) O'"T(r,O,z) 
+ + + 
or? r or r? дө? oz? 

О = {5/1 &r€l]x(9/0x0x08, <2л\х{2/0<25<1.} 


CL. a, T.(r,0,z)- Bal (ne). =-Jf„(0,2) a, T. (r,0, z) В.Т(.Ө,2) = f.,(0,z) 





=0 T(r,0,z) fini 


asT,(r,0,z)- Ba T (r.O,z), , -0 ai T, (r,O,z) + ВРТ ,Ө,2) о zi 
aiT.(r.0,2)-B;T(r.0,3)) ,=0  ajT.(r,0,z)+ B;T(r,0,z)) , =0 


que l'on sépare comme d'habitude en deux sous-problèmes : 
CL. a„T,(r.0,z)-B„T(r.0,z) 470 a„T,(r,0,2)+ B„T(r,0, z) Late) 


at, (r,8,z)- BS T(r.0,3)), , =0 ai T;(r.0,2)+ В'Т(,Ө,2) ER 
aiT.(r.0,2)-B;T(r.O,3)) ,=0  ajT.(r,0,z)+ BrT(r. 0.2) =0 


et 
CL. a„T,(r.0,z)-B„T(r,0,3)| , =-f„(0,2)  œ„T,(r.0,2)+ B„T(r,0,z)) __=0 


asT,(r,0,z)- BST(r,0,z)) sc af T,(r,0,z)+ BI T(r.0, 3), =0 
aiT,(r,9,z)- В;Т(,0,2) =0  afT.(r,0,z2)+ B;T(r.0,z)| =0 
Les fonctions propres axiales sont les suivantes : 
Z (2) = аА, Cos(A; z) + Dy Sin(2;.z) 
AL valeurs propres solution de sin LX. J aiai - B; B; )= A; Cos(Ai.L. \Bai + Bei 
Гексан.) Bii LJ + (ai Breet 
2.0) = 





1 
2| + e Cost 1) Bi Sina LO Nag SinCALL) — Bi Cos) eof 


Les fonctions propres angulaires sont les suivantes : 
6,,(8) = a£ Cos(22 0)+ Be Sinf? ө) 


Mg 


Ав valeurs propres solution de Sin(2 Ө, We, Jala? - pe pe )= м Cos(4^ 6, YB?a? + Blat) 


пө 


16 (аел Cos(4* 0, )+ B?'sin(4* 0, )) + (020 Sin(22 0, )- 8@Со5(4% ө, y) 
le, Gy =, 1 


= (a222 Cos(22.0,)+ B Sin(22 0, at? Sin(22.0,)- B2Cos(22 0, )+ a? BE 


ng 


Le développement en série de la solution prend la forme générale suivante : 

Р b» п (A? Cos (49 0) - Bo, Sin(49 0) AZ Cos(4;,z) - B; Sin(4;, z))x 

Т(7,0, 2) = íi i 
e x4 ue carey 


ng,nz^ ng по ,nz 
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Dans laquelle les parties radiales pour les deux sous-problémes introduisent un terme de division 
du fait de l'application successive des conditions aux limites inhomogénes en r: 

CL. a,,7,(r,0,2)* B. T(r.0,z) 4170 

m De (Ar) Ky (Ar) 

SUP Ga, (ls) * Biss (4,1, ) (Ек, Al) d B.K TUN) 

CL. oT, (r,0,z)- B„T(r,0, z) _ =0 


=> R 








72 Lo (Aer) К (Ar) 
> R,, „„ (r) = 2 1192 ; 2 2 1/92 : 2 
А.а, 11 (4) I Bal, (4) Аа K 39 (Ailn) = В.К (2,1,1) 
CL. a„T;(r,0,2)-B„T(r.0,2)) =—f(0,2) 





I, (ur) K o (ur) 
A nz '(r) = А. 2 ' 2 = 2 2 ' (27 - 2 
| WEEN (2,1,5 ) F BJ. (A5) defi (Anal a ) + BK, (2,.1,) 





rl аА. (rl) ж; Bal (Qa) à „An K о Al) ES В.К (2,1,1) 
v^ Rate ОД OL). (оК. Clu) + В.К G) 


CL. a,T.(r,0,2)+B,,T(r,0,2)… = f,,(0,2) 





Rtl A Me К Lei 
WE ш Gas (An Li) ES Bal TAR) Lea (la) ш В.К Ali ) 
2 ol, (ola) T D.I, (A715) а,,А,.К „ (Q7, L5) * BK, (2,.1,) 
D'E == 10 - = = 


Y À; a Lo "G-A = Bal (4,1,1) A a Ко (Orla) 7 В.К Al) 


nz ^rl nz ^ rl 
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Ce qui donne la solution du probléme inhomogène radial: 
Z (2) = асл, Cos(A;,z) + Bg Sin(A; 2) 


À;, valeurs propres solution de 5їп(А 1 (x. одо; - pipe )= A Cos( A; 1 (Во; + B;a;) 


nz 2 nz 2 


„асл Сов.) + BE Sin LD) + (оса sin) - 5 Cos) ) 


zop- i 
S 2| (ос Cos (45.1) B; Sin(A;.1.) асл: Sina L)- B; Cos (4;,1.)) ai B; 





nz z nz z 
nz 


6,,(8) =al22,Cos(22,0)+ Be Sin(22,0) 
M, valeurs propres solution de Sin(22,0, DÉI J ala? - pe p )= a? ,Cos(22 „Ө Le ai + В; Ва?) 


0,((a222,Cos(22,6,)+ B sin(a2,0,)f +(a2a2,Sin(22,0,)- B2Cos(22,0,)f) 








2 = — 
10,90) "D + = ag 22,Cos(42,9,)+ 0 Sin(2,0, lei 22, Sin(22,0,)- B Cost.) op 
E Do (A; „r) К (Ar) 
Кы Ga, L'U La) F ВЫТ TAR) tee. (ola) + В.К GL) 
T Ly (Ar) К (А.т) 
Gei B (a A. QA Lo, (A OR Bal. Zt) ae, ав Pak, (Жы) 


ge -| 406,,(0) | dz Z, (2) f„(0,2) c2, = |409, (0)ја 7„(г)/.(Ө,) 





T (r, 9, z) = m 53 S e nz К, nz (r) „ Erg Км (r) О (0)Z,. (z) 
4 AHL D D. JO OPIZ 


Il ne faut pas oublier les termes de valeur propre nulle, liés soit la valeur propre nulle sur la 
composante axiale z, soit à celle sur la composante angulaire. Comme on développe un 





2 


raisonnement assez similaire au cas cartésien, une valeur propre axiale nulle As. = Gë n'est possible 
que lorsque les conditions homogènes axiales sont de Neumann de part et d'autre et cela implique 
nécessairement une fonction propre axiale constante, soit 1. Dans ce cas les fonctions radiales ne 
sont plus valides pour représenter les solutions du problème aux limites, car soit elles s'annulent 
soit elles présentent une singularité : 


A, 0 Ds (4, „r) =0 
Lim, Kyo (Ar) =œ 


Aan 


Lim 


Cela coïncident avec le fait que l'on a en quelque sorte réduit le problème de la dimension axiale, 
avec un développement de la solution avec les mêmes fonctions propres angulaires de la forme : 

ng =+œ0 m E 

T: n p 4% 

T, (r,0) = > | „oY "^ + B,gr TMO 

ng-l 
Les fonctions radiales et coefficients associés s'établissent par le respect des conditions aux 
limites : 
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Aso =й 
" 
ш 79 R. (r) = (2, а, + ur) + (45, o, STE ) 
rl rl 


Xe -Ao 
r r 
kel = 0 R (r) = (2, d „o + ыбы]; Е +(4, a, As |; | 
r2 r2 


On doit même envisager le cas de la valeur propre angulaire nulle, qui conduit à une solution radiale en A+B 
Log(r) : 


y 
Cas(l a«,4Rj,(r)- В.К, (c) =0 > R... 9 (F) =а, + вло т) 
rl 


Cas (l) «4R,(r)-B,R 









































' 
r2^'rl r2^ rl 


l 
Cas (2) „Кв (P) + B,2R,10 GI. = 0 — Rao (r) = aü, + LaB,Log( =] 


2 


Par exemple les coefficients du développement lorsque À = 0 deviennent : 


Cas CL a„T,(r,0,2)+B„T(r.0,2)| _ = 300,2) 





2% -Aho 
9 9 [5 | % 1, 
(4459, 14f ES di + PAC epa Ang Qa Laba] É RE [def dz f.,(0,2)0, (0) 
, A"? к 1 0 0 


rl 
: 4| ^ 1 18,02, 
dë дей 
rl 


CasQ) СІ. a„T,(r,0,2)-B„T(r.0,2)| | =-f.(0,2) 








l 
+ Bats (2, Qa zk Al = 


rl 


6 


2% ag 
ө o Ly ` er 
i s (22. @„ DOE (2а а ЕЯ E = D + | | аө | а f4(0.2)0,,(0) 
_ 0 0 





1 r 
=> D, Ed 0 ө Ао = 2 
n mu Ko [өө D, 
+ Bata (2, C» + Lp n +4 (22, a no © mm Â 
r2 r2 


Et les coefficients du développement pour le terme A+B Log(r), soit lorsque ^^ ` 


Саз) CL a,T(r,0,z)+ B,:T(r,0,2) , = f2(6,2) 


$ L 

Lf , , faofa 6.0.2 

D,20 = 9, Nu F ва, + „бло 7 J > 4/2 = 33 9-70 
rl 0 fz 


r2 








D 
CasQ) СІ. a„T,(r,0,2)-B„T(r.0,2)| | -— fu (0,2) 


9, L 
и [def dz f.(@,z) 
> rl | 0 0 
| ER? 








С р 


rl 


Ï Ï 
Do =d Е + 16 + „Ва = 


r 1,0 
7 А XV. =0 et sf. . Y. gek 
On rappelle également que le dernier cas `” ne implique des conditions aux limites 
homogène de Neumann de part et d'autre dans les deux dimension axiale et angulaire. On peut maintenant 
écrire la solution complète dans le cas où les valeurs propres nulles doivent être envisagées , 
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nz T 


Si D = В; =0=> Z„ (z) = Со5(А„2) À, = j 











l 
ZG) = E 


6,6) - att, Cos(2^,0)« B'sin(2,0) 29, а Sint, Jat atat — BOB? )- 25, cos,6, (Beat + Bla!) 
Ө, (atat, cl al В8@5їл(д°®ө,))` + (e 29, sin(4,0,)- B'Coeläi éi y) 








1 
OOl == 
Bt", E (at 22, Cos(42,0, )+ Bt Sin(22,0, Na! 2° „Sin(22,0,)- BECos(a2,0,))+ a pe 
P nz (r) =T “a c Z z ы 2 
» Gate (la) + Brat jo, (1) (аск, (Qa) + В.К Al ) 
Le, (An d К», (An r) 





RE (r) = { 2 ' 2 T ) ( z ' 2 2 
Anz X Al, (4,1,1) z Palo, (4,1,1) s OAK о (Qa) E BK o, Qu) 


rl Gs po, (la) ш Bal, (4,1,1) GAS Ко (ol) tw BK y, (421,1) 


р, nz 2 ' 2 2 m z ' 2 2 
S An Qul (4,1,2) F Buol po, (42.1, ) ae, (A15) + BK (Al) )) 
a, 1 (ola) + Biol, G-A) Anh K S (Al) + B.K jo, (Ql) 





nz 19 
r2 > s Ano 


dd Аа, (Ql) = Bal, Al: H (алк, (al) = ВК о Ali )) 





по ,nz ng nz 


% L 9, 1. 
el. = | 40©„(Ө)| dz Cos(A,.2)f,(0.2). c?, = [d00,,(0)[ dz Cos(2,.2)f.,(0,7) 
0 0 0 0 


Ano -229 
r r r 
Ет) = (4%, a, +18,» E A (4%, quel li 
r2 r2 
e e 
P r r 
B (r) = (26, а, + "a + A, CAE "a 
rl rl 
A -483 
o, Aal Ao 0,5 +18 а {е =] В ы + 
r2F r2 1 nO ri r2F r2 1 d ; 
1 r2 r2 0 2 
Digo =T " „ү | Alo = [aof dz 7.00,2)Ә,„(9) 
rl 1 nO 1 —À,0 ^ 9 
To DN c FL D (45, d „o Lis [E 
r2 r2 


д9 E 
1 пө 1 nO 
A, ong D ара fe — (2, di DD E 


rl 


0 L 
А = | 40| dz 7,200,2)Ә,,(0) 
0 0 


0 =й 
1 пб 1 n0 
TE? (5 а, + а) + E di - la B D 
rl 


n Dj, Jee 


ng,0 





ng-l 


NC = 
rih N 97> PET 
L E SS s eR n) E c? R7 „ (r) \ Ө(Ө)Соз(А„) 

р" D? Io, CIR 


ng-0 nz-l по ,nz пө ,nz 
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Et lorsque toutes les conditions homogénes sont de Neumann, de part et d'autre : 


Si EE =0=> Z„(a)=Cos(Kz) EEE |Z, (2) == 


2 




















ng T | 


6,,(0) = Cos(42,0) À°, = А 





0 
e, GJ => 








Lo (Ar) E, (Ar) 
RE (r) = E ' 2 s 2 " Z ' 2 Т 2 
at, IEN + Bl, (Als ) аск (Al) x BK, Al ) 
De (Ar) Ku (Ar) 





RT far DB OLEO] aX 01 dE BK Ub) 
PARIES HU), GR TE RG) 
wm rate Gba) Bale 1) Жу Cala) EE Can) 
2 VISE US) GR НОИ) 


We Аа (Ala) Е Palo, TAA) (алк, (WM Em ВК о GA, )) 








Dg ,nz 


cu fao о,ө) dz Соѕ(А,2) (0,2) с? = fao Gala dz Cos(À, 2) f. .(0,z) 
0 0 0 0 


Déi —А% 

Е r \ ` r В 1, 

Ri o (r) = (45, „+ ih. |; | E (2, a, AE | A (r) = а, + Lä " | 
r2 r2 


Ano Ato 
Е r r Е r 
m (r) = (22, Qa Bs LiB, (=) + (X. Qa us LiB, |5) Ri (r) = QA * ibtd | 


La La rl 
0 0 5. is 0 La SCH 
Oe (2, Qo + Las "m n (4%, 9,2 Е LÉ T + à j 
E 1 r2 r2 . A d 
D, F L. р i ADS - [ao] dz f (0, z)Cos(22,0) 
r l nO 1 no 0 0 
+ Bla Ge а, +1,,В, D * DN e cs (2) 
r2 r2 
| Ato | -Aho 
UA ES ал tL, E E EN C4 JD + MP 
E 1 ri rl А 0 2 
m = L. ji E А» = [40| dz f. (6, z)Cos(4,0) 
r l n0 l пө 0 0 
+ 8,1, D Out а) + ES а, - Li, D 
rl rl 
LB | % L 
Ds = di Im + Ba e. + EE a) 4, = [do] dz ЛА (Ө, 2) 
rl rl 0 0 


E 1 1 А 90 l 
D = c Sei T B, H p Т) А = [do] dz fa (Ө, 2) 
rl 0 0 


r2 
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z oz (00 _ рө _ 
Ce qui donne la solution complète du cas : Bo = Bi = By = B, =0 ` 
1 (Ais Ros | Asso ), 





T.,(r,0,z)= 
didi c Э 
2 по =+% А"! СТА r NS n. r 
d Y cos(a2,9) ^ 2 o( ) 0 s ol ) 
1.9, пө=1 “рт 0 D,,0 
Rm 





NIERO уз ri rl r2 r2 

У у C Gi 0 (À ) Cno Кж (r) T Суу: Pong nz (r) 
ER OS, 9 OS VA „Z zi 2 

Oz | nat nz-l no ,nz ng nz 
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Solution du problème inhomogène angulaire pour le secteur de cylindre creux: 


Le probléme inhomogène angulaire est le suivant : 
©T(r,0,2) 10T(r,O,z) 1 O°T(r,0,z) OT(r,0,z) 
+ + + 
or? r Ór r? 00? oz? 
Q=fr/1,<r<l,}x{0/0<0<08,<2r}x{z/0<z<1} 


CL. a„T,(r.0,z2)-pB„T(r.0,z) „ =0 a„ T. (r,0,z)* B„T(r,0, z) „ =0 





=0 T(r,0,z) fini 


asT,(r.0,z)- Biet ai T,(r,0,2)+ BPT(r,0,2),, = 2) 
aiT,(r,0,z)- B;TQr, 0,2) =0  aÿT.(r,0,z2)+B;T(r,0,z). =0 
Problème qui se décompose en deux sous-problèmes : | 
CL. a, T.(r,0,z)- В„Т(т,Ө,2)] __ -0  a,T(r,z)- ËU = 0 
%Т,(т›Ө,2)— ByTQ.O,z), —0  afT;(r,0,z)+ BrT(G. sz), = ff (r,z) 
aiT,(r,0,z)- BT(r,0,7) =0  afT.(r,0,z)+ BiTQ.0.2) | =0 
еї 
CL. aT, (r,0,z)- Sitë, -0  a„T.(r,0,2)+ BT(r.0,3)) =0 
agT,(r,0,z)- BÈT(r.0,3)) Gei — af T,(r,0,2)+ ВРТ(,Ө,2) =0 
aiT,(r,0,z)- BT(r,0,7) | =0 аТ,(,0,2)+ B/T(r,O,z) | =0 
Les fonctions propres axiales sont les suivantes : | 
Z (Z) =a; A, Cos(X, z) + B; Sin(A, z) 
A valeurs propres solution de Snad aza; - ps p;)- д; Соѕ(А; L X оо + Por ) 


LL (asa cos: 1) + B; Sin(A,L)) + (ос Sin(A.L)- B; Cos(AL.I y) 


nz z nz 2 nz z 


D => (aj Ai, Cos (A; 1.) + BE Sin(A2.1.) Yo A; Sin(A.L)- BECos(A2.1.))+ a B 


nz 2 nz z nz 2 


20) = 





nz 


Comme a; A; (or; A; Sin(A.1.)- B; Cos(4;.1.)) Br (eA; Cos(Æ.L.)+ Sin(A;,L)B; ) 


2 = l| (zat. Cost.) B; Sins Lf | 1, + BE en, PC |+ ai Bz 
2 ai(x.) a; 








= 





ои 











Z; = в CL) SOL i (aJ S | 
1 1 
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Pour construire les fonctions propres radiales, revenons à la séparation des variables dans 
l'équation de Laplace en coordonnées cylindriques. La séparation des variables dans les trois 
coordonnées cylindriques donnent trois équations différentielles ordinaires : 


OT(r,0,z) 1OT(r,0,z) 1 O'T(r,O,z) ӘТ(Р,0, 2) 
2 + + 2 2 + 2 
Or r Or r 00 02 


нота (i) -EJ jun -e 90) (и foomo ZO (y J ze)-o 





=0 = T(r,0, z) = R(r)G(0)Z(z) 











dr? r dr r? do? z^ 
Les solutions radiales sont alors dans la catégorie des fonctions modifiées de Bessel d'ordre 
purement imaginaire, et dont la partie imaginaire de l'ordre est non entiére : 





2 
Posons u-2AÀ, et v=%, =r? 2 tr s + (v? — ur? )R(r.) =0 


Solution du type  R(r) = AJ, (iu r)+B Y. (iu r)=A I. (u r)+B К „(и r) 


riv r< iv кету 


Fonctions de Bessel modifiées d'ordre imaginaire 


Relation d'orthogonalité des fonctions propres 


La séparation des variables dans les trois coordonnées cylindriques donnent trois équations 
différentielles ordinaires : 


Orthogonalité des fonctions propres ÍR, (r), A, =X) ÍR (r) А, =} méme A. 


nr nr 


r| ERO, L eo. [69 Je "ei -o 





dr? r dr r 


s a. 1 = UE 


dr? r dr 


E Z2 AO _ в) e (x cl a) = 


r dr dr r 











> {r ARC) 8 (ys M =d DEER | ел = 


dr dr r 





Les conditions aux limites sur r, et r, => LE AC) -R (r) c all =0 
r r 


=> Si X A, alors ja BIRI о 
r 


ñ 


En d'autres termes il faut utiliser un poids 1/r pour démontrer la relation d'orthogonalité. 
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Fonctions de Bessel modifiées d'ordre purement imaginaire 


Posons-nous la question : est-ce que ces fonctions de Bessel modifiées sont à valeur réelle ou 
imaginaire, pour cela faisons appel à la propriété du miroir des fonctions de Bessel : 
Vv eC, set, I, (z)- 1,(z) et VveC,zeC, K,(z)- K,(z) 
D'autre part pour K: Vv eC, z eC, K ,(z)- K,(z) 
EtenfinvveC,zeC K,(z)=—7— ,(z)-L(z)) Sin(inv) =iSinh(xv) 
2Sin(zv) 
Il vient 
T T 
VveC,xeR K (х) = ———— Uart) L) == = = —— M j (X) Tj (X 
el sec Us (0-1, Ge ZAI )- 1,9) 
VveR,xeR K | (x)- К, x =K, „(x)= K, (x) à valeur réelle 
VveR,xeR І „(х)= 1, (x), 1, (x) à valeur complexe, comme suit 
T T 

K (x) (1 (x)- 1, x)=- 

2 em 4 2 ) Sinh(zv) 
On choisit donc d'utiliser les fonctions : 
K, (x)= K, (x)= Re(K,(x)) et Lis: Re. (x) 
T.M.DUNSTER dans un article « BESSEL FUNCTIONS OF PURELY IMAGINARY ORDER, WITH AN 
APPLICATION TO SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS HAVING A LARGE 
PARAMETER » introduit une autre fonction : 


Im(1 (z)) = I, (x) à valeur complexe 








+ T 
VveC,zeC L, (z)- Sinav) (1, (z)+ I. (z)) 
Ti T 
=VveR,xeR All ZL se) (en (+ La(s) 
T лі z л 
= LE se) Lo c Rt, (5) 


=T, (x)= веб, (x) = DL, (s) 


IV 


Le Wronskien des deux fonctions | et K est : 
1 
w(K, Gr (= À 


„су 


MR CITE (a), L, Gr) e Uc Dal All o — 





GA w(&, (x), 7 (x))= Sinh(zv) іл 1 


v 


л Sin(irv) x =x | 

Les formules de connexion en changeant le signe du paramètres sont les suivantes : 
K_„(x)=K„(x) L_„(x)=-L„(x) Vu 

=> K ,(x)=K,(x) T. (x) (ә) 

Les fonctions I et K d'ordre purement imaginaire ont toutes deux une infinité de zéros simple entre 
[O, v] et aucun zéro pour xèv. 
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K (x) K(&.)-0/Ak, Fir у> у>, >k,z>...>k,„>...>0 Lim, „„k,, =0 


V.S Js-] 5—00 "V,S 
TET f. J- iu cs ү ta >t, ig >...>t,„>...>O0 LiM pois 0 
De plus les zéros en I et К s'intercalent comme suit : 
K (k,  )=0/1k, + Li =0 
Ein, Belli ЫА АЮ | ch 
ARE | Lim, „„k,, =0 


La formule des dérivées pour les fonctions І et K d'ordre purement imaginaire : 
K,(x)- K,, (x) 


K,G)- к, (к), (к) X,'(s)=ZK(s)-K, della, (а) К, (8) 


> K,G)- к(а) К, (х) Ell кк) [К ele к(а) 


К.) = K (х)= К, (к Эко Jk Je, sl HK Jee, All 


Kj, Ax)- K x = Katy ) 
K,'(x) = - Re(K,,(x)) = - Re(&,, (x) 
1, (x)= Re(1, (x) 


POO dl 16) д (а), (а), (ж, (е 


2 
rx )- jl. (x )+ LN ES 1-7, (х )= Ке(1,, "(х )= 584, (x al )] 
Tus (x) = TO ë LZ) zd uu) 
ы х) „(= ке „(+ 1, All 


De plus Rell, (= Ref, i) et Re, (x]= кеј, ll 
Belt. (x) (x)]= ndr a) et ка (= RelZ 05) 


=> (к) = pl, elei, del: Rd Ge Lua) 


FON 
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Le caractére oscillatoire des deux fonctions Bessel modifiées d'argument purement imaginaire est 


confirmée par les formules asymptotiques suivante autour de x=0 : 
Ps = Агв\Г(1@+5+1у)} o, = Arg (Le iv)j 


K,(x)- ss t Log| | EN d об) 
Т (x)= sees) (eo 5) - e. | «ote ) 


L'amplitude des oscillations n'est plus bornée lorsque l'argument tend vers О. 


Au valeur x->c», la fonction K tend vers O, mais la fonction L et donc | diverge en exponentielle. 


soJe) 
о RUD 


Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 


Fonctions de Bessel d'ordre purement imaginaire 

Toujours dans l'article « BESSEL FUNCTIONS OF PURELY IMAGINARY ORDER, WITH AN APPLICATION 
TO SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS HAVING A LARGE PARAMETER » 
T. M.DUNSTER introduit également deux autres fonctions de Bessel : 


vreGzeC Ell — — (UG), 6) Gil 
aco "5 | 


>F ,(z)=F,(z) G,(z)-G(z) De plus Y,(z)= 


v 

















1 1 
нт) (Coste, = (Este 0-4.) 
Cos| 7” Co FY 
bau Dau a ennen dtt 
8 ORAO ) C). кє, (s) "t 
ES 2 LM sgle 2 == _ Ке „= 
acos Z2) acos Z. SE 
_ (Cor). 0), (Coste) LA) 
Lett Ѕіп(лу) R 25іп(лу) 


кл та E Pas 


Ѕіп(лу) 





Si =) od? ad? 


Cos| 7” 
"C Gr fit ZO- _ Y (#)=Y. kz) 


(J a(z) +J (z))= в, (к) Yale) X6). X G)- Y. (2) 


2Sinh| ZY 2 ï Sinh| Z 
2 2 
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Ces fonctions sont donc directement proportionnelles aux parties réelles des fonctions de Bessel J 
et Ү d'ordre purement imaginaire. On peut donc prendre comme choix de fonctions également les 
deux parties réelles des fonction J et Y : 


(r) 6,69 > 7,ь)=кеб (2) e E в, 0) 


2 


Е, 


iv 


G (= RE 7 (= relr, (x )= cos E be 
лу 2 
SC 
2 
Le Wronskien des deux fonctions С et F est le suivant : 


(к) 6,6) Z 


SEAN 
____\2/ 


=) cos Z rti e. (= — 


Les deux fonctions F et С affirment encore autour de la singularité x=0, un caractère oscillatoire : 
P, = Argo seiv)) o. = АГ) 


211 ani Z 


dei ol, SR = о.) «oz ) 


1 
2 
2Cotanh 7) 


GG» | — (suv nelle. kel 


L'amplitude des oscillations de F tend vers 1, celle de G tend +œ, vers lorsque l'argument tend vers 
О. Au valeur x-»ee, la fonction F et G tendent toutes deux vers 0, en oscillant d'amplitude 
décroissante. 


F. (x)=- 





Les fonctions F et С ont toutes deux une infinité de zéros simples respectivement entre [0, т] et 
[t >] : où la constante t est défini comme suit : 

т=т,—1 où T, racine positive de l'équation Cosh(r, )= Tg 

tT = 0.6627... 


bep VT > Дн > PAL > SE >..> Де >..>0 Мт, flore] =0 


F (x) F, (f...)=0/ 


а [vr, +0] [ve 399] [vr,+0] [vr, +0] . [ve] _ 
TA VTS RE "ego eque c <..<о Eu. er =%® 


Va s=1 


| los vr > girl» lvl lvl ail 20 Lim sl 
G, (x) G. (g, ,)=0/ M anm 
gd" оо ea с ж к = 


On note également que le système des zéros des deux fonctions F et G est entrelacé. 
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En revenant aux notations du problème les fonctions propres radiales sont une combinaison 
linéaire de la forme : : 


2 r D 2 
Ron) = An, „rI AF An LZ r )+ B AE K A nz r) 


Les fonctions radiales doivent respecter les conditions aux limites homogènes radiales sur la 
section creuse : 


GE, гав.) BR). = 

T, Qn) Ky (Ar) 
auis VE S VIENI ЕЛЕЕ Л NUES 
CE RSV) BR (01: 


> К, „(7 г) = 





~ 


Y a, 24% а | (д1 о) + "n A. EN (OAK, r Läd) E В.К. (42.1, ) 
— ^m, tq z Toa z z Y ï 
j a, 51 Dé (A l) В, T (A4) | aAR, r LA l) - В.К (AL 2 1) 


Comme А, dépend de A. on note А, > А, 


nr,nz 





On voit que les valeurs propres radiales représentent des ordres purement imaginaires et dont la 
partie imaginaire est non-entiére des fonctions de Bessel. Les normes se calculent hélas 
numériquement n fait que je n'ai pas trouvé de formules analytiques avec un poids 1/r : 


| RC rJ [ar Rau) 


Les fonctions HEH sont СА la forme : 
O ow (Ө) = A. А9 „Сов, ^ „0)+ Аб ,Sinh( x ,.0) 
as T,(r.0,2)- BéT(r,0,2). = 

= Ө (9) = a A, nr a COSh( Xu. 0) + В’ Sinh(À, 


ai T,(r,0,z) B; TG) a, zh 








nr,nz 0) 


= E, (0) 2 af A 


Le développement en série prend donc la forme : 


n, =+00 nr=+0 


Ty (r.0,2)= У Y (Courae Og rne (0) + Cus Orne (0))X В„ „ (F) X Z„ (2) 


n,=0 nr=0 


nr ,nz nr,nz 


Cosh(4' (Ө, -0))+ В’ Sinh(A; nr „ (9, -0)) 


0,nr „nz 
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Et les termes intégraux du développement en série prennent la forme : 





ÒO, nr, „2 (9) = ai A, nr „Сози e nz (6, -8))« В, ? пнд" Аһ „ (0, –0)) La Ё ci 1. 
[ar m | dz fé (r,2)Z,. (2) 
a! M E Sinh(4^,. nz (8, үт ө))+ 1 
Оо (9) AT 9 > Сон =>" 
+ BP Xa „a Cosh(X, „(бу -0)) Ape Rs Ze ON 














oT, (r.0,z)— Bs ed --f (r.z) 














O „„(8) = аА, ,. Cosh (o, 0)" PB Sinh X, 0) At 0 T а 200) 
Ө, r, „ (90) = ala ms Sinh(47,,.0)+ B; At OR 0) = С, nr,nz ( y [7 ( yf 
a nz m. nz r nz \ Z 
a; To(r.0,2)+ BI T(r,0,2)),_, = 702) 


avec d nz eh -" „Соһ(А u. „29, Ya? B? +U, ? BS )+ Sinh(X, m Ж, a, SE nr m, + Bí B; ) 
Ce qui permet de donner la solution du problème inhomogène angulaire dans le cas général : 
Z (Z) = aA, Coss z)* Bi Sinz)  K,(x)=K,(x)=Re(K, (x) et 1,(х)= Кеб, (x) 


X. valeurs propres solution de іп(А; 1 Yu асо, — Bo D; |- А Соѕ(А; 1 X í aj + SCH 


nz z nz nz z 


| caros Lo pismo + (ein) picos) ) 
2 
Z| 75 2 1 (a 





22 Соз(А А.) + BéSin(X.l Mp 1.)- BCos(4;.))- a; B; 


nz 


Ka D} 





K,()=K, (x) et 1(х)= Кеб„(х)) | 





zf -Ja £ 








R, (P) = —> гас -———= шс 
МЕ? ЖОО ОА) мал CLI BE; bo 
A "3 Qua) * Bab. Qul.) aA. АИА. NES 
| AE A) BT, TL) WADE. UM) F BaB a (АШ) 

e, (0) = af, „Cos, „(0, -0)+ Вот, „(9, ail 


(0) =a? A,  Cosbtä? 


nr,nz 


0)+ BL Sinh(A, 0 
A. = (Ж, „Соз, „0, Ya? B? +a? В). Sinha, „O ato? (A... вв") 


9, nr,nz nr,nz 


nr,nz 


eg r Raul) j de fr, DZ, C) e Ru.) | d: f. 22, C) 


O,nr,nz 








l,nr,nz 


А (=) Arr, nz 








А, „(Р 














2 nz 


R, | IZ, GY 


Toin (r, Ө, 2) = D У (C1270177) (9) zd Crus pi (0))х Ran (r) x Zu (z) 


n,-0 nr=0 


Il reste à considérer les cas particuliers où le développement en série nécessite la prise en compte des valeurs propres 
se 
nulles. Commençons par le cas de la valeur propre nulle ` z , comme d'habitude, elle ne peut s'envisager que 


2 2 2 
E А А = В =0 а= а zl 
lorsque les conditions homogènes axiales sont de Neumann de part et d'autre , B 0 B l 0 l 


implique nécessairement une fonction propre axiale constante , soit 1 (voir précédemment). 


, et cela 
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Dans ce cas les fonctions radiales ne sont également plus valides pour représenter les solutions du 
problème aux limites, car soit elles s'annulent soit elles présentent une singularité : 


Lim. T, (Ær)=0 


z r 
Anz 20 Za nz 


Lim. K, (Ui 


: r 
Anz 20 А Hz 


. Les fonctions propres radiales sont alors solutions de l'équation : 


r 2 £N 1 / А ` А А 
d’R(r) " 1 dR(r) A (a) R(r) - 0" qui ont déjà été exposées dans le cas à deux dimensions sur un 
dr? r dr r? 


secteur annulaire d'un disque creux : 


a = Ld =) or ZI 
La La 


Ro (r) = „Ао Cos ul 7) | T Pai Saal E) 


rl rl 








Cos(A; a ch BalaSin(a, до г) 


nr ,0 


Ко (т) = SEN 


nr ,0 


ТИ tq Sin(A^, lb o. UG, y SR E jÉ А, Cosi, up Yo, Bl, + Gu ba ) 


Les normes se calculent ainsi, en se rappelant que le poids utilisés dans le probléme de Sturm- 
Liouville est ici w(r)=1/r, de méme que pour le poids м(г)= 1/r pour le cas général qui conduit à des 
fonctions propres de Bessel d'ordre purement imaginaire : 





Ro 





l =a | dr (45, Cos(4^, „a c )+ B, Sin X, ост) 
0 


R 


nr,0 





[- gl [RB Jn Sinf ve 4 2 B; (ЖОООП 247 B: nä a) 


nr nr r пк nr 
2 Anr ‚0 


FA r f 
Comme Aw = „Ао et В v Bala 


| R = = (a, DG, y 3 B n эда) (о, (д y = Bau os(4^, да )+ 24,8, ola Sina, „a ) 
nr,0 


Les termes du développement de la série s'écrivent : 
Ou (0) =a? 45, Cosh(4;, (6, -0))+ B? Sinh(4, (6, -0)) 


nr 


Oro (8) = Ag Arr n Cosh(A;,0) + By Sinh(A;,0) 


nr,nz 


А =, Соза, 6, at B? +a? B2)+ Sinha, 0, ata? (x, + B2 B?) 








nr,0 








nr=+œ 


























Série de la forme D (С, E (0) * C, „Ө (Ө))х Ro (r) 
nr=0 
OR Bis ` ө "OR of) ` 9 
far = | dz fè (r. z) far == | dz ff (r. 2) 
Cond ; С EE 4 0 А А 0 
oordonnées r = Co „o 2 — j 
Anot. Raro O| Arola Rar O| 
1% 1. 1 1. 
a [dc К, (c) | dz fë (0.2) a [dc R, о) | dz f (r, z) 
Coordonnées т > Cono = E ° 5 ет = А 5 
А, ul. Rao O| A; JL. [Rao O| 
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Considérons maintenant la dégénérescence de ce méme cas avec une valeur propre radiale nulle, 
soit Ar = p alors la solution est de la forme A+B Log(r), et la fonction angulaire de la forme A+B 9, 


coincidant avec des conditions aux limites homogénes radiales de Neumann de part et d'autres, et 
une fonction radiale triviale : R(r)=1, dont on doit calculer la norme également avec le poids 


w(r)=1/r, car cette fonction propre est orthogonale avec les fonctions propres radiales . WË (09, 
comme il suit : 


rl 


Si ол=о0„=1 DB,-2p,-20 R(r)=1 R,(r) ЕВ] R, (т) = Со5(А,а т) 


À, tq Sin(1,a)=0= (Ry), R, (H = far O в, (r)2« far R (t) = a [dr Со5(А,а T) =0 
La r 0 0 


Les termes de valeurs propres nulles radial et axial, sont de la forme : 


A E= © r 
Termes de valeurs propres А, = А, =0 


Lala; В? + "i + пао) + B7 (0, -9))+C.. (ai + 88) 


L, 1 L L 
"dr 2 "dr 2 

Сосо = JÈ J dz 0.2 C, oo = [= Í dz f (r. z) 
ba 0 La 0 
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Ce qui permet de compléter la solution générale pour le cas particulier des conditions homogènes 
de Neumann axiales de part et d'autres : 


Si a =a =1 Ву = Вг = 























А 4 Sin. L)=0—= À = < E Z. (z) = dE | a ч = = 
D T К Hz X 
ROEK G) e Led, Jue jar es 
I, (Ar) K. Gi 
REE (r) = ze ` = Е gue ` T 
a „4,1 А, Mc = fal; EN 11) а „АК. А, MUS lr) - Bak. NC 14) 
д" f a, A, RC L5) + ps 7 Déi D (Al, 2) _ a, 2 nz 3 NE (Ql dei T BK yq (4.15) 
nr ,nz q z F ' = z D - ' 
à And gr Lil r ËD B, d „(= pe а as A. MC T4) D, Ky (ы r |) 





«= tod = | or d) Ry o (P) = 0A, ó Cm PE | + BL. sx Al 
rl rl rl rl 


Ratio, Cos(4^, o ch Bal in. „a г) 


nr,0 


Aso t4 Hd lk n „/ E ptus A „Cosa, a Ya 84. + a, Bal.) 














Sin\ 
Ко "= 2 af. G Аһ al +822) Set, (д, nr St -Bn L * leos( а) 20, B, m А апд, а) 
nr,0 
Ө, nr ,nz (9) = ai X, ,Cosh( m" ,nz (9, -0))+ B? $їлї(А, Ar, nz (Ө, -9)) 


Ө, ne (0) = 0 Xp me COSA y „„0) + Bo Sinh(A;, „Ө) 
A. =й, „Соз, 0 e B? +a? BE )+ Simi, „б, оде? (x... + BB") 
On „n. (0) = af X, Cosh(4^, (6, -0))+ B°Sinn(x;, (6, — Ө) 
Ө, „o (0) = AG Am Cosh(A;,0) + Bo Sinh(A7,0) 
Аг, i Cosh(4:,0, a9 B? +a? В )+ Sinha, 8, ata (x f + B? Be ) 
J Cas „| dr" zt D} — dz ff (r, LE: 


la 


Com = јао} "a ја уре, end? 


la 








x 


Coordonnées r = Com g = di = Ea fi (r.z) Ces Jar d EOD Ar noz) 
Li 


1 rl 
Coordonnées т = Co „ ъ= & | dt К„(т) j dzff(r,z) Cons = | dt К„(т) | dz f? (r, z) 
bi 0 la 0 


Solution développée ci-dessous : 


1 "ee Rs 
T, (650.2) = У (CO (0) С, (ED — Je) 
z nr-0 A, 0 К, 0 e 











Р, а DE 


2 n, =+90 nr=+00 


йа > > (Com. O nz Da nr,nz (0) + С, ‚пг SO (Ө))х 


L n,=0 nr=0 














EN ,nz Rr, nz 
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Et de compléter la solution générale pour le cas particulier des conditions homogénes de Neumann axiales et radiales 
de part et d'autres : 


Si aj aj QA 9,5 1 Bo Bí Ba В, 0 


A tq ЗША) EE %4, = LA Z. (z) = d ср J | 

















Ï 
2„(=)|' = 5 





Us nz 1.) d К eg? 10481,2) 
Al) E. А) 





x - P 
K, (x) = K; (x) et 1, (х) EE Re(1,, (x)) À, r tq Gë = 
nz°rl 
(Чы) E, Wr) | 
„ б Li) Em QS La) 


"nr ‚п 








1, б 
ат RC = Í dr Raa) 











a = Log 2 ат = 109 £ R,,,(r)= Cos Ш са R„ (1) = Соѕ(п,лт) 
La La | a La | 
Ao tq SinA,.0)-0e А, CC |е, „| = 2 
a 














Оо, „ (9) = aj À nr 2. Соз", „ (9, -0))+ p Sinh(4 А n -8) 
(8) = 0042, „„Соѕћ(А?, 0) + By Sinh(A,, 


4. („Сот „0, nt BE +a? Be) Sinh, „Ө, (обе? (a. + 828!) 
Oon (0) = 34^, Ee (0, -0))+ B? Sinn(x;, (6, Ө) 

6, m (0) = «А, ,,Cosh(4;,0) + By Sinh(A;,0) 

Ano = D Cosh( 2,0 al B? +a? Be )+ Sinh (A0, fata? (Y + pepe) 


PETS) T" п.л Ee ы 
Ciny É [ar [a fa (r. z)Cos| =z | Ci = lue 7 (r,z)Cos 
Ch fa F 0 L SE li La 0 


9, nr,nz nr ,nz Ө) 








N.T 
2 
x 
È 1 
r2 R 2 
Coordonnées r = Com g = Jar d Re fé fé(r,z) С, ло = [а 880 Las gea 
0 


l4 la 


1.5 l; 12 L 
Coordonnées т = Cono =d E Cos(n,z | dz ff (v,z) Со =X far Cos(n,z cl dz f? (т,2) 
e 0 


Gf aste. z) Go" f [oe (7,2) 


Avec la solution complète dans ce cas : 

_ Ï (Gros (o? 2s В’ (6, e ө))+ Co al ne pie) 
T3,(r,0,z) = 0 og D + 
al, (af pe + аб + BE BH.) 


Co 22 zeef) 
nr-oo ) a Is 


+ ЖЕЗ Z (Cour (0) "P C rei т» (0) x À" 


z = nr,0 











+ 


2 n, =+% nr=+œ0 Rp ,nz (r) x ES Д J 
c > 27 (Com. „O nz O, nr ,nz (9) + С, ‚пг mc (0))x 


l, n,=0 nr=0 














Aw ,nz К, nz 
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Solution du probléme de Dirichlet inhomogéne angulaire pour le secteur d'angle cylindrique 
creux et infini 


Il s'agit d'étudier la solution de l'équation de Laplace sur un corps illimité radialement comme sur 
la figure ci--dessous 











0 


Le problème inhomogène angulaire est le suivant : 
e^T(r,0,z) 1 OT(r,0,z) 1 e^T(r,0,z) eT(r,0,z) 
+ + + 
or? r Or r? 00° Oz? 
Q=fr/l, <r <+0}x{0/0<0<0,<27r}x{z/0<z<1,} 
CL. TU 0). =0 T(r,0,2},., = 0 T(r,0,z}, a = ff (r. z) 


T(r,0,z) = 0 T(r,0,z) " =0 





=0 T(r,0, z) fini 


z=0 z= 
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Les fonctions propres axiales dans ce probléme, simplifié sont les suivantes : 


Z. (ze Sin(4;.z) A. valeurs propres À. == |Z. c == 





i 
Les fonctions angulaires sont de la forme ( voir la séparation des variables ) : 
T(r,9,z), =0 = ©,„„(@Ө)= Sinh(r„ „.0) 





2 


Le développement en série prend donc la forme : 


n, =+90 nr=+00 


T(r, 0, z) = > Kä Cra Ore (8)x Кл: (r)x Zy (z) 


n,=0 т=0 


2 


Comme nous l'avons vu les fonctions propres radiales, sont des Fonctions de Bessel modifiées 


d'ordre purement imaginaire ^" (x)= К„(х)= ке(К„(х)) et 1,69) Re, (x) 


T.M.DUNSTER dans un article de 1990, « BESSEL FUNCTIONS OF PURELY IMAGINARY ORDER, WITH 
AN APPLICATION TO SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS HAVING A LARGE 
PARAMETER » introduit une autre fonction : 


vvécsec DG) our.) 


2Sin(zv) 
—VveR,xeR L,(x)- 35355 >С Lal sq (+ s 09) 
> Lx) RE x) LL, (x) TE Reli, (x) 


- 2Sin(izv) DRAN. 
>T (x)= веб, (x)) = 2609 7, (ху) 
Am 


Au valeur x->c», la fonction K tend vers 0, mais la fonction L et donc I diverge en exponentielle. 


(ze (est). tete] (ed) 


Cela conduit à ne pouvoir utiliser que la fonction K (de MacDonald) pour la construction d'une 
solution finie à grande distance radiale. La question est alors de savoir si la condition limite 
homogéne radiale conduit à un spectre discret de valeur propre. La réponse est oui car la fonction 
de MacDonald a bien un caractére oscillatoire selon le paramétre t, notamment lorsque ce 
paramétre est grand. Je n'ai pas trouvé d'argument théorique pour le prouvé, mais 
numériquement c'est bien ce que l'on constate. 


En revenant aux notations du probléme les fonctions propres radiales sont une combinaison 


linéaire de la forme : Ren NE Ks, (a. r) . Les fonctions radiales doivent respecter les conditions 


aux limites homogènes radiales sur la section creuse : 
2 2 r 2 
C.L. Kò. (Az. SR =0>7, tq K = l )= 0 Comme Tt dépend de A on note v,, Эт 


nz r nr ,nz 
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On voit que les valeurs propres radiales représentent des ordres purement imaginaires des 
fonctions de Bessel K et dont la partie imaginaire est non-entière. 


Relation d'orthogonalité des fonctions propres 


La séparation des variables dans les trois coordonnées cylindriques donnent trois équations 
différentielles ordinaires : 
O?T(r,0,z) 1 OT(r,0,z) 1 O?T(r,0,z) O?T(r,0,z) 
+ + + 
or? r Or r? 00? дг? 
=> T(r,0,z) = R(r)G(0)Z(z) 


d^n(r) А 1 dR(r) JS Hä Fo- 0 


=0 








dr? r dr r 








аӨ(Ө) ue 
20 - (v, f Ө(Ө)=0 
n.a 





Orthogonalité des fonctions propres ÍR, (rA, =X) ÍR (р), А, =2;) méme A. 


һу“, ) 0, [a] e "ei -o 





dr? r dr r 


©, 1 s. 8) Акко ] 











r dr r 
=> D vn (de) - R(r) ZO) 4-2 суру em (0, 


ағ dr 


=> LE оа (RO) Ï EE [И j реа, GËT 


ñ 1 


En d'autres termes il faut là encore utiliser un poids 1/r pour démontrer la relation d'orthogonalité. 
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Par la méme occasion, on redémontre une formule donnée par Prudnikov : 
SiR(r)-K,(r) R(r)-K,(r) ден Av 


f K, (r)K, (r) 1 dK, (0) dk. (r) | 
те : = => rr (009-9460 | 
Par la subsitution u = iu У — iv 


reu )K(r) 1 | 6 HEO к) ce) 


: r (и? =v?) 





1 









































fyc. Б =+œ=> Jar d e. pie = z x eO = | 
ja EE, OO m (e шу) Aën |) 
comme É Geiist 40) ик 0), (9) 





Lorsque u->v, alors il vient en développant autour de v en série de Taylor au premier ordre : 


fat пы. AO EEE dell 














u=v+dv > K (х)= KG. à 
_ jo ку " z| K, оја) ps a e) f K- (K, 2 | 








En développant се résultat, il vient : 





























ja Ele allt _ к к уде] ока) „о (К, =a 

x n ax Y 2v дА [ia дА |, ax 

Et de même lorsque les ordres sont imaginaires : 

je (К, (аг) = je (K, (»)) = ax K (a Ki (ax) ӘК, (ax) K (ax)- LK. (ax) 
с ы: у Qit | ^ OX Ua zeta où jog: E ax 
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Calculons une autre formule donnée également par Prudnikov, mais qui ne s'applique pas à notre 














probléme : 
+оо 2 
[ar "(К (ar = | 0) | = i t E |x, (ax)!  Prudnikov.vol.2.Formule.1.12.3.2 
x у у=ах 
Kx) акду) у 
De plus =K, (x)-K,„(x) => K (ax)=K,„(ax) 
dx dy |, ax 
Nn 2 x? у d 1 S 2 
> f arr (K lar} - ( LZ K (as)=K, (an) - [5*5 ка) 
2 \ ax 2 





Revenons aux normes des fonctions propres, elles se calculent maintenant comme suit grâce à la 
formule précédente, pour peu que l'on dispose de l'évaluation numérique de la dérivée 


paramétrique . 











R, (| = je 
L, 


























far а 
ш 4 l 217, nz 2 
с) 
nr „nz L 
2 E 
L 
Or K, ат -05 [al (ar) cet Dn к, 
^ l 5 2i Tes l ADE 


dnm: 


nr NZ 


-K, (5 
nr nz l 
a (22 











À TE Tar AZ 








AT pr nz 
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La solution du probléme aux limites prend la forme suivante : 


Z.()-Sin(z) Аза Sin(1)e02 x, = 


nz z nz 1 
2 


I 
Z (z)|^ = zi 





ICos(4;.1.] =1= 














nr ,nz 


Ò у: (9) = Sinh(r„. 9) Ru (r) = Ke. e r! (ЖУ tq K, e i| = 0 


2 










































































K, "л, = т 
“У.Е ` r 2i Tod L дА 
А=іт,, n 
Kei 
+00 İT nr nz l á 1, 
Cris [ar : faz HAE 2 
SE 1, É 0 L 
sie CSD Jr iz J 
T(r,0,z)= = > >. Z — 
L п,=1 el nm 
Sable, Al, ` E r) 
Pour un problème aux limites de Dirichlet de la forme : 
2 2 2 

Ô TOUS). 1 SZ Е Ô я Ô Të = (8,5) fini 

Or r Or r 00 Oz 
Q=fr/l << +о)х {0/0<0<0, «2z]x(z/0szxL] 
CL. T(r.0,z),., =0 T(r,0,z), , = ed LE T(r.0,z), = fo (r. z) 
T(r,9,z).., =0 T(r,9,z).., =0 
La solution s'écrit immédiatement : 
А m ci Tr nz tq K, са = 0 

1, | uz SÉ 
2 
2 ER Ka ar | =" | 

K, "=, Е [а І =_ "7 _ Ke, 77 L, á 

л L 1, r 2i Tu ,nz L Ze l, 0A 

А=іт nr ,nz 

















2 








N.T N.T 
+œ К. РМ | A r), +оо К... Ж | 1 r) L 
| z f D N.T | 2 ï 8 
Cu [ar S үй irch Í J Cimi far p |a felezi j 








п, =+00 n, =+œ d . Ka, * 
T(r 0 z) = 2 ` ` Cò. Sinh(r 1.0) Де Со» ne SERRE y ne (6, WE 0) 4 
5 Sinh(r,, 9 Sinh(r,,. „9, ) 


z nl п, =1 
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Supposons que les deux fonctions limites sont identiques, la solution est ainsi : 





x Е r) 
+оо ЇТ nr nz L 
f? (r, z)- fé (r, z)= ris zl С: — [ar CA dz TUE J 
r 0 z 


ir... C 


nr „nz ` nr ,nz 


ôK e l | 
À pn 
l 
4 n, =+0 n, =+% | И, к, AJ 1 Z x Sinh(r,...0,) 


лі. 


r hn; n 











A=i T nr nz 


H e dul 22 eat sit. (Ө, —Ө)) 











2 2 





Si les deux fonctions limites sont constantes, il vient : 


L 
fé(r.z)- ff(r.z)- f(r,z)=1= faz Sal? I ( -(=1)"- je 0 pour n, impair 
0 


nz 














A=i T nr „nz 














K. É r) 
C, nr,2nz+1 — C, nrj2nz4 — Cur nz — [ar = A. T SE T nz Iq K, [Enk 8 E 0 
vx uh i ; r L SI L 
| » (® m e) S (es +), | 
+00 IT 1 2 À r 
K. (2n, +1 far S Мы, ко (2n, +1 ] L 
i l, 1 r 2060 L д5 
Їн СЕ x 
à x [en lr i) 
gr "сел Cn +D Kus É SC a : x Sinh(r,, „0,) 
T(r,0,z)- = Ges L дА 
LT n,-0 n,=1 
А=іт,, nz 
| . ( Qn, +1)r (2n, +1)л 
x (Sinh(r,,....0)-- Sinh(r,,.,. (0, - 0)))$їл `= =Z Ks EI 
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Si les deux fonctions limites sont identiques et présentent un profil en palier, alors la solution 
s'écrit : 
































N.T 
: Z ïy 
iede 1 sus z}e [n.n |х [2,,2, | eg far LUI 
0 si {г} ë [n.n xl, z;] f z r 
ЇТ y Cu! COS SC z, |- Cos Se Zs 
2 2 х 
ӘК, е 1, 
A п. =+00 п, =+00 2 пл А | 
CERERI K. z | Sinh 9 
T(r.0, z) = 2 І, Z уу n; ir, as tl l, r SCH хып bk ,) 
Ait nr nz 
Xx Kien = r |Sin - 2 (Sinh(r „ „Ө )+ Sinh(r „. „. (6, -0))) 





A titre d'exemple voici la liste des valeurs propres pour un système de 30x60 termes, pour une 
valeur de hauteur l,=3, un rayon intérieur l,=0.1 et une ouverture du secteur d'angle de r2 : 


АПИЗ mu 
[(1.15851,2.06881,2.88156,3.6409,4.3647,5.0624,5.73971,6.40044,7.0473,7.68228,8.30693,8.92245,9.52982,10.1299,10.72 
32,11.3105,11.8921,12.4686,13.0402,13.6073,14.1702,14.7292,15.2844,15.836,16.3844,16.9295,17.4717,18.0109,18.5475, 
19.0813,19.6127,20.1416,20.6683,21.1927,21.7149,22.235,22.7532,23.2694,23.7838,24.2963,24.8071,25.3162,25.8236,26. 
3295,26.8337,27.3365,27.8378,28.3376,28.8361,29.3332,29.8289,30.3234,30.8166,31.3086,31.7994,32.289,32.7774,33.264 
7,33.7509,34.236] 

Au (2 n)/3 -»t,,,- 
(1.46832,2.52997,3.46526,4.33307,5.15643,5.94733,6.71301,7.45825,8.18644,8.90006,9.60103,10.2909,10.9708,11.6417,1 
2.3046,12.96,13.6086,14.251,14.8875,15.5186,16.1446,16.7658,17.3825,17.9949,18.6033,19.2079,19.8089,20.4064,21.000 
7,21.5917,22.1798,22.7649,23.3473,23.927,24.5041,25.0787,25.651,26.2209,26.7886,27.3542,27.9176,28.479,29.0385,29. 
596,30.1517,30.7056,31.2577,31.8082,32.3569,32.9041,33.4496,33.9937,34.5362,35.0773,35.6169,36.1551,36.692,37.2275 
:,31.7617,38.2947) 

п ->т 
{1.71919,2.8875,3.90875,4.8523,5.74493,6.6005,7.42733,8.23093,9.01517,9.78291,10.5363,11.2772,12.0068,12.7263,13.4 
367,14.1387,14.8331,15.5204,16.2011,16.8757,17.5446,18.2082,18.8666,19.5204,20.1696,20.8145,21.4554,22.0924,22.725 
7,23.3555,23.9819,24.6051,25.2252,25.8422,26.4564,27.0679,27.6766,28.2828,28.8865,29.4878,30.0868,30.6835,31.278,3 
1.8703,32.4606,33.049,33.6353,34.2198,34.8024,35.3832,35.9622,36.5396,37.1152,37.6893,38.2617,38.8326,39.402,39.96 
99,40.5363,41.1013) 
A, (4 m)/3 ->т „= 
(1.94145,3.19513,4.28504,5.28894,6.23665,7.14354,8.01884,8.86863,9.69719,10.5077,11.3025,12.0835,12.8523,13.61,14. 
3578,15.0964,15.8267,16.5493,17.2648,17.9735,18.6761,19.3728,20.0641,20.7501,21.4313,22.1078,22.7798,23.4477,24.11 
16,24.7717,25.428,26.0809,26.7304, 27.3767, 28.0198,28.66,29.2972,29.9317,30.5634,31.1926,31.8192,32.4434,33.0652,33 
.6847,34.302,34.9171,35.5301,36.141,36.7499,37.3569,37.962,38.5653,39.1667,39.7664, 40.3643, 40.9606, 41.5553, 42.1483 
‚42.7398, 43.3297} 
A (5 п) /3 -5t,., 
(2.14618,3.47231,4.62042,5.6754,6.66966,7.61988,8.53605,9.42475,10.2906,11.137,11.9666,12.7814,13.5831,14.3729,15. 
152,15.9214,16.6818,17.4339,18.1784,18.9158,19.6466,20.371,21.0896,21.8027,22.5106,23.2134,23.9116,24.6053,25.2947 
,25.98,26.6614,27.339,28.0131,28.6837,29.3509,30.015,30.676,31.334,31.9892,32.6415,33.2912,33.9383,34.5828,35.2249 
,35.8646,36.502,37.1372,37.77701,38.401,39.0297,39.6565,40.2813,40.9041,41.5251,42.1443,42.7617,43.3773,43.9913,44. 
6035,45.2142} 

Ag 2 Hog Vos 
(2.33882,3.72857,4.92775,6.0275,7.06251,8.05063,9.00252,9.92521,10.8236,11.7014,12.5614,13.4056,14.2359,15.0536,15 
.86,16.6561,17.4427,18.2206,18.9904,19.7526,20.5078,21.2564,21.9988,22.7353,23.4663,24.1921,24.9128,25.6289,26.340 
4,27.0476,27.7507,28.4498,29.1452,29.8369,30.525,31.2099,31.8914,32.5698,33.2452,33.9177,34.5873,35.2542,35.9184,3 
6.58,37.2392,37.8958,38.5502,39.2022,39.852,40.4996,41.1451,41.7885,42.4299,43.0693,43.7068,44.3424,44.9762,45.608 
2,46.2384,46.8669} 











X= (7 m/3 ->t 


{(2.52253,3.96938,5.2144,6.35429,7.42582,8.44787,9.43 


.505,17.3249,18.1348,18.9356,19.7279,20.5123,21.2892, 


73,10.3848,11.3124,12.2182,13.1052,13.9757,14.8315,15.6742,16 
22.0593,22.8229,23.5803,24.3319,25.078,25.8189,26.5548,27.286 


,28.0127,28.7351,29.4533,30.1675,30.878,31.5847,32.288,32.9878,33.6843,34.3777,35.068,35.7553,36.4398,37.1215,37.8 
004,38.4768,39.1506,39.8219,40.4909,41.1575,41.8218,42.4839,43.1438,43.8016,44.4574,45.1111,45.7629,46.4127,47.060 
7,47.7069, 48.3513} 
A (8 п) /3 ->т, „= 
(2.6993,4.19824,5.48504,6.66151,7.7663,8.81924,9.83218,10.8129,11.7669,12.6982,13.6098,14.5042,15.3832,16.2485,17. 
1013,17.9428,18.774,19.5955,20.4082,21.2126,22.0094,22.7989,23.5816,24.3579,25.1282,25.8927,26.6519,27.4058,28.154 
8,28.8991,29.6389,30.3744,31.1057,31.8331,32.5567,33.2765,33.9929,34.7058,35.4154,36.1218,36.8252,37.5255,38.223,3 
8.9176,39.6096,40.2988,40.9855,41.6697, 42.3515, 43.0309, 43.708, 44.3828, 45.0555, 45.726,46.3944,47.0607,47.7251,48.38 
76,49.0481,49.7068) 

z 3 п Tar, 

{2.87051,4.4175,5.74286,6.95307,8.08851,9.16991,10.2096,11.2158,12.1942,13.1489,14.0831,14.9994,15.8998,16.7858,17 
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.659,18.5203,19.3709,20.2115,21.043,21.8658,22.6806,23.488,24.2883,25.082,25.8694,26.6508,27.4266,28.1971,28.9624, 
29.7228,30.4786,31.2299,31.977,32.7199,33.4588,34.1939,34.9254,35.6533,36.3778,37.099,37.817,38.5319,39.2438,39.95 
28,40.659,41.3624,42.0632,42.7613,43.457,44.1502,44.841,45.5295,46.2157,46.8997,47.5815,48.2612,48.9389,49.6145,50 
-2882,50.96} 

A (10 m/3 ->т,, „= 
(3.03713,4.62886,5.99012,7.23173,8.39569,9.50356,10.5682,11.598,12.599,13.5754,14.5307,15.4674,16.3875,17.2929,18. 
1848,19.0646,19.9332,20.7915,21.6403,22.4802,23.3119,24.1358,24.9524,25.7621,26.5654,27.3625,28.1538,28.9395,29.71 
99,30.4953,31.2659,32.0319,32.7934,33.5507,34.3039,35.0531,35.7986,36.5404,37.2787,38.0135,38.7451,39.4734,40.1987 
,40.9209,41.6403,42.3568,43.0706,43.7817,44.4902,45.1961,45.8996,46.6007,47.2995,47.9959,48.6902,49.3822,50.0721,5 
0.76,51.4458,52.1296) 

A (11 п)/3 ern 
(3.19989,4.83359,6.22851,7.49955,8.69025,9.82292,10.9109,11.9628,12.985,13.9818,14.9567,15.9125,16.8511,17.7745,18 
.6841,19.5811,20.4666,21.3414,22.2064,23.0623,23.9096,24.749,25.5808,26.4056,27.2236,28.0353,28.8411,29.6411,30.43 
56,31.225,32.0094,32.789,33.5641,34.3348,35.1012,35.8637,36.6222,37.3769,38.128,38.8756,39.6198,40.3607,41.0985,41 
.8331,42.5647,43.2934,44.0193,44.7425,45.4629,46.1808,46.8961,47.609,48.3194,49.0275,49.7333,50.4368,51.1382,51.83 
74,52.5345,53.2296) 





z nr,nz 


(3.35934,5.03263,6.45931,7.75813,8.97403,10.1301,11.24,12.3129,13.355,14.371,15.3644,16.3381,17.2942,18.2345,19.16 
07,20.0739,20.9752,21.8656,22.7459,23.6168,24.4788,25.3327,26.1788,27.0177,27.8497,28.6752,29.4945,30.3079,31.1157 
;31.9182,32.7155,33.508,34.2958,35.0791,35.858,36.6328,37.4036,38.1705,38.9337,39.6932,40.4493,41.202,41.9514,42.6 
976,43.4408,44.181,44.9182,45.6527,46.3844,47.1134,47.8398,48.5637,49.2851,50.0041,50.7207,51.435,52.1471,52.857,5 
3.5648,54.2704} 

AQ 13 5)43 cT, = 
(3.51594,5.22675,6.68355,8.00868,9.24847,10.4267,11.5574,12.65,13.711,14.7452,15.7562,16.7468,17.7194,18.6758,19.6 
177,20.5462,21.4625,22.3676,23.2623,24.1474,25.0234,25.891,26.7507,27.6029,28.448,29.2865,30.1186,30.9447,31.765,3 
2.5799,33.3895,34.1941,34.9939,35.7789,36.5798,37.3663,38.1486,38.927,39.7015,40.4724,41.2397,42.0035,42.77639,43.52 
11,44.2752,45.0262,45.77742,46.5193,47.2616,48.0011,48.738,49.4723,50.204,50.9333,51.6602,52.3847,53.1069,53.8268,5 
4.5446,55.2602} 

M= (14 m)/3 —1,,- 
(3.67002,5.41655,6.90201,8.2522,9.51473,10.714,11.8646,12.976,14.0549,15.1064,16.134,17.1408,18.1291,19.1008,20.05 
75,21.0005,21.9311,22.8502,23.7586,24.6571,25.5464,26.427,27.2996,28.1644,29.022,29.8728,30.7171,31.5552,32.3874,3 
3.214,34.0353,34.8514,35.6626,36.4691,37.271,38.0685,38.8618,39.6511,40.4365,41.218,41.9959,42.77703,43.5412,44.308 
8,45.0732,45.8344,46.5926,47.3478,48.1001,48.8497,49.5965,50.3406,51.0822,51.8212,52.5577,53.2919,54.0237,54.77532, 
55.4804,56.2054} 
REDIÉ a 
(3.82187,5.60252,7.11538,8.4895,9.77375,10.9931,12.1626,13.292,14.3881,15.456,16.4995,17.5217,18.525,19.5112,20.48 
21,21.4391,22.3832,23.3156,24.2371,25.1485,26.0504,26.9435,27.8282,28.7051,29.5746,30.4372,31.293,32.1426,32.9862, 
33.824,34.6563,35.4834,36.3054,37.1226,37.9352,38.77433,39.5471,40.3467,41.1423,41.9341,42.77221,43.5065,44.2874,45. 
0649,45.839,46.61,47.3779,48.1427,48.9046,49.6637,50.4199,51.1734,51.9243,52.6726,53.4184,54.1617,54.9026,55.6412, 
56.3774,57.1115) 

A= (16 m/3 ->т,, „= 
{3.97174,5.78509,7.32419,8.72124,10.0263,11.265,12.4526,13.5991,14.7117,15.7953,16.8541,17.891,18.9086,19.9088,20. 
8933,21.8635,22.8207,23.7658,24.6998,25.6234,26.5374,27.4424,28.3388,29.2273,30.1082,30.9819,31.8489,32.7094,33.56 
38,34.4123,35.2552,36.0927,36.9251,37.77526,38.5753,39.3935,40.2073,41.0168,41.8223,42.6238,43.4215,44.2155,45.0059 
145.7928,46.5764,47.3567,48.1338,48.9078,49.6788,50.4469,51.2122,51.9747,52.7345,53.4916,54.2462,54.9983,55.7479,5 
6.4951,57.24,57.9826] 
A, (17 m)/3 tu 
(4.11982,5.96458,7.5289,8.94799,10.273,11.5302,12.7352,13.8983,15.0266,16.1254,17.1988,18.2499,19.2812,20.2947,21. 
2923,22.2752,23.2449,24.2022,25.1482,26.0836,27.0092,27.9255,28.8332,29.7327,30.6245,31.509,32.3866,33.25777,34.122 
4,34.9812,35.8342,36.6818,37.5242,38.3615,39.194,40.0218,40.8452,41.6642,42.4791,43.2899,44.0968,44.9,45.6995,46.4 
955,47.2881,48.0773,48.8633,49.6461,50.4259,51.2027,51.9766,52.77476,53.5159,54.2816,55.0446,55.8051,56.563,57.3185 
,58.0717,58.8225)] 

SES (Be Ibo. ре 
(4.26628,6.14129,7.7299,9.17021,10.5145,11.7895,13.0113,14.1902,15.3337,16.4471,17.5346,18.5993,19.6438,20.6702,21 
.6803,22.6755,23.6572,24.6263,25.5838,26.5305,27.4672,28.3945,29.313,30.2231,31.1254,32.0203,32.9081,33.7892,34.66 
39,35.5326,36.3954,37.2526,38.1045,38.9513,39.7931,40.6303,41.4628,42.291,43.1149,43.9347,44.77505,45.5625,46.3708, 
47.1755,47.9766,48.7744,49.5689,50.3602,51.1484,51.9335,52.7157,53.495,54.2715,55.0453,55.8164,56.5849,57.3509,58. 
1144,58.8754,59.6341} 

M= (I9 YA E ac 
(4.41126,6.31546,7.92752,9.38832,10.7512,12.0434,13.2813,14.4756,15.6337,16.7612,17.8622,18.94,19.9973,21.0361,22. 
0583,23.0654,24.0586,25.039,26.0076,26.9653,27.9128,28.8506,29.77795,30.6999,31.6123,32.5171,33.4148,34.3056,35.189 
9,36.068,36.9402,37.8068,38.6679,39.5238,40.3746,41.2207,42.0621,42.899,43.7316,44.5601,45.3845,46.205,47.0217,47. 
8347,48.6442,49.4503,50.2529,51.0524,51.8486,52.6418,53.432,54.2192,55.0036,55.7852,56.5641,57.3404,58.1141,58.885 
2,59.6539,60.4202} 

X (20 0/3 ->т,, „= 
(4.5549,6.48732,8.12205,9.60265,10.9835,12.2924,13.5459,14.755,15.9273,17.0684,18.1825,19.273,20.3425,21.3934,22.4 
272,23.4457,24.4501,25.4415,26.4208,27.389,28.3468,29.2949,30.2338,31.1641,32.0863,33.0007,33.9079,34.8081,35.7017 
:36.589,37.4703,38.3458,39.2157,40.0804,40.94,41.7947,42.6446,43.49,44.331,45.1678,46.0005,46.8292,47.654,48.4751, 
49.2926,50.1066,50.9172,51.7245,52.5286,53.3295,54.1274,54.9223,55.7143,56.5034,57.2899,58.0736,58.8547,59.6332,60 
.4093,61.1828) 

A es IB T p 
(4.69729,6.65705,8.31373,9.81352,11.2117,12.5367,13.8054,15.0289,16.2148,17.3691,18.4959,19.5987,20.6802,21.7426,2 
2.7878,23.8173,24.8325,25.8344,26.8242,27.8026,28.7704,29.7284,30.677,31.6168,32.5484,33.4722,34.3885,35.2978,36.2 
004,37.0966,37.9866,38.8708,39.7494,40.6225,41.4905,42.3535,43.2118,44.0654,44.9145,45.7593,46.6,47.4366,48.2693,4 
9.0982,49.9235,50.7452,51.5634,52.3783,53.1899,53.9983,54.8036,55.6059,56.4052,57.2017,57.9954,58.7863,59.5746,60. 
3603,61.1434,61.9241} 

X (22 п)/3 ->т,, „= 
(4.83854,6.8248,8.50279,10.0212,11.4363,12.7769,14.0603,15.2976,16.4969,17.6639,18.803,19.9177,21.0107,22.0844,23. 
1406,24.1808,25.2064,26.2187,27.2185,28.2068,29.1844,30.1518,31.1099,32.059,32.9997,33.9325,34.8578,35. 7759, 36.687 
1,37.5919,38.4904,39.383,40.2699,41.1513,42.0274,42.8985,43.7647,44.6263,45.4833,46.3359,47.1843,48.0285,48.8689,4 
9.7054,50.5381,51.3673,52.1929,53.0151,53.834,54.6496,55.4622,56.2716,57.0781,57.8816,58.6823,59.4803,60.2755,61.0 
681,61.8581,62.6456} 

M= (23 m/3 Sec 
(4.97873,6.99074,8.68942,10.2259,11.6574,13.0132,14.3108,15.5617,16.7739,17.9532,19.1043,20.2306,21.3348,22.4194,2 
3.4862,24.5368,25.5726,26.5948,27.6044,28.6023,29.5893,30.5661,31.5332,32.4914,33.441,34.3825,35.3164,36.243,37.16 
27,38.0758,38.9825,39.8833,40.7782,41.6676,42.5517,43.4306,44.3046,45.1738,46.0384,46.8986,47.7545,48.6062,49.4539 
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,50.2977,51.1378,51.9741,52.8069,53.6363,54.4622,55.2849,56.1044,56.9208,57.77342,58.5446,59.3521,60.1568,60.9588,6 
1.7581,62.5548,63.3489] 

Age B 5те 
(5.11792,7.15498,8.87379,10.4278,11.8753,13.2459,14.5573,15.8214,17.0461,18.2375,19.4002,20.5377,21.6528,22. 748,23 
.8252,24.8859,25.9316,26.9635,27.9826,28.9898,29.986,30.9717,31.9478,32.9146,33.8729,34.8229,35.77652,36.77001,37.62 
8,38.5491,39.4639,40.3725,41.2753,42.1725,43.0642,43.9507,44.8323,45.7709,46.581,47.4485,48.3117,49.1706,50.0255,50 
.8764,51.7235,52.5669,53.4067,54.2429,55.0758,55.9053,56.77315,57.5547,58.3747,59.1918,60.0059,60.8172,61.6257,62.4 
315,63.2347, 64.0352} 

A," (25 п) /3 ->т,, „= 
(5.2562,7.31764,9.05604,10.6272,12.0902,13.4751,14.8001,16.0769,17.3139,18.517,19.691,20.8394,21.9652,23.0708,24.1 
58,25.2285,26.2839,27.3252,28.3535,29.3698,30.3748,31.3694,32.3541,33.3294,34.296,35.2543,36.2048,37.1477,38.0836, 
39.0126,39.9352,40.8515,41.7619,42.6666,43.5658,44.4598,45.3486,46.2326,47.1118,47.9865,48.8567,49.7227,50.5845,51 
.4424,52.2964,53.1465,53.9931,54.836,55.6756,56.5117,57.3446,58.1742,59.0008,59.8243,60.6449,61.4626,62.27775,63.08 
96,63.899,64.7058} 
A," (26 n)/3 T aeu 
(5.39362,7.47881,9.23633,10.8242,12.3023,13.7013,15.0394,16.3287,17.575,18.7921,19.9771,21.1362,22.2724,23.3879,2 
4.485,25.5651,26.6298,27.6803,28.7176,29.7427,30.77565,31.7596,32.77526,33.77363,34.7111,35.6774,36.6358,37.5866,38.5 
302,39.4669,40.3971,41.3209,42.2387,43.1508,44.0573,44.9584,45.8544,46.77454,47.6316,48.5133,49.3904,50.2632,51.131 
9,51.9964,52.8571,53.7139,54.567,55.4165,56.2625,57.1051,57.9444,58.7805,59.6134,60.4432,61.27,62.0939,62.915,63.7 
333,64.5488,65.3617] 
A; 9 m tL. = 
(5.53023,7.63861,9.41476,11.0189,12.5118,13.9244,15.2754,16.5769,17.8373,19.0631,20.2588,21.4283,22.5745,23.6999,2 
4.8065,25.896,26.9698,28.0293,29.0754,30.1091,31.1313,32.1427,33.144,34.1358,35.1185,36.0927,37.0589,38.0173,38.96 
85,39.9126,40.8502,41.7813,42.7064,43.6256,44.5391,45.4473,46.3502,47.2481,48.1412,49.0296,49.9135,50.7929,51.6682 
192.5393,53.4065,54.2698,55.1293,55.9852,56.8376,57.6864,58.532,59.3742,60.2133,61.0493,61.8822,62.7122,63.5392,64 
.3635,65.185,66.0038} 

X (28 n)/3 ->T„ n= 
{5.66607;,7.79709,9.59146,11.2115;12.7188,14.1448,15.5083,16.8217,18.0935;19.3301,20.5363;,21.716,22.8721,24.007,25. 
123,26.2216,27.3043,28.3725,29.4272,30.4693,31.4998,32.5194,33.5286,34.5283,35.5188,36.5006,37.4743,38.4403,39.398 
8,40.3503,41.295,42.2333,43.1654,44.0915,45.012,45.927,46.8367,47.7414,48.6411,49.5361,50.4265,51.3125,52.1942,53. 
0718,53.9453,54.8149,55.6807,56.5428,57.4013,58.2563,59.1079,59.9563,60.8014,61.6433,62.4822,63.318,64.151,64.9811 
,65.8084,66.633)] 

M= (29 m)/3 >т, à 
(5.80121,7.95435,9.76651,11.4021,12.9235,14.3626,15.7383,17.0633,18.3462,19.5934,20.8099,21.9995,23.1652,24.3096,2 
5.4346,26.5421,27.6336,28.7103,29.7734,30.8237,31.8623,32.8898,33.9069,34.9142,35.9123,36.9016,37.8827,38.8559,39. 
8217,40.7803,41.732,42.67773,43.6163,44.5493,45.4765,46.3982,47.3145,48.2257,49.132,50.0334,50.9302,51.8226,52.77106 
,93.5944,54.4741,55.3498,56.2217,57.0899,57.9545,58.8155,59.673,60.5273,61.3782,62.226,63.0707,63.9123,64.751,65.5 
868,66.4198,67.25} 

A 10 п ^n 
(5.93567,8.11045,9.94002,11.5909,13.126,14.5779,15.9656,17.3019,18.5957,19.8533,21.0799,22.2791,23.4543,24.6078,25 
.7417,26.858,27.958,29.043,30.1143,31.1727,32.2191,33.2543,34.2791,35.2939,36.2994,37.2961,38.2844,39.2648,40.2375 
,41.2031,42.1617,43.1138,44.0595,44.9992,45.933,46.8612,47.7784,48.7017,49.6143,50.522,51.4251,52.3236,53.2178,54.1 
077,54.9934,55.8752,56.77531,57.6272,58.4976,59.3645,60.2279,61.0879,61.9446,62.7981,63.6484,64.4957,65.34,66.1814, 
67.0199,67.8556} 





A titre d'exemple voici le profil reconstitué de la fonction limite par le développement en série 
bidimensionnelle, en prenant 30 termes еп z et 60 termes еп r, pour une valeur de hauteur 1,=3, un 
rayon intérieur 1,=0.1, une ouverture du secteur d'angle de л/2 et un palier fixé à r:=21,, 
r2-5l, 21=172/3, z:=21z/3 : 
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Voici une coupe en r, du profil de solution du méme probléme aux limites symétrique, à mi angle 8 
-n/4: 





0.1 0.2 0.3 0.4 0.5 0.6 D€ 0.8 
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Et une coupe en г,2 du profil de solution du méme problème aux limites symétrique, à mi-hauteur : 


N 


N 











0.1 0.2 0.3 0.4 0.5 о.е 0 


La vue en coupe de la section cylindrique en coordonnées cartésiennes à demi-hauteur se présente 
ainsi : 


0.7 





D.0 0.2 0.4 0.8 0.8 
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Quelques intégrales des fonctions de Bessel 


Les formules les plus couramment utilisées : pour calculer les intégrales des fonctions de Bessel on 
utilise les formules suivantes : 


@ Гага) = [о] (dz anoo] 
D'autre part faz zJ,'(z)= [zJ, (2)]- [а J (z) intégration par partie 
J,'(z) 2 J, ,(z)- А J (z) = faz zJ,'(z)- [а TE — E о) = [а zJ, az) -vÍ dz J (z) 
2 2 
= [27,02)]- fd J, (2) = [dz at, €) - v | dz Ate 
=> faz zJ, (2) =[zJ, (3)]+(v - Df dz J,(2) 
(2) (dz zJ Da] [az J, C) 
On obtient aussi avec J ,'(z) = Е J (2)- tx) 
2 


(3) [а 27,102) = BI GO (v + )f dz J (2) faz zJ,(z)= -[2/,(@]+у{ dz TRA A) 


(4) [а J (2) = 2) Jaaa) > [а J, (2) = D Jari (x) 


> [а Л, (z) = 25, Janar Œ) faz Jon (2) = 25. Јо (5) 
0 к=0 0 к=0 


(5) [а 2 J (z) = Za (2)|+ (v +1- a) de STI GH А>0 méme procédé intégration par partie 
(6) [а Sekt Jon rb (v -1+ DI dz z^ J,.(z) А>0 méme procédé intégration par partie 


(7) faz Ju ()- [а J, (z) = dit: S dus 0) = Р Jari (X) - 3m H = ay Jara (X) 
=> je J,,(z) = je 12-25 Лк (x) 


Í d IAEI at 
k=0 
(8) ` 
У и Уа) 
0 k=0 k=n 


* * n-l 
^ x n-l [& J, (2) E [а Јол (2) is 2» Jaara) 
k=0 
= | dz J,(z)- | dz Ja 0) 24, Q1 А 


b 


fa J, (z) = faz J aan (z) + ў Dis Cf 


+00 


(9) j dz Jani (2) | dz J (z) = À Jua) - Y о H m 29 з (x) = 29 Ох) 


k=0 


> (dz Jama (2) [d 40-2944 0071- 4,007244 
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L'application des formules (5) et (6) donne : 
(5) fd ID= |27, +y -A)| d zS, 


(9 fero- ,(@)]+(v-1+4)f dz z7J, Œ) 

(6) = [dz z^ J, (2) =- [50,0| 1+) )[ dz z^?J, (z) 

(52 [az z^?J (2) = |22, 2) +3-А)[ dz Ss 

[а 272,0 - - "o, cone (и-1+4) (2-27, (2)] (+3) )f az z^?J. 40) 
IER DEE и +3) = [40 12s 727,0) 
L'application des formules (2) et (4) donnent : 


(2) Ја 2091.) о [а,ә 46 € 1.0 avec = Аг 


(4) јело DANCE jé ane; [arro 2 > ыйы (Аз) 


ege innato] spi en nr aan) 


De méme : 


Posons I(u,v,x)= Је z" A EH 





ajaaa (4 2)" J,(Az)- Je z" J (2) 


Am 


dz z" Ј (А2) = =a 
j А 


I(u,v,Ax) 


dz z" Ј (А2) = т 
je taa е 


Problèmes aux limites de Laplace en coordonnées cartésiennes et polaires cylindriques 3D — p 212 


Des formules plus générales d'intégration des fonctions Bessel d'ordre quelconque sur des 
fonctions polynomiales sont les suivantes et peuvent être utile à connaître : 

















| (een . (eeu) 
Ја 2, (s) х MA +2k+1) SES Jl Ушу tq Re(v-u+1)>0 
o r үү ле = 
S 2 

: r ei - (ena) 
—u-2, |а PJ (dex >> M +2k+1) Jobi) 

0 (H ) nt (a +k+2) 

2 
Comme T| =Â YA. +1 и а. Г ү un = de x E Г шүү = 
2 2 2 2 2 2 


_ v+1 pud v+l cr. 
2. 2 2 


J k+1 
jé zum: 2x кыр о; 


В гэ) D do 
mul (dzzJ(23)=x=5F- +Y (у 21) aal) 
0 


rz) ч d +k+ 2) 

2 2 

Comme T| -«1|--T| - | T| к+2 |= Le kl кв |=[LZ+k+1 (Z+k4 LE 
2 2\2 2 2 2 2 2 2 








(v -- 2k +1) 
dz zJ,(z) 3 
| з олт с BEE 
e) 
=> u =0, EE (z)- Y (v +2Е+1) Jana (X) 
md (eet 


Comme (e +k+ 1) = e 4 dE + d 
2 2 2 


PERAS UE DIT] Sei nr -2WJ 
| z zJ,(z) d * 1520 viaa (X) H уо (X) 


L'inconvénient d'utiliser ce type de formules pour évaluation numérique peut résider dans la 
lenteur à obtenir la convergence de l'évaluation de l'intégrale, c'est la raison pour laquelle on 
cherche parfois à calculer ces intégrales à l'aide d'autres fonctions connues de la physique 
mathématique, comme par exemple les fonctions hypergéométriques, les fonctions de Struve oü 
les fonctions de Lommel. On peut alors bénéficier d'une implémentation efficace de ces fonctions 
dans un logiciel de calcul mathématique. 
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Utilisation des fonctions hypergéométriques 


On a également à l'aide des fonctions hypergéométriques, les formules suivantes : 














x А+у+1 2 
[dz zJ, = - 2 VE ЗАВИН ee si Bell 4v)» -1 
d 2" (А +у +1)Г(у +1) 2 2 4 
b А+у+1 2 \ P? 
1+4 À 
= [ 277,02) = i T E mp ы ДУ. si Re(2+v)> 
А 2' (A v - DI(v +1) 2 2 il 





b v+l 2 ? 
1 
SCHER id d SE | si Re(v)>-1 











2'T(v+2) 4 2 2? 4 
^ grs 24v iw 
= [427,(2)= E, pp e si Re(v41)» -1 
d 2" (v +2)Г(у +1) 2 2° 4)| 
fe) 
+œ А+у+1 2 
(224) (3) =? 2 € x F, ERP A aE si 
J =” 2" (Av +1)Г(у +1) 2 2 4 
2 


. (e 
= faz z^J (2) 22^ EE C. 
0 


=” 
2 


si -1- Ré) е) « == [de 2.) -1 
0 





-1 


Re(A)« 
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Utilisation des fonctions de Lommel (cas particuliers des fonctions de Struve) 


Dans l'ouvrage de A. Erdelyi. H. Bateman « HIGHER TRANSCENDENTAL FUNCTIONS VOL Il », les 
fonctions de Lommel de premiére et deuxiéme espéce sont décrites par les formules suivantes 
page 40 : 

Fonctions de Lommel | s, et S,, 





dz? 2 


dw 1а В 
s, (z) solution de l'équation inhomogéne de Bessel uan E + К т. 5 |» 2g 
2 Z az 
Toute fonction s, ,(z)+ АЈ, (2) + BY,(z) ои s,,(z)* AJ ,(z) + BJ,(z) est également solution 


_ > e _ i z* R k=% E i SEH 
Siv (2) =Z Kë 1) m=k = ( 1) m-k 
= — [[T(4+2m-1f-v2) = T(4+2m+1f =v?) 


т=1 


DR (+ Zei (еч) 

2 А 2 

= A-1 SE 

äs 2 UT Ze. 
2 




















x TE 


ge Ai ni KEEN d 
Ау Ta) = 2142 2 

L'expression présente un pôle lorsque : 
(4+2m-1f wien Än Am & A£v =1-2m 


< À +v entier positif ou négatif impair 





9 9 


2 4 


Elle n'est donc pas définie sur ces valeurs. C'est une des raisons pour lesquels on construit les 
fonctions de Lommel de deuxième espèce. 
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La fonction de Lommel de deuxième espèce est définie comme suit avec les deux formes 
équivalentes : 





ау ld Ё 
S, „(2) solution de l'équation inhomogène de Bessel ы К Sd № = 74 
Z 





(ee 


2 2 | el Zr r, @)-cof л, o 


Sin(v л) 
а Avl Av Y ((A-v) (a-v) 
S,,(2)=5,,(2)+2 T S Г Sin жа: J,(z)-Cos Ta | Ү (2) 


Pour v eN mais aussi v eR 
2 


Mais cette représentation ne marche pas aux pôles des fonctions Gamma, soit lorsque l'écart entre 
les paramètres À et v est un nombre impair négatif (voir cas précédent), ce qui a pour conséquence 


on le verra par la suite d'être impropre pour le calcul d'intégrales de la fonction de Bessel: 


Eee EE 


UE VV быйыл =-(2k +1) 
2 2 

Impropre par exemple pour 

[dz Ja) [dz A 

[dz z 4,0) [dz z J,(2) 





zarf 
Say (2) = Say (z) + 

















Convenable pour 
[dz Jy) [dz Ais 
faz z J, (z) faz z J.(z) 
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Lorsque l'écart entre les paramètres À et v est un nombre impair négatif alors on utilise plutôt les 
formules suivantes qui permettent de construire la solution dans ce cas. 

А-у 2-Qk41)2 4 =v - (2k 41) 

k-1 (— 1)” 7 2k+2m (— 1)“ S (2) 


S = 
ve юрт ишы) 





к. гт 


21 тГ(у+т+1) 
(a) в Symbole de Pochhammer 
Lorsque v =0 


2 — S su (Z) 1 8 í a, 
On part de $, (z) = = a = 5 10(2) = 5 mi i -Sua 2} 


А=-1 


_ gân 4 +3 = k G) asp A3 ? À TA 
Spa) =2 (22 ŽD ans) +2 (222) ed Ze o-s Zx oC) 


Lë ER z SEN 1 
>S ()=< Y =Â =Z EE aV (т+1) л 


2, (т) 


Dans ce cas, il est donc possible de construire une fonction de Lommel valide et donc de l'utiliser 
dans la détermination de certaines intégrale de Bessel. 





Les formules de récurrence sur les fonctions de Lommel de première et deuxième espèces sont : 


AH 
s 81,5, (2) 


DEE ч] 4,0227 - ay -v J.) 


Dum S A, (z)= (4 ER 1), (8) 
SE =(4-v-1)5, аба) 
dz ДА, | 
z" _ Sia) PES 2 y 
2) = i des "Wé Ё É 
5, (z) quae) cen ) Guth 
Teu И (ЕУ, (z)= (4 +V- 1), 4, (2) 





dS, (z) ze 


єй у, (2) = (A -V- DS asd (z) 


2 
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Avec le développement en série hypergéométrique des fonctions de Lommel et Bessel, on obtient : 
su 4-V+1\( 4+vV+1 
2) T д 


2 2 Lues) y.o- od KZ, ku 


Sin(v л) 


1 Ca z? 1 zy z? 
от) |у) WH X |) 


y 3+1-v 3+4+v 7° 
(cy | Fa)" 





SÉ (2) = Hä (2) + 














S, (z) = 2 2 


4 
{А-у Äis @-у+) _ "n (ev). | | 
келесе ы 




















onm d PEN T" (4-v+1)_ E л Lond A+v+1 Si (4+v+1) _ R л 
2 2 ]1-A-v 2 2 1-A-v 
d 2 d 2 





Et Sale zf ul —5ї(ул)Г(1+ъ)= 














r(v) r(-v) 
3+2-v m z 
Sar (2) = [ы x } 
4 
y pan fit) 


+ 











Siva) ey 1227 däi "n eset o) t 


N 

















gum 3+A-v 3+А+у z? 
S; (z) = 1; у= 
OQ) Guy 21 ( 2 2 4 
fete) rfi tie) | 
2 z 2 
+ | | Р 
2 4 








(2 
2 


z^" | 3-A-v 3+4+v d 
2| 15 y^ T 








9 H 


2 2 4 


24-1 „у (47 +1 
2 














bo a 
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Hors du développement hypergéométrique, cela donne : 
eee 





2 


of ed) 
zi (аы 


3 2 | {su = " +1) x) J.,(2)- s (etr d o 


Sin(v л) 


| {= 80 +, v) 
2^ 2 2 2 
Sin(v л) (а (2) 
2 


à T 
S (z) = $av (z)+ 











= Sir (z) 5 s, „(2)+ 











= Sir (z) = Ss, „(2)+ 
































2 2 
ү (25 2 LH D ) 
2^ x 2 
= S,,(2)=5,,(2)+ Р e? 
| | Sin(v л) =” ==” 
2 2 
dl {= 
25$,0)75,,C) «2^ T()r(-v foie)” Eca (2) 
Г D 
2 2 
(^ r SC 
2 2 


< 


> $,,()2s,, C) 2 2^ vr(v)r(- 








On ne peut utiliser cette représentation de la fonction lorsque le degré v est entier du fait de la 
divergence de la fonction Gamma. On utilise plutót l'expression suivante hors du développement 
en série hypergéométrique : 


EE y (E s AE y. (z)- od BH v) ch ol 


2 
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Calcul des intégrales des fonctions de Bessel à l'aide des fonctions de Lommel 


Les formules pour calculer ces intégrales, à l'aide des fonctions de Lommel de deuxième espèce, 
voir l'ouvrage de A. Erdelyi. H. Bateman « HIGHER TRANSCENDENTAL FUNCTIONS VOL Il », page 
90, formule , sont : 


[dz z2*J o = (й+у As Ale, dell ele (eil 
[а 27У, (2) = + 1) Osina O-z as] 

de même 

[ez o = (й+у As eg, deisde eil 
[az 21, (2) =[(4+v-1)z tie, del Adel, dell 


Pour montrer ces formules , оп va se baser sur les calculs présentés dans l'ouvrage de références sur 
les fonctions de Bessel de G.N.Watson « A treatise on the theory of Bessel functions », Cambridge 
University Press 1944 en page 350, du chapitre 10, section 10.74 « Functions expressible in terms 
of Lommel's functions ». Mais avant on rappelle également que les fonctions de Bessel possède les 
propriétés de dérivation suivantes : 


ACUTE ao AC) NES 
Г = Е - PN 
ША a(z) SM 

3 EE 


A l'aide des formules de récurrence, calculons les expressions suivantes : 


de, ©) E д5, ©) +v Z Se (z)=z (0), ы ©) 4 js Ak 27 "(A +V- 1)у,_ „—1(2) 





Oz Oz Oz 


Zeil An), c) e _ RH sez (A-v -1)s, „а (2) 


LA 
= rie =z  A-1-(v-1)-ljs,, (2)=2 70 -v -1)s,, (z 
2 D e) 6734 -1- (v —1)—1), (о) 2 z "(A -v -1)s, ,, (2) 
z 


Partons maintenant des expressions : 
v 1-v v 1-v 
(1) olz J, (z)z бау (2)) Sch e J, (2) T 2 (z) EL zu (z) SÉ S1-Lv-1 (2) 
д2 д2 i Oz 
-2'J, д a) z 7,22 (A -v - Ms, Xd 


o(z J, (2)2""5, eil 

















= (1) = = ZJ 1(2)в, 11(2)+27, (2)(2 =v UR LS 
alz” Cat 2 Бу 2 002°, „(2 
(2) | = Sa D 2 2 а 1. (oh, 9 ll Gë ) 


--z^*J,(z)z"s, 04277, (2) 2 (Av 1), 3,02) 
ez" "Je 3:463) - 


Oz 


= (2) 





-zJ (z)s,, (2) 3J, (204+ 1), (2) 
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Au passage remarquons que dans l'ouvrage de G.N.Watson, il s'est glissé une petite erreur de calcul 
dans l'expression 2. D'autre part il est à remarquer que les expressions sont tout autant valables avec 
les fonctions de Lommel de premiére espéce que de deuxiéme espéce. 


Combinons maintenant les expressions (1) et (2) en éliminant un inconnu dans le systéme 
d'équation, il vient : 


T ките) 





=2/, 4 (Gig, (2)+2/, GU =V- Tot (z) 


alJ Dz s. ell 


> zJ, |(2)5,_ Lv- (2) = & 


270, s 0)) _ 





-zJ,(z KA =v -l)s; ,, (2) 

















(2) = —2/,(®)з, „(®) +27, GOYA +v -1)s, 1,102) 
„=. т озан ЕЕ ааа) иду, чу, „д 
4" J: 222 e, (H An =; 


=—2/,(2)5, (2) - zJ, (2(A +v -1A -v -1)s, ,, (2) 
д", (адз "в, (0) - (9v AE J, 02", ill 
Oz 
Comme s,,(z)-z^— (a -1y -у? Ж (2) 
a, 267 a0 s, 0-0 ev 17, 025,1) 


-- (e)s C) - 24, CR 1 ech, C) 








Oz 
=- (208 (а 1) ach, s c)- zr оа) и), о) 
alz J, (22's, (0 - (A * v -1 J, (25, 1,0) 
Oz 
=>2J (2)z^7 + zJ, GA Ech, (о) zr (2-1 и), ›„() 
2,987 a0 5, 0- uc IA". 
lä ev -DI (о), (2), (2), (2)]) 
Oz 








=—J, (zz^ 


< z*J (2) = 





Au passage comme les formules de récurrence sont identiques avec les fonctions de Lommel de 
deuxième ol il E 


apa Quis „> „382, (2) 
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Et de plus les formules sur les fonctions de Bessel de deuxième espèces sont également les mêmes, 
il vient donc aussi : 








zY (z) = KB +V - 1)Y, Pi (2) = а (2)5,, (2)) 
гҮ (2) = e(za +v -1)Y, (2)5 d (2)-7,.,(2)8,, 2)) 


Ce qui permet bien de définir les intégrales indéfinies suivantes : 
[dz (= (й+у 1), (а), del deis, )] 


[ez (= ETOO RAONO) 

[ezY.o = a +v Reie, (2) У, (од, „(2)] 

[dz z^Y,(2) = dv RS, (2) - У, (298,02) 

Concrètement l'utilisation de ces intégrales indéfinies dépend de l'intervalle d'intégration, de la 


présence d'une singularité du membre à droite des expressions dans l'intervalle d'intégration ou dans 
ses bornes, ainsi que de la valeur des paramètres 


Considération sur les intervalles d'intégration et les paramètres des intégrales des fonctions de 





Bessel en utilisation des fonctions de Lommel 
Soit donc maintenant à calculer les intégrales définies : 


x2 

Í dz z^ J (z) 

Xi 

Si l'intervalle d'intégration comporte l'origine z=0 alors on utilise la fonction de Lommel de 
première espèce car sa valeur est définie en z=0, et où la fonction Bessel Y n'est pas définie pour la 
valeur 0., et et à condition que l'écart des paramètres ne soit pas un entier impair : 


faz zJ (2) = Х(А +v -19,095,4,400 72,4605, 09]. si Atv #-(2k+1) 


Dans le саз où l'intervalle d'intégration ne comprend pas la valeur 0 sur ses bornes alors on peut 
indifféremment utiliser les fonctions de Lommel de première et deuxième espèce, avec les 
restrictions évoquées sur les paramètres : 


faz, (2) =[x](2 +v 19, Ов, dl tele, tellt si 4 +v #-(2k+1) 


x 


(dz zJ, (2) = Bd +v 10), OS 1,400)-J4,.0)S5,, Ga Va,v>0 


х 
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Utilisation des fonctions de Struve 


Les fonctions de Struve et modifiées de Struve sont solutions respectivement des équations de 
Bessel inhomogéne : 








2 2 
1 
H,,(2),K, (z) solutions de l'équation inhomogène de Bessel ` = E + К De |» 








d'w ld Á 
L (z), M (2) solutions de l'équation inhomogène modifiée de Bessel ` e [ i É № 


Les fonctions de Struve possèdent le développement en série suivant : 


2 
ES d + H: ty d 
2 2 


2k 
B 
am. Vt a 
L (2)=-ie Haah) 2 = 
P (a + Hu ty | 
2 2 


K,(2)=H,(2)-Y,(2) М,(2)=1,(2)-1,(2) 

Ces fonctions ont le méritent d'être souvent proposées dans les logiciels de calcul mathématique. 
Parmi les intégrales que l'on peut calculer avec les fonctions de Struve, nous avons par exemple en 
intégrales indéfinies: 


Гало) = хло) + | он) ләнд] [агл = хло) 








[dz лб) = 40)- лән, О) н, QT [42 49) = xr 09+ (4) 6] 
[dz z* J (z) -G* -9x* (x) + 39 J (x) + SLE Lon, (х) 4, 69H, ©] 

Í dz I(2) =x I, (x) + = BOLD- 0) (х)] Í dz z I,(2) = х1(х) 

[dz z? 12x a- = LOLO-L0L0)] [dz 1,07 xo? +4), (х)—2х1,(х)] 


[de 2 I (2) =(x*+9x?)IH,(x)-3x* 1,099 E Drei, 69 - 1,691,0)] 
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Ainsi que : 

Јата) = хт) + TEKO- 

Jd KG) = хк) RL 0 (0) д] [dk (= bus 

[dz Кб) = х KG)» 5 [K.L + KL] [dz KG) = den 4) 00 2x, 00] 


[dzz*K,(3)= -(x*+9x?)K,(x)-3x3 K,(x)+9 - [K, GOL, (x) - K (х) (х)] 


Les intégrales avec les puissances impaires de z ne dépendent pas des fonctions de Struve, tandis 
que celle avec les puissances paires peuvent s'exprimer à l'aide de ces fonctions. 


Calcul d'intégrales dans des problèmes aux limites inhomogéne de Robin (voir système 
paraboloidal de révolution) 


Dans certains problémes aux limites rencontrés plus tard dans ce document, il apparaít nécessaire 
de calculer des intégrales de la forme : 


( xV1+ x? ( 
dx—————— J,(ax) et |аххүі+х?Ј (ax 
| І+еу1+ x? S | S 


On peut calculer ces intégrales à l'aide du développement en série des fonctions de Bessel, comme 
suit : 


l=+œ0 C 1) vc 


J,(x)= X 


CDU +v +1)! 2"? 


b 1=+% _ Y v+21 b 
= (dx xvV1+ x*J, (ax) = >. Ci a [ax e? 
0 0 


ZG TU zw +1)! 27? 


xix [=+ю (- 1) а”?! à х x 


b 
NP e уш d 
| “Hd Jl bier] ET 


b b v22L | 2 
Posons Ly(b,v +21+1)= [dx 1+ x? 1 (b,c,v +21+1)= [а= Lex 
0 


9 еу" 














b 1=+00 ee? 1 у+21 
[ахі xJ, (ах) = > (ei og 1, (b,v +21+1) 
0 


ELE eva T2 
— 


b / 2 l=+0 En v42l 
fax XVI+ x ) (-1) a 


————J = 
о lteXlex “(a Zentre Oz 





1,(b,c,v+21+1) 
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Soit donc à calculer ces deux intégrales : 





Iy(b,u)- |dxx"N1-x^ 1„(Ь,с,и)= | dx £ 

(b. u) | y x b.c H) | кы 

b +1 

[ат =L r 1,231 #33, EM 
РЕЧІ 2571 3 














u+ 2. 2 2,2 
[a мМ х2 KT Jens! 1 Hla t3, es VEY N ньня Во) uuo | 
J TE (N-e) 5° 2 2 pes с Z” D pe j 


Lim Аррей = : PELA AAS: -b?; Z => F, EU Duc m 
2 2 2 1-c 2,420 2 























2 2, 
FL EE C jJ] ell? p c г) Lis 
2 2 1-с u+l 2 2 1-c 
Et Lim 5 =1-b? 2 L = 0 
c0 C u + c20 C 
No / 2 url 
= Lim |ахх“ brest (pe) 
ат 1+сл+х?° (u+1) 2 2 2 


: ANS Г(а +т+ п) T(b + m) CU + п) r(e) x” y” 
H A IIF (a,b ,b,;c; x; y ) = 1 2 
Ж С „е y) 2, per: T(a) Г(Ъ,) T(b,) T(c+m+n) m! n! 
Les deux cas c-1 et c--1 se déduisent facilement la premiére intégrale et d'une intégrale connue. 


En effet : 


b url 
1,(b,u)= [dx x" V1+x° = S as A M 
0 


àv d i 
^o x" l+ x ` тыз 
1„\Ь1,и)= | ах = NU -1 dx ————— 
„ (b,1, u) | те 4 E | - 
b b 
Lb a) 1, (5-1, )= [dx ar d Nee d 2 dx i ee 
0 0 


1+V1+x? Idi 


p 1 u-1 u+1 , 
--2I,lb,u-2)--2 x, F | d ;=b 
GE H ) 1-1 2 | 2' 2 2 


b b 
ban) fase) | mee 
0 X EZ X 0 



































и-1 Së и-1 ux 
BELI: b «n 1 uA uel ch Ki Е: | 
H-l u+l 


и-1 "m и-1 ul 
Ld b «A lu ] gil, gd p^ b 
"m SET 











H-l u+l 
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D'où le développement en série des deux intégrales : 


b 1=+00 S d у+21 
| &ху1+ xJ, (а) = у = ЕЛИ. 1, (b,v +21+1) 
0 


EST +v um 
b | 2 l=+0 y v«21 
| PES 1+х s (-1) a 


J (ax) L(b,c,v - 2141 
A pre Je LG 22 ENT 


p 1 u+l u+3 
1, (b.u)= E DE E 














u+l X 2: 1877 
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